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ABSTRACT 


Quantum  algorithms  and  complexity  for 
certain  continuous  and  related  discrete 

problems 


Marek  Kwas 


The  thesis  contains  an  analysis  of  two  computational  problems.  The  first  problem 
is  discrete  quantum  Boolean  summation.  This  problem  is  a  building  block  of  quantum 
algorithms  for  many  continuous  problems,  such  as  integration,  approximation,  differ¬ 
ential  equations  and  path  integration.  The  second  problem  is  continuous  multivariate 
Feynman-Kac  path  integration,  which  is  a  special  case  of  path  integration. 

The  quantum  Boolean  summation  problem  can  be  solved  by  the  quantum  sum¬ 
mation  (QS)  algorithm  of  Brassard,  Hpyer,  Mosca  and  Tapp,  which  approximates 
the  arithmetic  mean  of  a  Boolean  function.  We  improve  the  error  bound  of  Brassard 
et  al.  for  the  worst-probabilistic  setting.  Our  error  bound  is  sharp.  We  also  present 
new  sharp  error  bounds  in  the  average-probabilistic  and  worst-average  settings.  Our 
average-probabilistic  error  bounds  prove  the  optimality  of  the  QS  algorithm  for  a 
certain  choice  of  its  parameters.  The  study  of  the  worst-average  error  shows  that 
the  QS  algorithm  is  not  optimal  in  this  setting;  we  need  to  use  a  certain  number  of 
repetitions  to  regain  its  optimality. 

The  multivariate  Feynman-Kac  path  integration  problem  for  smooth  multivariate 
functions  suffers  from  the  provable  curse  of  dimensionality  in  the  worst-case  deter¬ 
ministic  setting,  i.e.,  the  minimal  number  of  function  evaluations  needed  to  compute 
an  approximation  depends  exponentially  on  the  number  of  variables.  We  show  that  in 
both  the  randomized  and  quantum  settings  the  curse  of  dimensionality  is  vanquished, 


i.e. ,  the  minimal  number  of  function  evaluations  and/or  quantum  queries  required  to 
compute  an  approximation  depends  only  polynomially  on  the  reciprocal  of  the  desired 
accuracy  and  has  a  bound  independent  of  the  number  of  variables.  The  exponents 
of  these  polynomials  are  2  in  the  randomized  setting  and  1  in  the  quantum  setting. 
These  exponents  can  be  lowered  at  the  expense  of  the  dependence  on  the  number  of 
variables.  Hence,  the  quantum  setting  yields  exponential  speedup  over  the  worst-case 
deterministic  setting,  and  quadratic  speedup  over  the  randomized  setting. 
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Chapter  1 


Introduction 


This  thesis  contains  an  analysis  of  two  computational  problems.  One  of  them,  Boolean 
summation,  is  discrete  and  the  other,  multivariate  Feynman-Kac  path  integration,  is 
continuous.  As  we  shall  see,  they  are  closely  related  if  we  study  them  in  the  quan¬ 
tum  model  of  computation.  Quantum  Boolean  summation  is  a  building  block  of 
quantum  algorithms  for  many  continuous  problems,  such  as  integration,  approxima¬ 
tion,  differential  equations  and  path  integration.  The  quantum  Boolean  summation 
problem  was  previously  studied  in  the  worst-probabilistic  setting.  We  improve  the  ex¬ 
isting  results  and  extend  the  analysis  to  two  more  settings — average-probabilistic  and 
worst-average.  This  extension  is  especially  important  for  the  analysis  of  multivariate 
Feynman-Kac  path  integration,  in  which  the  concept  of  randomized  queries  is  used. 
The  thorough  knowledge  of  quantum  Boolean  summation  properties  is  necessary  to 
establish  sharp  complexity  bounds  in  the  quantum  setting  with  randomized  queries. 

This  chapter  is  organized  as  follows.  Section  11.11  gives  an  overview  of  research  in 
the  area  of  quantum  computation  for  continuous  problems.  The  quantum  model  of 
computation  for  continuous  problems  with  deterministic  and  randomized  queries  is 
presented  in  Section  11.1.11  The  main  results  of  this  thesis  are  outlined  in  Section  11.21 
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CHAPTER  1.  INTRODUCTION 


1.1  Background 

Feynman  H3  was  the  first  to  suggest  computational  devices  based  on  the  quantum 
mechanical  principles.  He  conjectured  that  the  intrinsic  difficulties  of  simulating 
quantum  phenomena  on  classical  computers  might  be  overcome  by  using  quantum 
computational  devices.  Indeed,  the  computational  intractability  of  some  quantum 
phenomena  caused  by  exponential  requirements  of  time  and  memory  in  the  num¬ 
ber  of  simulated  components  makes  classical  computers  useless  for  simulating  such 
phenomena. 

Deutsch  [8]  first  studied  rigorous  models  of  quantum  computation.  Initially,  ad¬ 
vantages  of  quantum  computation  over  classical  computation  were  shown  for  some 
discrete  problems,  e.g.,  the  Deutsch  and  Deutsch- Jozsa  problems  of  distinguishing 
between  constant  and  balanced  Boolean  functions  [H  '9]  and  the  Simon  problem  of 
checking  Boolean  function  invariance  [43]. 

Shor  and  Grover  made  significant  contributions  that  initiated  the  explosion  of 
research  in  the  area  of  quantum  computation.  Shor  [42J  found  a  quantum  algorithm 
that  solves  the  problem  of  factorizing  a  composite  n-bit  integer  with  cost  of  order 
n 2  log  n  log  log  n,  which  provided  an  exponential  speedup  over  the  best  classical  algo¬ 
rithm  known.  Grover  [IS]  discovered  a  quantum  algorithm  that  solves  the  problem  of 
searching  an  unstructured  database  consisting  of  N  elements  with  cost  of  order  a /N, 
which  yields  a  quadratic  speedup  over  any  classical  algorithm. 

Novak  and  Heinrich  initiated  the  study  of  continuous  problems  in  the  quantum 
setting  in  the  framework  of  information-based  complexity.  Novak  [33]  studied  the 
quantum  complexity  of  the  multivariate  integration  problem  for  Holder  classes  and 
proved  exponential  speedup  over  the  worst-case  deterministic  setting  and  roughly 
quadratic  speedup  over  the  randomized  setting.  Heinrich  [15]  extended  the  quantum 
computation  model  to  continuous  problems  and  dealt  with  integration  in  Lp  spaces, 
showing  results  similar  to  [33] .  A  selection  of  continuous  problems  has  been  considered 
since  then.  A  partial  list  includes: 
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•  approximation  for  Korobov  spaces  [35]. 

•  approximation  for  Sobolev  spaces  mm, 

•  the  eigenvalue  problem  mm  and  the  Sturm-Liouville  eigenvalue  problem 

•  integration  for  Sobolev  spaces  unnmEi], 

•  ordinary  differential  equations  [23]. 

•  parametric  integration  [5T], 

•  path  integration  [4B!|- 

The  quantum  setting  gives  an  exponential  speedup  over  the  worst-case  deterministic 
setting,  and  a  quadratic  speedup  over  the  randomized  setting  for  integration  over 
Sobolev  spaces  of  multivariate  functions.  The  same  is  true  for  path  integration. 
For  parametric  integration,  quantum  algorithms  turn  out  to  be  much  faster  than 
the  classical  ones  for  some  choices  of  the  problem  parameters.  However,  there  is  no 
improvement  for  some  other  choices  of  these  parameters.  A  similar  situation  occurs  for 
approximation  of  functions.  In  particular,  for  Korobov  spaces  we  obtain  polynomial 
speedups  over  the  known  algorithms  for  the  classical  settings.  For  Sobolev  spaces  the 
complexity  in  the  classical  and  quantum  settings  is  known.  For  some  choices  of  the 
problem  parameters  we  have  polynomial  speedup;  for  some  other  choices  there  is  no 
speedup.  For  ordinary  differential  equations  we  again  have  polynomial  speedup. 

The  power  of  the  quantum  setting  for  the  eigenvalue  problem  and  the  Sturm- 
Liouville  eigenvalue  problem  depends  on  what  kind  of  quantum  queries  are  permit¬ 
ted.  For  bit  queries,  which  are  used  for  the  other  problems  mentioned  in  the  previous 
paragraph,  we  have  at  most  polynomial  speedups  over  the  classical  randomized  set¬ 
ting.  However,  if  power  queries,  see  [3J  (37,  38],  are  permitted  we  have  exponential 
speedup  over  even  the  randomized  setting  for  the  eigenvalue  problem.  We  recall  that 
power  queries  are  used  in  the  phase  estimation  algorithm.  In  this  thesis  we  deal  with 
several  variations  of  bit  queries  and  we  do  not  consider  power  queries  at  all. 
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CHAPTER  1.  INTRODUCTION 


1.1.1  Quantum  computation  model 

We  now  give  a  brief  overview  of  a  simplified  quantum  model  of  computation  for  con¬ 
tinuous  problems  with  deterministic  and  randomized  queries.  Bit  queries  for  discrete 
problems  are  also  outlined  in  Section  12.21  see  |i30j .  The  model  with  deterministic 
queries  for  continuous  problems  is  thoroughly  described  in  [T5].  Randomized  queries 
are  studied  in  [52]  •  We  shall  use  the  framework  outlined  below  in  the  remaining  part 
of  this  chapter,  as  well  as  in  Chapter  [3]  where  we  study  multivariate  Feynman-Kac 
path  integration. 

We  start  with  a  general  computational  problem  formulation.  For  a  given  class  F 
of  input  functions  /  :  D  — »  C  we  want  to  approximate  the  solution  operator 

S  :  F^G, 

with  G  being  a  norrned  space  whose  norm  is  denoted  by  ||  •  ||g.  We  approximate  S(f) 
by  a  quantum  algorithm  as  described  below. 

First,  we  transform  a  given  input  function  /  G  F  by  using  a  classical  algorithm  Ps 
with  s  classical  function  evaluations  and  obtain 


/  —  Ps(f)  ■  D  —>  C. 


The  goal  of  this  first  classical  step  is  to  prepare  an  input  for  actual  quantum  com¬ 
putation.  This  can  be  used,  for  instance,  to  achieve  variance  reduction  of  the  input 
function  /  through  approximation,  as  it  is  done  for  multivariate  Feynman-Kac  path 
integration  in  Section  13.5.31  Afterwards,  we  use  the  transformed  function  /  as  the 
input  to  a  quantum  algorithm. 

We  recall  that  the  main  part  of  a  quantum  algorithm  is  a  sequence  of  unitary 
operators  acting  on  the  space  Tik  =  C2  ®  ®  C2,  which  is  a  tensor  product  of  k 
copies  of  the  two  dimensional  complex  Hilbert  space  C2,  i.e. , 


Un(D  =QnQfQn-l-"QlQjQo. 


(1.1) 
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Here,  Qo,  Qi,  ■  ■  ■  ,Qn  are  unitary  operators  and  Qj  is  a  quantum  query  for  /  G  PS(F). 
Quantum  queries  are  used  to  collect  information  about  an  input  function  /  and  play 
a  role  analogous  to  the  use  of  function  values  in  the  worst  case  and  randomized 
settings.  In  this  thesis  we  use  the  most  commonly  studied  bit  quantum  queries  in  the 
deterministic  and  randomized  forms. 

For  a  Boolean  function  g  :  {0, 1, . . . ,  2k~1  —  1}  —>  {0,1},  the  bit  query  Qg  is 
defined  as 

Qg\j)\y)  =  \  j)\y@g(J))- 

Here,  | j)  €  TLk-i  and  | y)  G  7di,  with  ©  denoting  the  addition  modulo  2. 

For  a  real  function  h  :  [0, 1]  — >  [0, 1]  the  bit  query  is  constructed  by  taking  the 
most  significant  bits  of  the  function  h  evaluated  at  some  points  tj.  More  precisely,  as 
in  [15],  the  bit  query  Qh  is  of  the  form 

Qh\j)\y)  =  \j)\y®  p(h(T{j)))), 

where  the  number  of  qubits  is  now  k  =  m!  +  m"  and  |  j)  G  Ttrn/  and  | y)  G  Ttm"  with 
some  functions  f3  :  [0, 1]  — ►  {0, 1, ... ,  2m"  —  1}  and  r  :  {0, 1, ... ,  2m'  —  1}  — ►  [0, 1]. 
Hence  we  compute  h  at  the  points  tj  =  r(j)  G  [0, 1],  and  then  take  /3(h(tj ))  which  is 
the  m"  most  significant  bits  of  h(tj)  .  The  randomized  quantum  query  Qh  =  Qh, u, 
defined  as  in  [52],  depends  on  a  random  element  u  G  H,  which  indicates  that  we 
compute  the  values  of  h  at  randomized  points  t]W. 

Therefore,  if  we  use  randomized  queries,  the  unitary  operator  (11.111  depends  on  a 
random  element  u  and  is  of  the  form 


Un,u)\f )  Qn  Q Qn—1  '  '  '  QlQf,u>  Q 0 

with  n  being  the  number  of  randomized  queries.  As  usual,  we  assume  that  the  initial 
state  is  |0)  and  we  compute  the  final  state 


W’f)  =  Un(f)  |0)  =  QnQfQn- 1  '  '  '  QlQfQo\0) 
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for  deterministic  quantum  queries,  and 

|'0/,w)  =  Un}W(f)  |0)  =  QnQf,uiQn- 1  ’  '  ‘  QlQf,u>Qo\fy 

for  randomized  ones.  Then  we  perform  a  measurement  of  the  final  state  and  obtain 
an  outcome  j  e  {0, 1, . . . ,  2fc  —  1}  with  probability 

Pfti)  =  \{i’f\j)\2  or  Pf.uti)  =  \&f,Jj)\2- 

Knowing  the  outcome  j  we  compute  the  final  result  on  a  classical  computer  as  <f>(j) 
and  c p(j )  for  some  mappings  <f>,ip  :  {0,1,...,  2k  —  1}  — ■>  G.  Thus,  the  quantum 
algorithm  An  yields 

KUO)  =  4>{j)  or  An,u>(f,j)  = 

The  most  commonly  used  error  criterion  for  quantum  algorithms  with  determin¬ 
istic  queries  is  the  worst-probabilistic  error.  This  is  defined  for  p  G  (|,  1]  by 

eWOT~pro(An,Ps,  S,p )  =  sup  min  max  || S(f)  -  An(Ps(f),j) ||G,  (1.2) 

f£F  J:/i(J,f)>p  j&J 

with  /i( J,  /)  =  J2j£jPfU)-  We  consider  here  the  probabilistic  error  with  respect 
to  the  quantum  algorithm  outcomes.  In  words,  for  a  fixed  /  G  F,  we  take  the 
maximal  error  for  the  best  outcomes  with  probability  at  least  p.  Then  we  maximize 
the  probabilistic  errors  over  all  input  functions  from  F.  In  addition  to  (1 1 . 2 p .  we  shall 
consider  a  wider  selection  of  error  criteria  for  quantum  algorithms  with  deterministic 
queries  in  Chapter  [2]  for  the  Boolean  summation  problem. 

For  quantum  algorithms  with  randomized  queries  we  use  the  worst-average-average 
error  criterion 

,  2fc-l  \  1/2 

=“^’‘(1  n,P„  S)  =  sup  (  E  V  p/,„(j)||S(/)  -  Av(P,(!),j)  11/  ,  (1.3) 

p;  / 

where  E  is  the  expectation  over  the  probability  space  hi.  For  a  fixed  input  function  /  6 
F  we  measure  the  L2  average  error  over  the  probability  space  hi  and  with  respect  to 
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the  quantum  algorithm  outcomes.  Then  we  maximize  the  average  errors  over  all  input 
functions  from  F.  Observe,  that  if  we  neglect  the  expectation  E  in  (11.3)1.  then  the 
sum  denotes  the  error  of  the  quantum  algorithm  with  respect  to  the  probability  of 
all  the  outcomes,  i.e. ,  we  consider  the  error  when  the  algorithm  succeeds  and  when  it 
fails  according  to  (11. 2ft .  This  error  criterion  is  a  reasonable  choice  for  problems  where 
the  result  of  an  algorithm  cannot  be  verified.  Besides,  it  is  a  stronger  error  criterion 
than  (11.21),  and  therefore  by  Chebyshev’s  inequality  the  results  obtained  for  (11. 3|)  can 
be  extended  to  (11.2)1 . 

As  usual  we  are  also  interested  in  the  query  complexities  for  the  respective  error 
settings,  i.e.,  the  minimal  numbers  of  quantum  queries  (deterministic  or  randomized) 
and  the  minimal  number  of  the  classical  function  evaluations  that  are  needed  to 
guarantee  that  the  error  does  not  exceed  e,  i.e., 

nwo F)  =  min{s  +  n  :  3Ps3An  ewor-pro(An,  Ps)  <  e}, 

nwor-avg-avg(£)  pj  =  min|s  +  R.  3Pg3An  e  wor-avg-avg  ^  Ps)  <  e} . 

The  quantum  setting  with  randomized  queries  is  considered  in  Chapter  [3]  for  the 
multivariate  Feynman-Kac  path  integration  problem. 

1.2  Research  results 

We  now  outline  the  results  of  this  thesis,  which  are  proved  in  Chapters  [2]  and  [3l  The 
results  of  Section  11.2.11  are  published  in  [19,  [3(1].  The  results  of  Section  11.2.21  are 
published  in  [29]  and  included  in  [27J  [28] . 

1.2.1  Quantum  Boolean  summation  in  various  settings 

For  all  the  problems  mentioned  in  Section  11.11  except  the  eigenvalue  problems,  the 
results  are  based  on  the  quantum  Boolean  summation  (QS)  algorithm  of  Brassard, 
Hpyer,  Mosca  and  Tapp  [4].  The  optimality  of  this  algorithm  was  proved  by  Nayak 
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and  Wu  ED,  using  the  polynomial  method  for  quantum  query  lower  bounds  obtained 
by  Beals,  Buhrman,  Cleve,  Mosca  and  de  Wolf  PJ.  Because  of  the  significance  of  the 
QS  algorithm  for  the  quantum  complexity  of  continuous  problems,  it  is  important  to 
study  this  algorithm  in  greater  detail.  The  performance  of  a  quantum  algorithm  can 
be  studied  in  various  error  settings,  depending  on  how  the  input  and  output  of  the 
algorithms  are  treated.  We  may  have  worst  and  average  behavior  with  respect  to  the 
input,  and  probabilistic  and  average  with  respect  to  the  output.  The  QS  algorithm 
in  the  worst-probabilistic  error  setting  was  studied  in  pE].  It  was  shown  that  its  error 
is  of  order  1/M  with  M  queries. 

To  define  the  Boolean  summation  problem  in  the  framework  of  Section  11.1.11  we 
set  D  =  {0, 1, . . . ,  N  —  1}  and  C  =  {0, 1}.  The  class  F  of  input  functions  is  just 
the  class  Bat  of  all  Boolean  functions  defined  on  an  N  element  set.  The  solution 
operator  S  computes  the  mean  of  a  Boolean  function 

i  N~ 1 

S(/)  =  ^£/(i)  V/eF  =  B„.  (1.4) 

i= 0 

The  QS  algorithm  of  Brassard,  Hpyer,  Mosca  and  Tapp  [4]  computes  an  approxi¬ 
mation  to  the  solution  operator  S.  Information  regarding  the  Boolean  function  is  sup¬ 
plied  by  quantum  queries.  Suppose  that  we  use  M— 1  quantum  queries.  Obviously,  the 
only  case  of  interest  is  when  M  is  much  smaller  than  N.  Brassard  et  al.  [4j  proved  that 
the  worst-probabilistic  error  of  the  QS  algorithm  is  at  most  7r/M  +  7r2/M2,  with  prob¬ 
ability  8/ 7r2  =  0.81 ....  Nayak  and  Wu  [31J  showed  that  with  probability  p  e  (|,  1) 
the  error  of  any  quantum  algorithm  that  uses  no  more  than  M  —  1  quantum  queries 
is  bounded  from  below  by  the  quantity  proportional  to  M-1.  Therefore,  the  QS 
algorithm  enjoys  the  smallest  possible  error  modulo  a  factor  multiplying  M-1.  Since 
the  QS  algorithm  has  so  many  applications,  as  mentioned  also  in  Section  11.11  it  is 
important  to  check  whether  the  error  estimate  of  Brassard  et  al.  is  sharp  and  how  the 
error  decreases  if  we  lower  the  probability  p  =  8 / 7r2  to  p  >  It  also  seems  reasonable 
to  study  whether  the  QS  algorithm  retains  its  optimality  with  different  error  criteria. 
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We  recall  that  the  estimate  of  Brassard  et  al.  corresponds  to  the  worst-probabilistic 
setting,  which  is  most  frequently  used  in  the  study  of  quantum  algorithms.  We  also 
study  the  average  performance  of  the  QS  algorithm  with  respect  to  a  measure  on 
Boolean  functions.  This  is  the  average-probabilistic  setting.  In  the  worst-average 
setting,  we  study  the  worst  performance  with  respect  to  Boolean  functions  and  the 
average  performance  with  respect  to  all  outcomes  of  a  quantum  algorithm.  We  add 
in  passing  that  the  worst-average  setting  is  usually  used  for  the  study  of  the  classical 
Monte  Carlo  algorithm.  This  setting  is  also  used  to  establish  a  lower  bound  for  the 
quantum  summation  and  integration  algorithms  that  use  randomized  queries,  see  E23- 
We  study  error  bounds  in  the  worst-probabilistic  and  average-probabilistic  settings 
with  probabilities  p  G  (1/2,  8/ 7r2] .  If  we  want  to  obtain  error  bounds  with  higher 
probability,  it  is  known  that  it  is  enough  to  run  the  QS  algorithm  a  number  of  times 
and  take  the  median  as  the  final  result,  see  e.g.,  [15]. 


1.2. 1.1  Worst-probabilistic  error 

The  worst-probabilistic  error  of  the  QS  algorithm  is  defined  according  to  (11.21) 


3wor-pro 


(QS,  p)  —  sup  min  max \S(f)  -  QS(/,  j)\ 


f SBjv  J:n(J,f)>p  i&J 

with  S  given  by  (11.11)  and  p  G  (1/2, 1],  Note  that  /  =  /,  i.e. ,  we  do  not  use  a  classical 
computer  to  prepare  the  input  for  the  QS  algorithm. 

Our  results  improve  the  worst-probabilistic  error  bound  ewor"pro(QS,  8 / 7r2)  <  i r/M 
+  7 r2/M2  of  Brassard  et  al.  Namely,  the  worst-probabilistic  error  of  the  QS  algo¬ 
rithm  ewor-pro(QS,  8 / 7T2)  is  at  most 

3  7 r  3 

- -  with  probability  —  =  0.81 .... 


4  M 


7T 


Furthermore,  the  last  estimate  is  sharp  when  M  and  N/M  are  large.  Indeed,  for 
p  G  (1/2,  8/7 r2],  the  error  of  the  QS  algorithm  is  at  most 

(l  -u 1(p))  7T 


M 


with  probability  p, 
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where  v  1  is  the  inverse  of  the  function  w(A)  =  sin2(7rA)/(7rA)2.  We  have  1—v  l(p)  € 


(2,  §]  and  it  is  well  approximated  by  2p  +  In  particular,  for  the  most  frequently 
considered  values  of  p  we  have 


(1  -  =  1.75  . . . ,  (1  -  ))tt  =  2.23 

(l  —  V-1(8/7 T2))7T  =  |7T  =  2.35 _ 


The  detailed  analysis  is  given  in  Section  12.3.11 

1.2. 1.2  Average-probabilistic  error 

To  define  the  average-probabilistic  error  of  the  QS  algorithm  we  need  to  equip  the 
input  function  class  MN  with  a  probability  measure  p.  Then  for  p  e  (|,  1]  we  define 


p(/)  ;  min  max  \  S(f)  —  QS(/,  j)| 
J-pW)>p  jsJ 


/eBjv 


We  are  then  interested  in  the  average  error  that  holds  with  a  certain  fixed  probabil¬ 
ity  p. 

I11  the  average-probabilistic  setting,  we  shall  consider  two  measures  on  the  set 
of  Boolean  functions.  The  first  measure  px  is  uniform  on  Boolean  functions,  while 
the  second  measure  p2  is  uniform  on  arithmetic  means  of  Boolean  functions.  These 
measures  have  different  properties.  Although  the  mean  element  of  the  arithmetic 
means  is  |  for  both  measures,  the  first  central  moment  is  of  order  A"1/2  for  the  first 
measure  pi,  and  about  |  for  the  second  measure  p2.  The  first  central  moment  is 
exactly  equal  to  the  error  of  the  constant  algorithm  that  always  outputs  This 
explains  why  we  can  obtain  the  error  of  order  iW1/2  without  any  quantum  queries  for 
the  first  measure.  This  provides  the  motivation  for  us  to  check  whether  the  error  of 
the  QS  algorithm  enjoys  a  similar  property.  Our  investigation  shows  that  that  this  is 
indeed  the  case  iff  M  is  divisible  by  4.  More  precisely,  the  average-probabilistic  error 
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of  the  QS  algorithm  for  the  measure  pi  satisfies 


0(min{M  1,N  1//2})  for  M  divisible  by  4, 


0(M-1) 


for  M  not  divisible  by  4. 


(1.5) 


Moreover,  using  the  lower  bounds  of  [3HJ  when  M  is  divisible  by  4,  and  since  we 
have  eavg_pro(QS,  p,  pi)  =  for  M  not  divisible  by  4  and  large  N/M  the  esti¬ 

mates  (ll.5[l  are  essentially  sharp. 

For  the  measure  P2,  the  average-probabilistic  error  of  the  QS  algorithm  is  of 
order  M^1  for  all  M.  For  both  measures,  the  upper  bounds  mentioned  earlier  match 
lower  bounds  that  were  obtained  by  Papageorgiou  [36] .  Hence,  the  QS  algorithm 
enjoys  minimal  error  bounds  also  in  the  average-probabilistic  setting  if  we  choose  M 
divisible  by  4  for  the  first  measure  and  with  no  restriction  on  M  for  the  second 
measure. 

The  detailed  analysis  is  given  in  Section  12.3.21 
1.2. 1.3  Worst-average  error 

In  this  setting,  we  take  the  worst  case  performance  over  all  Boolean  functions  and 
the  average  performance  over  all  outcomes  of  the  QS  algorithm,  so  that 


wor-avg 


The  average  performance  is  measured  in  the  Lq  norm,  q  E  [l,oo].  Since  we  do  not 


use  a  classical  computer  to  prepare  the  input  for  the  QS  algorithm  we  have  /  = 
/  and  Pf(j)  =  Pf(j)-  This  setting  is  analogous  to  the  randomized  (Monte  Carlo) 


setting  used  for  algorithms  on  a  classical  computer.  The  worst-average  setting  seems 
natural  for  the  analysis  of  quantum  algorithms  for  the  same  reasons  that  motivates  the 
worst-average-average  error  criterion  of  (11 .3j).  The  results  depend  on  the  choice  of  q. 
Obviously,  for  larger  q,  the  effect  of  the  average  behavior  becomes  less  significant. 
In  fact,  the  limiting  case,  q  —  oo,  leads  to  the  deterministic  case  (modulo  sets  of 
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measure  zero).  Not  surprisingly,  for  q  =  oo,  the  results  are  negative,  i.e. ,  the  error  is 
of  constant  order. 

We  shall  study  error  bounds  for  large  M.  Without  loss  of  generality  we  as¬ 
sume  N  >  M,  so  that  if  M  tends  to  infinity,  then  so  does  N.  To  make  error  bounds 
independent  of  N,  we  take  the  supremum  over  N  >  M  in  the  corresponding  defini¬ 
tions  of  the  errors.  When  we  speak  about  the  sharpness  of  error  bounds,  we  usually 
take  a  large  M  and  select  a  still  larger  N  and  a  Boolean  function  for  which  the 
presented  error  bound  is  sharp. 

We  now  indicate  our  results  for  q  G  [l,oo).  The  worst-average  error  e^or_avg(M) 
of  the  QS  algorithm  with  M  quantum  queries  satisfies: 


wor-avg 


(QS) 


{@(ln  M/M)  for  q  —  1, 
0(M-1/9)  for  q  G  (l,oo). 


Comparing  with  [36],  we  conclude  that  these  error  bounds  do  not  match  the  com¬ 
plexity  lower  bounds.  We  shall  consider  a  slight  modification  of  the  QS  algorithm  to 
address  this  issue. 


Quantum  Boolean  summation  with  repetitions 

The  error  bounds  of  the  QS  algorithm  can  be  improved  by  the  use  of  repetitions. 
Namely,  we  repeat  the  QS  algorithm  2n  + 1  times  and  take  the  median  of  the  outputs 
obtained  as  the  final  output.  This  procedure  boosts  the  success  probability  of  the 
approximation  at  the  expense  of  the  number  of  quantum  queries.  We  show  that 
with  n  independent  of  M  and  depending  linearly  on  q,  we  decrease  the  QS  algorithm 
error  to  be  of  order  M _1.  Hence,  the  use  of  repetitions  is  particularly  essential  for 
large  q,  since  we  change  the  error  bound  0(M~1//<?)  without  repetitions  to  the  error 
bound  0(M_1)  with  repetitions. 

The  error  bound  of  order  is  optimal.  This  follows  from  the  use  of,  for 

instance,  Chebyshev’s  inequality  and  the  fact  that  the  lower  bound  is  sharp 
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in  the  worst-probabilistic  setting,  see  [36].  Hence,  the  QS  algorithm  with  repetitions 
is  also  optimal  in  the  worst-average  setting. 

The  detailed  analysis  of  the  QS  algorithm  in  the  worst-average  setting  including 
repetitions  is  given  in  Section  12.3.31 

1.2. 1.4  Quantum  Boolean  summation  simulation 

Simulating  arbitrary  quantum  algorithms  on  a  classical  computer  is  very  difficult  due 
to  the  exponential  time  and  memory  requirements  on  the  number  of  qubits.  Hence, 
such  simulations  can  be  run  only  for  input  data  of  moderate  size.  We  designed  a  mat- 
lab  procedure  simulating  the  QS  algorithm.  This  simulation  computes  amplitudes  of 
the  QS  algorithm  final  state  depending  on  the  Boolean  input  function  and  a  desired 
accuracy  e.  The  results  of  the  simulation  can  be  used,  for  instance,  to  visualize  the 
final  state  distribution.  Detailed  studies  of  the  QS  algorithm  show  that  this  task  can 
be  accomplished  with  cost  of  order  (N/e)  log2(l/e).  Here  N  stands  for  the  cardinality 
of  the  domain  of  a  Boolean  function.  We  recall  that  the  QS  algorithm  requires  of 
order  log2  N  +  log2(l/£)  qubits  so  the  memory  requirement  for  the  simulation  is  N/e. 
Thus  the  cost  of  our  simulation  is  optimal  modulo  a  logarithmic  factor  log2(l/e).  The 
presence  of  this  factor  is  due  to  the  use  of  the  Fast  Fourier  Transform.  Since  we  do 
not  know  any  faster  algorithm  for  the  discrete  Fourier  transform,  which  is  a  part  for 
the  QS  algorithm,  this  simulation  achieves  the  best  possible  cost  order  using  existing 
computational  tools.  The  details  are  given  in  Section  12.41 

1.2.2  Multivariate  Feynman-Kac  path  integration 

Path  integration  is  one  of  the  continuous  problems  that  has  been  studied  studied 
in  the  quantum  setting.  Traub  and  Wozniakowski  [48]  studied  general  path  inte¬ 
grals  with  smooth  integrands.  We  consider  a  special  case  of  path  integrals,  namely, 
multivariate  Feynman-Kac  path  integrals.  These  path  integrals  arise  in  many  prob¬ 
lems  of  quantum  physics,  quantum  chemistry,  and  even  in  financial  mathematics.  In 
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particular,  the  multivariate  Feynman-Kac  path  integral  expresses  the  solution  of  the 
initial  value  problem  for  the  diffusion  equation.  Studies  of  this  problem  in  the  worst- 
case  deterministic  setting  are  in  [29,  39].  We  study  multivariate  Feynman-Kac  path 
integration  in  the  randomized  and  quantum  settings. 

Path  integrals  are  defined  as  integrals  over  an  infinite  dimensional  space  equipped 
with  a  probability  measure.  A  path  integral  is  called  a  Wiener  integral  if  the  respective 
measure  is  the  Wiener  measure  w  on  the  space  C  of  continuous  functions  from  M+ 
to  Md.  In  the  case  d  >  2,  we  add  the  term  “multivariate”.  The  multivariate  Feynman- 
Kac  path  integral, 


(1.6) 


is  the  solution  of  the  initial  value  problem  for  the  multivariate  diffusion  equation 


We  assume  that  the  input  functions,  initial  value  v  and  potential  V  belong  to  a  certain 
function  class  F  for  which  the  integral  (11 ,6j)  is  well  defined. 

Many  algorithms  have  been  developed  for  the  univariate  case,  d  —  1,  and  the 
Feynman-Kac  path  integral  is  the  solution  of  the  diffusion  equation  with  one  space 
variable.  Most  of  these  algorithms  are  randomized.  They  are  obtained  as  follows. 
First,  the  path  integral  is  approximated  by  a  multivariate  integral  over  Mn,  with 
large  n,  and  then  this  integral  is  approximately  computed  by  using  a  randomized 
algorithm  such  as  Monte  Carlo. 

A  new  approach  was  proposed  in  [39]  for  the  univariate  case  d  —  1.  A  deterministic 
algorithm  based  on  Z/2-approximation  of  v  and  V  was  constructed.  This  approach 
was  modified  and  generalized  for  the  multivariate  case  in  [261129],  and  is  also  the  basis 
for  the  randomized  and  quantum  algorithms  presented  in  Chapter  [3J 

We  want  to  check  how  much  the  power  of  randomization  and  quantum  computa¬ 
tion  helps  in  solving  the  multivariate  Feynman-Kac  path  integration  problem.  This 
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question  has  been  addressed  for  general  path  integrals  with  smooth  integrands  in  [08j , 
where  it  is  shown  that  this  problem  can  be  solved  on  quantum  computers  exponen¬ 
tially  faster  than  in  the  worst-case  deterministic  setting  and  roughly  quadratically 
faster  than  in  the  randomized  setting.  One  of  our  questions  is  how  the  special  form  of 
the  Feynman-Kac  path  integrals  can  be  exploited  and  how  it  can  improve  the  general 
results  of  [08]. 

1.2. 2.1  Worst-case  deterministic  setting 

Although  we  focus  mainly  on  the  quantum  and  randomized  settings,  we  first  formu¬ 
late  the  computational  problem  of  finding  a  local  approximate  solution  of  the  path 
integral  (13.31)  at  a  given  point  in  the  worst-case  deterministic  setting.  Assume  that 
we  have  a  deterministic  algorithm  An  that  uses  n  function  evaluations  of  v  and  V. 
We  define  the  worst-case  deterministic  error  of  the  algorithm  An  as 

ewor(An,F)  =  sup  \zvy{u,t)  -  An(v,  V)\, 

v,V&F 

where  z(u,  t)  in  the  exact  solution  of  the  diffusion  equation  at  a  given  fixed  point 
(u,  t)  €  !d  x  (0,  oo). 

In  [26j  [29j  we  extended  the  algorithm  of  [39]  to  the  multivariate  case.  We  pre¬ 
served  its  original  structure,  but  we  used  uniform  approximation  instead  of  L2- 
approximation.  This  permits  the  algorithm  to  be  used  for  arbitrary  d,  while  an 
algorithm  based  on  ./^-approximation,  such  as  that  of  [39],  can  be  only  used  for  d  —  1. 
The  cost  of  computing  an  ^-approximation  is  roughly  of  order  e~a^U  for  a  certain  pos¬ 
itive  ol(F)  which  depends  on  a  given  class  F  of  input  functions.  The  upper  bound  is 
derived  from  the  complexity  of  uniform  approximation  for  the  class  F.  When  F  is  a 
class  of  r  times  continuously  differentiable  d-variate  functions,  we  have  ot{F)  =  d/r. 
So  in  the  worst  case  deterministic  setting  the  problem  suffers  from  the  curse  of  dimen¬ 
sionality.  Moreover,  this  algorithm  and  the  one  in  |J9j  require  the  precomputation  of 
a  large  number  of  coefficients  which  is  quite  difficult. 
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The  (information)  complexity  of  multivariate  Feynman-Kac  path  integration  in 
the  worst-case  setting  is  defined  by 

nwor(e,  F)  =  min{n  :  3  An  such  that  ewor(An)  <  e}. 

The  complexity  is  bounded  from  below  by  complexity  of  multivariate  integration,  and 
from  above  by  complexity  of  uniform  approximation.  Furthermore,  the  algorithm  of 
[2D]  is  almost  optimal  for  classes  F  for  which  the  complexities  of  the  multivariate 
integration  and  the  uniform  approximation  are  of  the  same  order.  This  holds,  for 
instance,  for  the  class  F  of  r  times  continuously  differentiable  multivariate  functions. 

1.2. 2. 2  Randomized  and  quantum  settings 

A  randomized  algorithm  An  depends  on  a  random  element  u  chosen  from  some  proba¬ 
bility  space  12,  with  n  denoting  the  number  of  function  evaluations.  In  the  randomized 
setting,  the  computed  approximation  of  the  exact  solution  z(u,  t)  is  then  a  random 
variable  An(v,  V‘,u)  that  depends  on  a  random  element  uj  e  12.  We  measure  the  error 
of  the  algorithm  An  with  respect  to  the  L2  norm,  so  that 

erand(An)  =  sup  (E(z„jv(u,  t)  —  An(v,  V ;  a;))2) 1/2  .  (1.7) 

v,veF 

We  are  also  interested  in  estimating  the  (information)  complexity  nrand(e,  F)  of  mul¬ 
tivariate  Feynman-Kac  path  integration,  i.e. ,  the  minimal  expected  number  of  func¬ 
tions  v  and  V  values  needed  to  compute  an  ^-approximation  in  the  randomized  setting, 
which  is  given  by 

nrand(e,F)  =  min{n  :  3  An  such  that  erand(An)  <  e}. 

We  now  turn  to  the  quantum  setting.  We  met  some  technical  difficulties  dealing 
with  deterministic  (bit)  queries,  so  we  left  this  case  as  one  of  the  open  problems  in 
Chapter  U  Our  results  are  only  established  for  randomized  queries. 

We  now  outline  the  results  obtained.  We  consider  quantum  algorithms  with  ran¬ 
domized  queries.  The  error  of  a  quantum  algorithm  An  is  defined  analogously  to 
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Section  11.1.11  as 

equant(An,  Ps)  =  sup  (EEq  \zVjV(u,  t)  —  An>ul(Ps(v,  V),j;  o;)|2)1/2  ,  (1.8) 

v,veF 

where  E  is  the  expectation  over  the  probability  space  fl,  and  Eq  is  the  expectation 
with  respect  to  the  distribution  of  the  quantum  algorithm  outcomes.  We  recall  that  Ps 
is  a  classical  algorithm  using  s  function  evaluations  that  prepares  the  input  for  the 
quantum  algorithm,  see  Section  11.1.11  Similarly  to  the  randomized  setting,  we  also 
want  to  know  the  minimal  number 

nquant(£,  =  min{s  +  n  .  3  Ps  3  An  such  that  equant(Hn,  Ps)  <  e}. 

of  randomized  quantum  queries  and  classical  function  evaluations,  which  are  needed 
to  guarantee  that  the  error  does  not  exceed  e. 

We  present  algorithms  that  compute  an  ^-approximation,  i.e.,  with  the  errors 
(11.71)  and  (jl.8|)  at  most  e,  and  provide  their  cost  analyses.  These  algorithms  are 
also  based  on  uniform  approximation.  However,  the  power  of  randomization  and 
quantum  computation  permits  us  to  improve  the  worst-case  deterministic  complexity 
bound  0[e~a ^).  The  number  of  function  evaluations  required  by  the  randomized 
algorithm  is  roughly  of  order  £~2aT)/(a(F)+2) ;  and  the  number  of  function  evaluations 
and  queries  required  by  the  quantum  algorithm  is  roughly  of  order  £~a(F)/(.a(F)+1) . 
We  stress  that  the  exponent  of  e^1  in  the  randomized  setting  is  at  most  2,  and  in  the 
quantum  setting  is  at  most  1. 

We  also  study  the  complexity  of  multivariate  Feynman-Kac  path  integration  in 
the  randomized  and  quantum  settings.  The  complexity  is  bounded  from  below  by  the 
complexity  of  multivariate  weighted  integration,  just  as  in  [3Sj.  The  upper  bounds 
are  provided  by  the  costs  of  the  algorithms  mentioned  above. 

For  the  class  F  of  r  times  continuously  differentiable  d-variate  functions  we  have 
ot(F)  =  d/r.  In  the  randomized  setting,  the  complexity  is  roughly  @(e_2/(1+2r'/rf)j-  jn 
the  quantum  setting  it  is  roughly  ©(£^1/b+r/rf)).  Furthermore,  we  know  algorithms 
that  use  0(e~2)  function  values  in  the  randomized  setting,  and  0{e~l)  function  values 
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in  the  quantum  setting,  with  the  factors  in  the  big  O  notation  independent  of  the 
number  of  variables  d.  In  both  cases,  the  curse  of  dimensionality  is  vanquished.  We 
thus  have  exponential  speedup  over  the  worst  case  setting.  For  d  r,  we  have 
quadratic  speedup  of  the  complexity  in  the  quantum  setting  over  the  complexity  in 
the  randomized  setting.  We  refer  the  reader  to  Chapter  |3]  for  a  detailed  discussion. 


Chapter  2 
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Chapter  2 

Quantum  Boolean  summation 

2.1  Introduction 

The  quantum  summation  (QS)  algorithm  (also  known  as  the  amplitude  estimation 
algorithm)  of  Brassard,  Hpyer,  Mosca  and  Tapp  [TJ  computes  an  approximation  to 
the  arithmetic  mean  of  all  values  of  a  Boolean  function  defined  on  a  set  of  IV  =  2n 
elements.  We  denote  the  class  of  such  functions  by  B/y.  Information  regarding  the 
Boolean  function  is  supplied  by  quantum  queries.  Quantum  queries  play  a  role  similar 
to  that  of  function  evaluations  in  the  classical  worst-case  and  randomized  settings. 
Suppose  that  we  use  M  —  1  quantum  queries.  The  only  case  of  interest  is  when  M  is 
much  smaller  than  N.  It  was  proved  in  [4]  that  the  error  of  the  QS  algorithm  is  at 
most 

Jj  +  Jp  With  Probabilit^  ^2  =  °-81  •  •  •  •  (2-1) 

Nayak  and  Wu  [31]  showed  that  for  any  p  e  (|,  1)  the  error  of  any  quantum 
algorithm  that  uses  no  more  than  M  —  1  quantum  queries  must  be  proportional 
to  M-1  with  probability  p.  Therefore,  the  QS  algorithm  enjoys  the  smallest  possible 
error  modulo  a  factor  multiplying  M-1. 

The  minimal  error  estimate  of  order  M~l  in  the  quantum  setting  should  be  com¬ 
pared  to  the  minimal  error  estimates  in  the  worst  case  and  randomized  settings  of 
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algorithms  using  M  —l  function  values.  It  is  known,  see  [33] ,  that  in  the  worst  case 
setting,  the  error  bound  is  roughly  |(1  —  M/N ).  This  means  that,  as  long  as  M  is 
much  less  than  N,  the  error  is  almost  and  is  therefore  0(M )  times  larger  than 
that  in  the  quantum  setting.  In  the  randomized  setting,  the  classical  Monte  Carlo  is 
almost  optimal,  and  the  error  bound  is  roughly  1/(2 y/M),  see  again  pi].  Hence,  the 
error  0(y/rM)  larger  than  that  in  the  quantum  setting. 

We  check  whether  the  estimate  (12.  ip  is  sharp  and  how  the  error  decreases  if  we 
decrease  the  probability  p  =  8 / 7r2  to  p  >  We  also  study  the  error  of  the  QS  algo¬ 
rithm  in  various  settings.  The  estimate  (12.11)  corresponds  to  the  worst-probabilistic 
setting,  which  is  most  frequently  used  for  quantum  algorithms.  The  average  perfor¬ 
mance  of  the  QS  algorithm  with  respect  to  a  measure  on  Boolean  functions  is  studied 
in  the  average-probabilistic  setting.  In  the  worst-average  setting,  we  study  the  worst 
performance  with  respect  to  Boolean  functions  and  the  average  performance  with 
respect  to  all  outcomes  of  a  quantum  algorithm.  This  setting  is  usually  used  for  the 
study  of  the  classical  Monte  Carlo  algorithm.  We  study  error  bounds  in  the  worst- 
and  average-probabilistic  settings  with  probabilities  p  G  (1/2,  8/7 r2].  If  we  want  to 
obtain  error  bounds  with  higher  probability  it  is  enough  to  run  the  QS  algorithm  a 
number  of  times  and  take  the  median  as  the  final  result,  see  e.g.,  mm- 

In  the  worst-probabilistic  setting,  we  show  that  (12. ip  can  be  slightly  improved. 
Namely,  the  error  of  the  QS  algorithm  is  at  most 

3  vr  .  ,  ,  ,  8 

-  —  with  probability  — . 

4  M  7 t2 


Furthermore,  for  large  M  and  N/M  we  prove  that  the  last  estimate  is  sharp.  I11 
particular,  for  p  G  (1/2,  8/7 r2]  we  prove  that  the  error  of  the  QS  algorithm  is  at  most 


(l-u  1(p))tt 


M 


with  probability  p, 


where  v  1  is  the  inverse  of  the  function  u(A)  =  sin2(7rA)/(7rA)2.  We  have  1—v  l{p)  G 
(2,  §]  and  it  is  well  approximated  by  -^ir2p+  For  the  frequently  considered  values 
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of  p  we  have 

(1  -  y-1(l+))7r  =  1.75  . . . ,  (1  -  tT^f  ))tt  =  2.23 

(l  —  V-1(8/7 T2))7T  =  |7T  =  2.35 _ 

In  the  average-probabilistic  setting,  we  consider  two  measures  on  the  set  of  Boolean 
functions.  The  first  measure  is  uniform  on  Boolean  functions,  while  the  second  mea¬ 
sure  is  uniform  on  arithmetic  means  of  Boolean  functions.  The  results  for  these  two 
measures  are  quite  different.  The  mean  element  is  |  for  both  measures.  However, 
the  first  moment  is  of  order  iV_1//2  for  the  first  measure,  and  about  |  for  the  second. 
The  first  moment  is  exactly  equal  to  the  error  of  the  constant  algorithm  that  always 
outputs  which  explains  why  we  can  obtain  the  error  of  order  IV-1/2  without  any 
quantum  queries  for  the  first  measure.  This  provides  the  motivation  for  us  to  check 
whether  the  error  of  the  QS  algorithm  enjoys  a  similar  property.  It  turns  out  that  this 
is  indeed  the  case  iff  M  is  divisible  by  4.  That  is,  for  M  divisible  by  4,  the  average- 
probabilistic  error  of  the  QS  algorithm  is  of  order  min{M-1,  IV-1/2},  and  if  M  is  not 
divisible  by  4,  then  the  error  is  of  order  M”1.  For  the  second  measure,  since  the  first 
moment  is  not  small,  the  average-probabilistic  error  of  the  QS  algorithm  is  of  order 
M-]  for  all  M.  For  both  measures,  the  upper  bounds  presented  in  this  paper  match 
lower  bounds  from  [36] .  Hence,  the  QS  algorithm  enjoys  minimal  error  bounds  also 
in  the  average-probabilistic  setting  if  we  choose  M  divisible  by  4  for  the  first  measure 
and  with  no  restriction  on  M  for  the  second  measure. 

In  the  worst-average  setting,  we  take  the  worst  case  performance  over  all  Boolean 
functions  and  the  average  performance  over  all  outcomes  of  the  QS  algorithm.  The 
average  performance  is  measured  in  the  Lq  norm,  where  q  €  [l,oo].  This  setting  is 
analogous  to  the  randomized  (Monte  Carlo)  setting  used  for  algorithms  on  a  classical 
computer.  Recall  that,  for  a  number  of  reasons  discussed  in  Sectionfl.2.1.31  the  worst- 
average  setting  is  a  natural  choice  for  the  analysis  of  quantum  algorithms.  As  we  shall 
see,  the  estimates  depend  on  the  choice  of  q.  Obviously,  for  larger  q,  the  effect  of  the 
average  behavior  becomes  less  significant.  In  fact,  the  limiting  case,  q  =  oo,  leads  to 
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the  deterministic  case  (modulo  sets  of  measure  zero).  Not  surprisingly,  the  results  are 
negative  for  q  =  oo.  In  what  follows,  we  indicate  error  bounds  for  large  M.  Since  we 
always  assume  that  M  <  N,  this  means  that  for  M  tending  to  infinity  we  also  let  N 
tend  to  infinity.  To  make  error  bounds  independent  of  N,  we  take  the  supremum 
over  N  >  M  in  the  corresponding  definitions  of  the  errors.  When  we  speak  about  the 
sharpness  of  error  bounds,  we  usually  take  a  large  M  and  select  a  still  larger  N  and 
a  Boolean  function  for  which  the  presented  error  bound  is  sharp.  The  worst-average 
error  of  the  QS  algorithm  with  M  quantum  queries  satisfies  the  following: 

•  For  q  —  1,  the  worst-average  error  is  ©  (In  M/M).  Furthermore,  the  asymptotic 
constant  is  2/-7T  for  M  —  2  divisible  by  4. 

•  For  q  G  (1,  oo),  the  worst-average  error  is  ©(M~1//<?).  Furthermore,  the  asymp¬ 
totic  constant  is  approximately  ( JJ1  sin9  2(x)clx/ 7r)1/q  for  M  -  2  divisible  by  4 
and  q  close  to  1. 

•  For  q  —  oo,  the  worst-average  error  is  constant  and  equals  1. 

The  error  bounds  of  the  QS  algorithm  are  improved  by  the  use  of  repetitions. 
Namely,  we  repeat  the  QS  algorithm  2n  + 1  times  and  take  the  median  of  the  outputs 
obtained  as  the  final  output.  This  procedure  boosts  the  success  probability  of  the 
approximation  at  the  expense  of  the  number  of  quantum  queries.  We  show  that 
with  n  independent  of  M  and  depending  linearly  on  q,  we  decrease  the  QS  algorithm 
error  to  be  of  order  M-1.  Hence,  the  use  of  repetitions  is  particularly  important  for 
large  q  since  we  change  the  error  bound  0{M~l^q)  without  repetitions  to  the  error 
bound  0{M~l )  with  repetitions.  The  constant  in  the  last  big  O  notation  is  absolute 
and  does  not  depend  on  q. 

The  error  bound  of  order  M_1  is  optimal.  This  follows  from  the  use  of,  for  instance, 
Chebyshev’s  inequality  and  the  fact  that  the  lower  bound  ^(ilW1)  is  sharp  in  the 
worst-probabilistic  setting,  see  also  [Sb].  Hence,  the  QS  algorithm  with  repetitions  is 
optimal  in  the  worst-average  setting. 
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We  now  outline  the  contents  of  this  chapter.  In  Section  12.21  we  define  the  QS 
algorithm.  In  Section  12.31  we  define  precisely  the  error  settings  discussed  above  and 
analyze  the  performance  of  the  QS  algorithm  in  these  error  settings. 

2.2  Quantum  summation  algorithm 

We  consider  the  most  basic  form  of  the  summation  problem,  i.e. ,  the  summation  of 
Boolean  functions.  Let  Bjv  denote  the  set  of  Boolean  functions  /  :  {0, ...  ,N  —  1}  —*■ 
{0, 1}.  Let 


denote  the  arithmetic  mean  of  all  values  of  /.  Clearly,  a/  G  [0,1]. 

Problem.  For  f  G  By,  compute  an  e- approximation  df  of  the  sum  aj  such  that 


\aj  -  “/I  <  E- 


(2.2) 


We  are  interested  in  the  minimal  number  of  evaluations  of  the  function  /  that  is 
needed  to  compute  df  satisfying  (12.211.  It  is  known  that  in  the  worst  case  setting,  we 
need  roughly  IV (1  —  e)  evaluations  of  the  function  /.  In  the  randomized  setting,  we 
assume  that  (if  is  a  random  variable.  We  replace  (j2.2j)  by  the  requirement  that  the 
expected  value  of  \df  —  a/|  is  at  most  e  for  any  function  /.  It  is  known,  see  e.g.,  [33], 
that  in  the  randomized  setting  we  need  roughly  min{IV,  e^1/2}  function  evaluations. 
In  the  quantum  setting,  we  want  to  compute  a  random  variable  df  such  that  (12.211 
holds  with  high  probability  (greater  than  |),  either  for  all  Boolean  functions  or  on 
the  average  with  respect  to  a  probability  measure  defined  on  the  set  By.  These  two 
error  criteria  in  the  quantum  setting  will  be  precisely  defined  in  Section  12.31 


In  this  section  we  describe  the  quantum  summation  algorithm,  which  is  also  called 
the  quantum  amplitude  estimation  algorithm.  This  algorithm  was  discovered  by 
Brassard,  Hpyer,  Mosca  and  Tapp  [4],  and  uses  Grover’s  iterate  operator  as  its  basic 
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component,  see  [13J.  We  use  standard  notation  of  quantum  computation,  see  e.g., 

m- 

For  simplicity  we  assume  that  N  =  2n.  Let  TLn  denote  the  tensor  product  C2  ® 
•  •  •  0  C2  of  n  copies  of  C2,  with  C2  the  2-dimensional  complex  vector  space.  Unit 
vectors  from  C2  are  called  one  qubit  quantum  states  (or  qubits).  Let  |0)  and  |1)  be 
an  orthonormal  basis  of  C2.  Then  any  qubit  \ip)  can  be  represented  as 

\f>)  =  ^o|0)  +  "01 1 1)  with  ifk  G  C  and  |'0o|2  +  l'0i|2  —  l- 

For  j  =  0, 1, . . . ,  TV  -  1,  we  have  j  =  JfVo  2n“1_fcJfc,  with  jk  e  {0, 1}.  Let 

n—1 

I  j)  =  0 \a)- 

k= 0 

The  set  { \j)  :  j  —  0, . . . ,  N  —  1}  forms  an  orthonormal  basis  of  TLn  and  any  unit 
vector  |  -0)  e  Tin  can  be  represented  as 

N- 1  N—\ 

=  with  e  C  and  \ipj\2  =  1. 

3=0  3=0 

Unit  vectors  from  Tin  are  called  n  qubit  quantum  states  (or  quantum  states  or  just 
states,  whenever  n  is  clear  from  the  context). 

The  only  transformations  that  can  be  performed  on  quantum  states  are  defined 
by  certain  unitary  operators  on  TLn.  We  now  define  the  six  unitary  operators  that  are 
basic  components  of  the  summation  algorithm.  Since  unitary  operators  are  linear,  it 
is  enough  to  define  them  on  the  basis  states  | j). 

1.  Let  So  :  Tin  — 1 >  Tin  denote  the  inversion  about  zero  transform 

s0\j)  =  (-i)^Ij), 

where  is  the  Kronecker  delta.  Hence,  S0|0)  =  — 10)  and  S0\j)  =  \ j)  for 
all  j  0.  This  corresponds  to  the  diagonal  matrix  with  one  element  equal 
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to  —1,  and  the  rest  equal  to  1.  We  claim  that  the  operator  So  can  be  also 
written  as  the  Householder  operator 


So  =  I  ~  2 1 0)  (0 1 . 


Here,  for  a  state  |X),  we  let  |X)(X|  denote  the  projection  onto  the  space  span(X)} 
given  by 

(|X)(X|)  X)  =  (XX)  X), 

where  (XX)  is  the  inner  product  in  Hn,  (XX)  =  Yk=o  The  matrix  form 

of  the  projector  |X)(X|  in  the  basis  (|j)}  is  (XfcXXjvXo-  One  can  also  view 
the  matrix  form  of  the  projector  |X)(XI  as  the  matrix  product  of  the  iV  x  1 
column  vector  X)  and  the  iV  x  1  row  vector  (^|  =  |X)X  which  is  the  Hcrmitian 
conjugate  of  X).  To  prove  this  claim,  note  that  for  any  X)  —  Y^o  xj\J)  e  ^ n 
we  have 


X|(/-2|0)(0|)X)  =  {k\x)  -2{0\x){k\0) 


{Xk  —  2xk  =  —  Xk  for  k  =  0, 
Xk  —  0  =  Xk  for  k  X  0. 


Hence,  /  —  2|0) (0|  =  S0,  as  claimed. 


2.  Let  Wn  ■  TLn  — ►  TLn  denote  the  Walsh- Hadamard  transform 

..  n—  1 

WN\j)  =  0  (|0)  +  (-1)“ |1))  . 

That  is,  the  Walsh-Hadamard  transform  corresponds  to  the  matrix  with  entries 


(i\WN\j) 


1 

Xiv 


n— 1 


n 


(i*|(|0)  +  (-l)»|l)) 


1 

Xiv 


n—  1 

fc=0 


1 

_ ( _ 1  ')52fe  =  0  ^kjk 

vXv  J 
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The  matrix  ((i|VFjv|j))ij=o  is  symmetric.  Furthermore, 

n—1 

w%\j)  =  7=W„<g)(|0>  +  (-l)*|l» 

V  ^  k= 0 

1  X  /  1  (-RL  \ 

=  7N®(c2m+m)+-Trm-m)) 

1  n—1 

=  75®^)  =  U>- 

VJV  fc=0 

Thus,  Wjy  =  /  and  IF  ^ 1  =  Wn  is  orthogonal.  This  means  that  the  operator  Wn 
is  symmetric  and  unitary. 

3.  For  K  =  1,2,...,  2n,  let  Fx.n  :  Hn  — >  TLn  denote  the  quantum  Fourier  trans¬ 
form 

(  K-l/2  EK-1  e2irijk/Kjky  for  ./_().! . A'  ~  1 ,  (i  =  ^T) 

Fk,u\3)  =  \ 

{  \j)  for  j  =  K,  ...,2n-l. 

Hence,  Fx.n  corresponds  to  the  unitary  block-diagonal  matrix 

Fk  0 
0  I 

where  Fk  =  {K~1^2  e27ntfc/fc)^  ^  is  the  matrix  of  the  quantum  Fourier  trans¬ 
form.  For  K  =  2 n  =  N  we  have 

N- 1  JV-1  /JV-1  \ 

=  E =  tsE  E i*r>. 

1=0  V  W  k=0  \j= 0  / 

The  coefficients  of  FN,n\if)  in  the  basis  (|j)}  are  the  quantum  Fourier  transforms 
of  the  coefficients  of  the  state  \if).  Note  that  Wn  and  FNn  coincide  for  the 
state  |0),  i.e., 

i  N~ 1 

|o>  =  fn,„  10}  =  -=  E  li>- 

v  v  i=0 
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4.  Let  Sf  :  TLn  —>  TLn  denote  the  quantum  query  operator 

S,\j)  =  (-1)/Wli>, 

This  again  corresponds  to  the  diagonal  matrix  with  elements  ±1  depending 
on  the  values  of  the  Boolean  function  /.  This  operator  is  the  only  one  that 
provides  information  about  the  Boolean  function  /.  This  is  analogous  to  the 
concept  of  an  oracle  or  a  black-box  which  is  used  in  classical  computation  and 
which  supplies  information  about  the  function  /  through  its  values. 

The  standard  definition  of  the  quantum  query  Sf  is 

Sf  :  Hn  ®  C2  -»•  Hn®C2,  Sf\j)\i)  =  \j)\i  © 

where  ©  means  addition  modulo  2.  We  can  simulate  Sf  by  Sf  if  we  use  an 
auxiliary  qubit  (1/a/2)(|1)  —  |0)),  namely, 

=  (-1)«yi>L©o)  =  (5/b))W©o). 

5.  Let  Qf  :  Hn  — >  Hn  denote  the  Grover  operator 

Qf  —  ~WN  5'0  WNl  Sf. 

This  is  the  basic  component  of  Grover’s  search  algorithm,  see  [T3J.  As  we  shall 
see,  Qf  also  plays  a  major  role  for  the  summation  algorithm.  The  eigenvectors 
and  eigenvalues  of  Qf  will  be  useful  in  further  considerations.  Let 

i  N~ 1 

W’>  =  WW|0)  =  ^Ypk) 

and  |^o)>  IV’i)  denote  the  orthogonal  projections  of  j-0)  onto  the  subspaces 
span{|j)  :  / (j )  =  0}  and  span{|j)  :  / (j )  =  1},  respectively.  That  is, 

I i’j)  =  -7=  •/  (i.  i. 

ViV  k:f(k)=j 
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Then  |  -0)  =  \i/;0)  +  \ijji)  and  (V’olV’i)  =  0-  Furthermore,  ('^jl'fpj) 

=  N_1  Yjk-.f{k)=j  !>  for  3  =  0,1,  so  that  (-0i|-0i)  =  a  and  (^ol^o)  =  1  -  a, 
where  a  =  aj  is  the  sum  we  want  to  approximate. 

From  [3],  we  know  that 

Q/IVh)  =  (1  -  2a)|Vh)  +  2(1  -  a)|V>i), 

Q/|^i)  =  -2a|^o)  +  (1  -  2a)|^i). 

For  the  sake  of  completeness,  we  provide  a  short  proof  of  (I2.3[).  By  the  definition 
of  the  operator  Sf  we  have 

=  (-1)JlV’j),  3  =  0,1, 

and 

IF/v^PF^1  =  kF/v(/  —  2|0)(0|)  PF^1  =  /  -  2(Wjv|0)(0|Wjv). 

Since  (0|VFat  =  (VFjv|0))^  =  (|-0))^  =  {tp\,  we  obtain 

WNS0W^j)  =  |^>  -  2(|^)(^|)  |^-) 

=  iVh)  -  2(V#j)|V>)  =  IV’j)  -  2(Vh#i)l# 


for  j  =  0, 1.  From  this  we  calculate 

Qt\*i)  =  (-1  Y+iWuScW^Wi) 

=  (IV’j)  -  2(tij,ia  +  <5,.0(1  -  «))(|Vo)  +  IVO)) , 

for  j  =  0, 1,  which  is  equivalent  to  (12. 3 j). 

Thus,  the  space  span{|^0),  |t/h)}  is  an  invariant  space  of  Qf  and  its  eigenvectors 
and  corresponding  eigenvalues  can  be  computed  by  solving  the  eigenproblem  for 
the  2x2  matrix 

1  —  2a  —2a 
2(1 -a)  1  —  2a 
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If  a  G  (0, 1),  then  the  eigenvalues  of  Qf  are 

A±  =  1  —  2a  ±  2%\J a{\  —  a)  =  e±2lda,  9a  =  arcsin  >fa, 
and  the  corresponding  orthonormalized  eigenvectors  are 

lv,±)  =  ^  (Wi =«  +  ^w)' 

Moreover,  it  is  easy  to  check  that 

W  =  =^(e«^+)-e-^.)).  (2.4) 

If  a  G  {0, 1},  then  we  have  span{|'0o))  |Vi)}  =  span{|V)}  and  IV7)  is  the  eigen¬ 
vector  of  Qf  with  eigenvalues  ±1,  respectively.  For  a  G  {0, 1},  we  define 

|V+)  =  V"“\/2|V)  and  |V-)  =  0. 

Then  it  is  easy  to  check  that  (12.4jl  is  valid,  and  A-t  =  e±2*0a  =  (—1)“  is  an 
eigenvalue  of  Qf  for  all  a  G  [0, 1]. 

6.  The  next  unitary  transform,  called  the  Grover  iterate  operator ,  is  defined  on  the 
tensor  product  of  TCm  ®  Ttn  and  uses  m  +  n  qubits.  The  number  m  is  related  to 
the  accuracy  of  the  quantum  summation  algorithm,  whereas  n  is  related  to  the 
size  of  the  problem.  The  Grover  iterate  operator  A m{Qf)  '■  'Rm®‘Hn  'Hra®'Hn 
is  defined  by 

Kn{Qf)  I j) \y)  =  I j)  Q}  I y)  for  | j)  \y)  Un. 

Hence,  the  power  of  Qf  applied  to  the  second  component  depends  on  the  first 
one.  Note  that  j  may  vary  from  0  to  2m  —  1.  Therefore  A m(Qf)  rnay  use  the 
powers  of  Qf  up  to  the  (2m  —  l)st. 

We  need  one  more  concept  of  quantum  computation,  that  of  measurement.  Sup¬ 
pose  s  is  a  positive  integer  and  consider  the  space  Tis.  Given  the  state 

2s  — 1 

IV)  =  Vfclfc)  e  Hs, 

k= 0 
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we  cannot,  in  general,  recover  all  the  coefficients  We  can  only  measure  the 
state  \^j)  with  respect  to  a  finite  collection  of  linear  operators  {Mj}j=0,  where  the  Mj  : 
Hs  — >  7is  satisfy  the  completeness  relation 

3=0 

After  performing  the  measurement,  we  obtain  the  outcome  j,  and  the  state  \ij))  col¬ 
lapses  into  the  state 

1  =Mj\t/i); 

this  occur  with  probability  Note  that  for  =  0  the  outcome  j 

cannot  occur  with  positive  probability.  Hence,  with  probability  1  the  outcome  j 
corresponds  to  ^  0. 

The  most  important  example  of  such  a  collection  of  operators  is  {|j)(j  Ijylo1- 
Then,  the  measurement  of  the  state  \ij))  with  respect  to  this  collection  of  operators 
gives  us  the  outcome  j  and  the  state  |^)  collapses  into  the  state 

ML\j) 

i  mr 

with  probability  |^j|2,  j  =  0, 1, . . . ,  2s  —  1. 

Another  example  is  a  variation  of  the  previous  example  and  will  be  used  in  the 
quantum  summation  algorithm.  We  now  let  s  —  rn  +  n,  as  for  the  Grover  iterate 
operator,  and  define  Mj  :  Ttm  ®  TLn  — >  Hm  (8)  Tin  by 

=  |j)0'|®/ 

for  j  —  0,1, ,  2m  —  1,  with  /  denoting  the  identity  operator  on  Hn.  That  is, 

{\j)(j\  ®J)  \x)\y)  =  (j\x)  \j)\y) 
for  \x)  G  7im  and  | y)  G  Hn. 

Since  ®  -0  \x)\y)  =  \x)\y)  f°r  all  basis  states  |a;)  of  7 -Lm  and  | y)  o!TLn, 

the  completeness  relation  is  satisfied.  Consider  now  the  probability  of  the  outcome  j 
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for  a  special  state  | i/j)  of  the  form  | -0)  =  |'0i)|t/;2)  with  |-0i)  G  7Tm,  \^2)  £  7Yn,  where 
(V;fe|V;fc)  =  1  for  k  =  1,2.  Since  \j)(j\®I  is  self-adjoint,  the  outcome  j  and  the  collapse 
of  the  state  | if})  to  the  state 


OlV’i) 


b')  l^2> 


I  Wi>  I 

occur  with  probability  |(i|V,i)|2-  Hence,  this  collection  of  operators  measures  the 
components  of  the  so-called  first  register  1^)  of  the  quantum  state  \^). 


Following  [3],  we  are  ready  to  describe  the  quantum  summation  (QS)  algorithm 
for  solving  our  problem.  The  QS  algorithm  depends  on  a  Boolean  function  /  and 
on  an  integer  parameter  M  that  controls  the  number  of  quantum  queries  and  the 
accuracy  of  the  algorithm.  We  perform  computations  in  the  space  Tim  <S>  Hn,  with 
m  =  [log 2M],  so  we  use  m  +  n  qubits.  As  we  will  see  later,  the  accuracy  of  the 
algorithm  is  related  to  the  dimension  of  the  space  Tim. 


Algorithm  QS (/,  M) 

Input  state:  |0)|0)  G  Ttm  ®  Tin  with  m  =  [log2  M]  and  n  =  log2  N. 
Computation: 

1-  \vi)  =  FM,m®WN\0)\0), 

2.  \rj2)  =  A m(Qf)  \vi), 

3-  \V3)  =  (F^m^I)\r]2). 

Measurement: 

Perform  the  measurement  of  the  state  I773)  with  respect  to  the  collection 
{  ( |i )  (,'j | )  ®  I  }2j=o  *  •  Denote  the  outcome  by  j. 

Output:  af(j)  =  sin2  (7 rj/M). 


We  briefly  comment  on  the  QS  algorithm.  The  input  state  is  always  the  same  and 
does  not  depend  on  /.  Step  1  computes  (771)  =  (N ,  which 
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is  the  equally  weighted  superposition  of  the  basis  states.  Step  2  computes  (772)  by 
using  the  Grover  iterate  operator.  During  this  step,  we  use  the  successive  powers  of 
the  Grover  operator  Qf ;  this  is  the  only  step  where  information  about  the  Boolean 
function  /  is  used.  We  shall  see  that  the  QS  algorithm  uses  M  —  1  quantum  queries. 
Step  3  computes  (773)  by  performing  the  inverse  quantum  Fourier  transform  on  the 
first  m  qubits,  and  prepares  the  system  for  measurement.  After  Step  3,  we  perform 
the  measurement,  obtain  the  outcome  j  and  compute  the  output  a/(j)  on  a  classical 
computer.  We  stress  that  the  distribution  of  the  outcomes  j  depends  on  the  Boolean 
function  /,  and  this  is  the  only  dependence  of  the  output  a/(j)  on  /. 


2.3  Performance  analysis 

In  this  section  we  analyze  the  error  of  the  QS  algorithm.  As  we  have  seen  in  Sec¬ 
tion  |2]2l  the  output  a/(j)  of  the  QS  algorithm  is  a  random  value  chosen  according 
to  a  certain  distribution  that  depends  on  the  input  function  /.  I11  this  sense,  the  QS 
algorithm  is  a  randomized  algorithm.  Various  ways  of  measuring  the  performance 
of  randomized  algorithms  are  commonly  used  in  the  analysis  of  algorithms  and  com¬ 
putational  complexity;  these  all  correspond  to  various  error  criteria.  In  particular, 
we  consider  three  error  criteria:  worst-probabilistic,  average-probabilistic  and  worst- 
average. 

Worst-probabilistic  error 

We  start  with  the  error  criterion  that  is  used  in  most  papers  dealing  with  quantum 
computation.  We  are  interested  in  the  worst  case  error  of  the  QS  algorithm  that 
holds  with  a  given  probability  p.  Here  pg  [0,1]  and  1  —  p  measures  the  probability 
of  QS  algorithm’s  failure  and  usually  p  is  set  to  be  I11  our  analysis,  however,  we 
will  allow  an  arbitrary  p  G  (1/2,  8 / 7r2] .  The  choice  of  the  upper  bound  8 / 7r2  =  0.81 . . . 
will  be  clear  from  the  analysis  of  the  QS  algorithm.  The  QS  algorithm  outputs  a/(j), 
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j  —  0, 1, . . . ,  M  —  1,  with  probability  Pf(j),  which  is  given  in  (12.61)  of  Theorem  12.3.21 
The  worst-probabilistic  error  of  QS  is  formally  defined  as  the  smallest  error  bound 
that  holds  for  all  Boolean  functions  with  probability  at  least  p,  i.e. , 

ewor-pro  (MjP)  =  inf  ja  :  Y  PfU)>P  V/  G  Bjvj. 

It  is  easy  to  see  that  ewor_pro(M, p)  can  be  rewritten  as  follows.  Let  J  C  {0,1 ,M  — 
1}.  For  /  e  Mn  define  the  measure  of  J  as 


p(jj)  = 

16J 

Then 

ewor-pro(M,p)  =  max  min  max \at  —  af(j)\ 
f&N  J-n(J,f)>p  j£J 


(2.5) 


Average-probabilistic  error 

The  worst-probabilistic  error  ewor"pro(M,  p)  of  the  QS  algorithm  is  defined  by  the 
worst  performance  with  respect  to  Boolean  functions.  It  is  also  natural  to  consider 
the  average  performance  of  the  QS  algorithm  with  respect  to  Boolean  functions.  Let  p 
be  a  probability  measure  on  the  set  MN,  so  that  each  Boolean  function  /eBjv  occurs 
with  probability  p(/).  Obviously,  p(/)  >  0  and  Xl/eiw  P(/)  =  1-  The  average- 
probabilistic  error  is  defined  by  replacing  the  first  max  in  (12.51)  by  the  expectation, 
i.e., 

eavg-pro (M,p)  =  ^  P (/)  min  max  \af  -  af(j)\, 
f&N 

Hence,  we  are  interested  in  the  average  error  that  holds  with  a  certain  fixed  proba¬ 
bility. 


Worst-average  error 

The  worst-average  error  corresponds  to  the  worst  case  performance  with  respect  to  all 
Boolean  functions  from  IB tv  and  the  average  performance  with  respect  to  all  outcomes. 
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This  average  performance  is  measured  by  the  expectation  in  the  Lq  norm,  q  £  [1,  oo], 
with  respect  to  the  probability  measure  of  all  outcomes  provided  by  the  QS  algorithm. 
As  mentioned  before,  we  make  the  worst-average  error  independent  of  N  by  taking 
the  supremum  over  N  >  M.  That  is,  the  worst-average  error  is  defined  as: 

•  for  q  £  [1,  oo), 

(M—  1 

£  (i)r 

r, 

J=0 

•  for  q  —  oo, 

ewor-avg(M)  =  sup  |0/  _ 

N>M  /£®iv  J'-Pf  (j)>0 

It  is  easy  to  check  that  for  q  =  oo,  the  QS  algorithm  behaves  badly.  Indeed,  if  M  is 
odd,  we  can  take  /  with  all  values  one,  and  then  a/  =  1,  pf( 0)  =  1/M 2  and  df( 0)  =  0. 
Hence  e™or~avg(M)  =  1.  If  M  is  even,  we  take  /  with  only  one  value  equal  to  1,  and 
then  cif  =  1/N,  Pf(M/ 2)  >  0  and  df(M/ 2)  =  1.  Hence,  | a /  —  df(M/ 2)|  =  1  —  1/Af 
and  e^or-avg(M)  =  1. 

Therefore,  in  the  rest  of  the  chapter,  we  consider  only  q  £  [1,  oo).  As  we  shall  see, 
the  cases  q  >  1  and  q  —  1  will  require  different  analyses  and  lead  to  quite  different 
results. 


2.3.1  Worst-probabilistic  error 

We  begin  by  citing  a  theorem  from  [4]  for  which  we  shall  propose  a  number  of  im¬ 
provements. 


Theorem  2.3.1.  For  any  Boolean  function  f  £  Btv7  the  QS  algorithm  uses  ex¬ 
actly  M  —  1  quantum  queries  and  outputs  a  that  approximates  a  —  a/  such  that 


2vr 


a  —  a 


AJi  / — - T  71  71  71 
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with  probability  at  least  S/n2  =  0.81 ....  Hence, 


E 

3-  \af(j)-af\<(2n/M)ffaf{l-af)+iv2/M2 


Pfti)  >~o  V/e 


7 r 


»iV, 


and,  therefore, 


3 w°r-pro  I  M  <  iL  I 

’W  -  M  +  M2' 


Using  ideas  from  the  proof  of  Theorem  12.3.11  from  |1J  we  present  the  following 
theorem  and  the  subsequent  corollaries. 

Theorem  2.3.2.  For  any  Boolean  function  f  E  Bat,  denote 

M 

u a  =  (7af  =  — arcsin  V a  E  [0,  \  . 

1.  The  QS  algorithm  uses  exactly  M  —  1  quantum  queries,  and  log2  N  +  |~log2  M] 
qubits. 

2.  For  j  =  0, 1, . . . ,  M  —  1,  the  outcome  j  of  the  QS  algorithm  occurs  with  proba¬ 
bility 


Pfti)  = 


sin2(7r(j  -  a a  )) 


2  M2  sin2  (irti  -  aaf)/M) 


1  + 


sm" 


sin' 


NixfiiLFl)  (2  6) 
Wi  +  f  ['b> 


(If  sin  (t t(J  ±  °af)/M)  =  0  we  need  to  apply  the  limiting  value  of  the  formula 
above.)  For  j  =  M ,  M  +  1, . . . ,  2l"log2  —  1,  the  outcome  j  occurs  with  proba¬ 
bility  0. 


3.  If  a a f  is  an  integer,  then  the  QS  algorithm  outputs  the  exact  value  of  a /  with 
probability  1.  This  holds  iff  aj  =  sin2(/c7r/M)  for  some  integer  k  E  [0,  \M\.  In 
particular,  this  holds  for  aj  =  0,  for  a/  =  1  and  even  M ,  and  for  aj  —  )  and 
M  divisible  by  4. 


4-  Letx  =  7r(|"(ja]  —aa)/M  andx  =  tt (<ra  —  [aa\)/M.  If  <jaf  is  not  an  integer,  then 
the  QS  algorithm  outputs  the  same  value  a  =  a/([cra])  =  df(M  —  [cra])  for  the 
outcomes  [cra]  and  M  —  \aa~\  such  that 


a— a  = 


sin(x)  ^2 a/ a(l  -  a)  cos(x)  +  (1  -  2a)  sin(x)^  <  -^(  [cra]-cra)  (2.7) 
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with  probability 


sin2(7r(  |~cra]  -  cra)) 

M2  sin2  (vr(  [cxa]  -  <Ja)/M ) 

.  o 
sir 

> 


,  /i  jt  a  sin2  (7r(  [cra] 

1  +  (I  -  0\aa\,M/V^^2 - 


0/M)\ 


sin2  (7r(  [cra]  +  aa)/M)J 


Sm"  |Vr(r<Ta]  —  <7aJJ  ,  7T 


=  1  -  T(M  - + (2-8) 


and  outputs  the  same  value  a  =  a/(|_cra J)  =  a/((  1  —  <5|_t7oj)o)M  —  L^aJ)  for  th e 
outcomes  [aa\  and  (1  —  —  [aa\such  that 


■a  = 


sin(x)  ^2 a/ a(l  -  a)  cos(x)  +  (1  -  2a)sin(x)j  <  -^(cra-  |yv|)  (2.9) 


with  probability 


sin2  (7r(era  -  |aaj)) 

M2  sin2  (7r(cra  -  [<Ja\)/M) 
>  sin2  (7r(cra  -  |aaj)) 

_  7T2(0-a-LCTa  J)2 


n_c  \  sin2  (7r(<ra-  [craJ)/M)\ 

kj,°  sjn2  +  L<jaJ)/M)y 


LcxaJ)2  +  0((aa  -  L^aJ)4). 


(2.10) 


Proof.  As  before,  let  6a  =  arcsin  y/a  and 

1  M— 1 

|S„M)  =  -j=  J2  <?~k\k),  i  =  V=i, 

for  arbitrary  u6l.  Note  that 


Fftf,m\j) 


1 1 SmU/M))  for  j  =  0, 1, . . . ,  M  —  1, 

I  |  j)  for  j  =  M,  M  +  1, . . . ,  2m  —  1. 


The  steps  1-3  of  the  QS  algorithm  are  equivalent  to  the  application  of  the  opera¬ 
tor  {F^m  <g)  I)  A m{Qf)  FM,m  ®  WN  to  the  state  |0)|0)  G  Hm  <g>  Hn.  Now  \rji)  can  be 
written  as  M~1^2^2jI=01  \ j)\ij)),  where  | ijj)  =  Wjv|0)  is  given  by  (12.41).  Hence 

.  M—l 

\vi)  =  -J=  {ei9a\^+)  -  e~i9a\^-))  ■ 

V  3=0 
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Applying  A m(Qf)  in  Step  2  and  remembering  that  Q}|0±)  =  A±|0±),  we  obtain 

.  M—l 

M  =  A „(<?,)  M  =  \i)  -  e-2«8*e-"-|V--» 

V  3=  0 

=  A  (eie‘\SM(<rJM))\-i,+)  -e-‘e‘\SM(-^/M)W-)) . 

Since  j  =  0, 1, . . . ,  M  —  1,  the  largest  power  of  Qf  is  M  —  1.  Hence,  we  nse 
exactly  M  —  l  quantum  queries  to  compute  \rj 2).  The  remaining  steps  of  the  QS 
algorithm  do  not  use  quantum  queries.  This  means  that  the  total  number  of  quantum 
queries  used  by  the  QS  algorithm  is  M  —  1,  and  obviously  we  are  using  n  +  m  qubits. 
This  proves  the  first  part  of  Theorem  12.3.21 
Step  3  yields  the  state 

M  = 

We  are  ready  to  analyze  the  probability  of  the  outcome  j  of  the  QS  algorithm. 
Observe  that 

l«±)  :=  (Ij)O'I  0-0  F^m\SM(±^a/M))\ip±)  =  O'  I  \  SM  (±cra/M) )  |j)|0±) 

[  (SM(j  /M)\SM(±CTa/M))  |j)|0±)  for  j  =  0, 1,...,M-  1, 

I  0  for  j  =  M,  M  +  1, . . . ,  2m  —  1, 

and  therefore 

(\(SM(j/M)\SM(±<Ta/M))\2  (0±|0±)  for  j  =  0, 1, . . .  ,M  -  1. 

<a±|a±)  =  < 

lo  for  j  =  M,  ...,2m  — 1. 

Observe  that  for  a  G  (0,1),  we  have  (0±|0±)  =  1)  whereas  for  a  G  {0,1},  we 
have  (0+I0+)  —  2  and  {0— 10—)  =  0- 
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For  Wi,w26Rwe  have 


|(5'm(^i)|5'm(co’2))|2  — 


1  M—l 

—  V' 

VmR,' 


M—l 


3— 2  niujij 


01 


E  e2niU2j\  j) 


M2 


M- 1 

E 

3=0 


2tti(uji—  uj2)j 


If  uJi  —  lv2  is  an  integer,  then  the  last  sum  is  clearly  M,  and  the  whole  expression  is  1. 
If  U\  —  uj 2  is  not  an  integer,  then 


1 

M 


M- 1 

E 

3=0 


e—2niM(u>i—LU2)  _ 

M(e~2niO i_a,2)  _  iy 


which  holds  for  all  wi,W2  GKif  we  take  0/0  as  1.  Therefore 


M 


Hence 


M—l 

E' 

3=0 


,-2iri(u>i-W2)j 


1  —  cos(27 tM  (loi  —  0J2))  sin2  (77  M  {uj\  —  0*2)) 


M2(  1  —  cos(27r(cni  —  cn2)))  M2(sm2(n(uJi  —  uj2))) 


2  sin2(M7r(o;i  —  u2)) 


M 2  sin2(7r(o;i  —  cn2))  ’ 
which  holds  for  all  Wi,W2  €  1  if  we  take  0/0  as  1.  Applying  this  we  conclude  that 

sin2(7r  (jt^o)) 


(2.11) 


(a±|a±)  =  {  M2sm2(n(jT<Ja)/M) 
0 


(0±I0±)  for  j  =  0, 1,...,M-  1, 

for  j  —  M,M  +  1, . . . ,  2m  -  1. 


The  outcome  j  occurs  after  the  measurement  and  the  state  (773)  collapses  to  the 
state 

( (h3 1 Mj  Mj 1 773 ) ) “ 1  Mj 1 773 ) ,  where  Mj  =  \j){j\  <8>  /.  For  j  =  0, 1, . . . ,  M  —  1,  we  have 


Mite)  =  \])(-y(eie‘(SM(,IM)\SM(aJM))\^) 

-  e-"“{SM(j/M)|SM(-^/M)>  |V-»  V 

whereas  Mj  |r/3)  =  0  for  j  =  M ,  M  + 1, . . . ,  2m  —  1.  Since  |^+)  and  |^_)  are  orthogonal 
we  have 


(^3 1  AfjMj- 1 773)  =  |((q!+|q!+)  +  («_!«-))• 
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We  now  claim  that  for  any  j  G  {0, 1, . . . ,  M— 1},  the  outcome  j  occurs  with  probability 

n  _  1  f  sin2(7r(j  -  ora))  sin2(7r(j  +  aa)) 

^  2  \M2sin2(7r(j  ■—  aa)/M)  M2  sin2  (n(j  +  aa)/M) 

Indeed,  for  a  G  (0,1),  we  have  =  1  and  (I2.12p  follows  from  the  form  of 

(a±|a±).  For  a  G  {0, 1},  we  have  (^+|'0+)  =  2  and  =  0.  Since  the  two  terms 

in  (12. 12(1  are  now  the  same,  the  formula  for  (a+|a+)  again  yields  (12. 12(1 . 

Since  sin2(7r(j  —  aa ))  =  sin2(7r(j  +  cra)),  the  last  formula  is  equivalent  to  (12.611. 
Obviously  for  j  =  M,  M  +  1, . . . ,  2m  —  1,  the  probability  of  the  outcome  j  is  zero. 
This  proves  the  second  part  of  Theorem  12.3.21 

Assume  now  that  aa  G  Z.  If  aa  =  0  or  <ja  =  (if  M  is  even),  then  the 
probability  Pf(cra)  of  the  outcome  aa  is  1.  For  <ja  =  0  we  have  a  =  0  and  the  output 
is  a/(0)  =  0.  For  <ra  =  \M  we  have  a  =  1  and  the  output  is  a/(|M)  =  1.  Hence,  in 
both  cases  the  QS  algorithm  outputs  the  exact  value  with  probability  1. 

If  a a  G  Z  and  aa  ^  (0,  t,M},  then  the  probability  of  the  distinct  outcomes  aa 
and  M  —  aa  is  These  two  values  of  the  outcomes  yield  the  same  output 

sin2  {nOa/M)  =  sin2  (n(M  —  aa)/M)  =  a. 

Hence,  the  QS  algorithm  outputs  the  exact  value  with  probability  1.  This  proves  the 
third  part  of  Theorem  12.3.21 

We  now  turn  to  the  case  when  aa  ^  Z.  It  is  easy  to  check  that  the  third  part 
of  Theorem  12.3.21  holds  for  M  —  1.  Assume  then  that  M  >  2,  which  implies  that 
\\M\  <  M  —  1.  Since  aa  is  not  an  integer,  we  have  [cra]  >  1,  [cra]  <  | \\M\  <  M  —  1 
and  M—  \aa~\  <  M—  1.  This  means  that  both  \aa~\  and  M—  \aa~\  may  be  the  outcomes 
of  the  QS  algorithm.  Obviously,  these  two  outcomes  are  different  iff  [cra]  \M. 
Similarly,  both  [aa\  and  (1  —  8\pa\ ,o)M  —  [aa J  may  be  also  the  outcomes.  They  are 
different  iff  [<ja J  ^  0. 


We  show  that  the  outputs  for  the  outcomes  |~cra~|  and  [aa\  satisfy  (12.71)  and  (12.91) 
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with  probabilities  (12 ,8ft  and  (12.10 jl.  respectively.  We  focus  on  the  output  for  the 
outcome  |~<ra]  and  its  probability.  The  proof  for  the  outcome  \_aa J  is  similar. 

We  estimate  the  error  of  the  QS  algorithm  for  the  output  a  =  sin2(7r  \aa~]  /M). 
Recall  that  x  =  7r(|~cra]  —  oa)/M.  We  have 


a  —  a |  =  |  sin2(7r  [<ra]  /M)  —  sin2(7rcxa/M)|  =  |  sin(x)  sin(x  +  27rcra/M)| 

=  |  sin  (a;)  (  sin(27rcra/M)  cos(a;)  +  cos(27rcra/M)  sin(x))  | 

<  ?r(  \&a  \  -  cra)/M. 


Since  sin(27rcra/M)  =  2yja(l  —  a)  and  cos(27t aa/M)  =  1  — 2a,  this  proves  the  estimate 
of  the  error  of  the  QS  algorithm  in  the  fourth  part  of  Theorem  12.3.21 

We  find  the  probability  of  the  output  a.  Since  t  h- >•  sin2  (nt  /  M)  is  injective  for  t  6 
[0,  |M],  the  output  a  occurs  only  for  the  outcomes  |"cra]  and  M  —  \ ~aa].  If  |"cra]  =  |M, 
then  these  two  outcomes  are  the  same  and  a  occurs  with  probability  Due  to 

(I2.12h  we  have 


(1  ,  _  sin2  (tt(| M-ga/)) 

Pf[2  ’  M2sin2  (7t(|M  —  cra/)/M)  ’ 
which  agrees  with  the  claim  in  Theorem  12.3.21 

If  [<7aj  7^  |M,  then  |"cra]  j -  M  —  \aa~\  and  a  occurs  for  exactly  two  distinct  out¬ 
comes.  The  probability  of  a  is  now  equal  to  the  sum  of  the  probabilities  £>/([cra])  + 
Pf(M  —  |~cra]),  where  pj  is  given  by  (12.121).  Since  both  terms  are  equal,  the  prob¬ 
ability  of  a  is  2pf(\cra})  which  also  agrees  with  the  claim  in  Theorem  12.3.21  Since 
sin  (7r([<7a~|  -  cra)/M)  <  7r([cra]  -  cra)/M  we  have 

sin2(7r(  [cra]  -  a a))  sin2(7r  [cra]  -  ga) 

M2  sin2  (7r([cra]  ^  aa)/M)  ~  n 2(\aa]  -  aa)2 

We  hnish  proving  (12.81)  by  using  the  standard  expansion  of  the  sine.  This  completes 

the  proof.  □ 


Based  on  Theorem  12.3.21  we  present  simplified  estimates  of  the  error  of  the  QS 
algorithm  and  of  the  corresponding  probability. 
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Corollary  2.3.3.  The  QS  algorithm  outputs  a  such  that 


7 r 


\a  ~  <*l  <  T7  max{  \aa]  -  aI/:  aa  -  [<?„ J } 


(2.13) 


with  probability  at  least 

sin2(7r(  |"<ra]  -  cra)) 


+ 


sin2(7r(aa  -  |aaJ)) 


>  — . 


(2.14) 


M 2  sin2  (7r(  \aa]  -  cra)/M )  M2  sin2  (7r(cra  -  [oa\  )/M)  7r2  ' 

Proof.  It  is  enough  to  prove  Corollarvl2.3.3lif  ua  is  not  an  integer.  The  estimate  of  the 
error  of  the  QS  algorithm  by  the  maximum  of  the  estimates  (12.71)  and  (I2.9j)  holds  with 
probability  that  is  the  sum  of  the  probabilities  (12. 8j)  and  (12.101).  Moreover,  [cra~|  —  aa  = 
1  —  (cra  —  |_craJ).  Observe  that 

=  ,  sin2(7r(l  -  A)) 

9[  ’  vr2A2  tt2(1-A)2 

is  a  lower  bound  of  the  left  hand  side  of  (I2.14p  with  A  =  | Va]  —  cra.  Since  the 
function  g  attains  the  minimum  8/7 r2  on  the  interval  [0, 1]  for  A  =  see  also  j0],  this 
completes  the  proof.  □ 


Corollary  12.3.31  guarantees  that  the  estimate  ()2. 13jl  holds  with  high  probability. 
Unfortunately  this  estimate  does  not  preserve  the  continuity  of  the  estimates  (12.71) 
and  (12. 9p  with  respect  to  [cra~|  —  aa  and  oa  —  [_craJ .  The  continuity  of  the  estimates  is 
present  in  the  next  corollary,  at  the  expense  of  the  probability  of  the  outcome.  This 
corollary  will  also  play  an  essential  role  in  the  study  of  the  average-probabilistic  error 
of  the  QS  algorithm. 


Corollary  2.3.4.  The  QS  algorithm  outputs  a  such  that 


|a  —  a|  <  min{  (" cra] 


C  a.  1  &  n 


M} 


with  probability  at  least 

sin2(7r([aa]  -  aa))  sin2(7r(cra  -  |oa_|))  )  4 

M2  sin2  ( ( \(Ja  \  -  *a)YM2  sin2(§K  -  KJ))  I  “  ^2' 


(2.15) 


max 


(2.16) 
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Proof.  We  may  again  assume  that  cra  is  not  an  integer.  Let  us  define 


w(A) 


sin2(7rA) 

M 2  sin2  (ttA/M) 


for  A  G  [0, 1]. 


Then  w{\c ra]  —  aa)  is  the  probability  that  (I2.7|)  holds  and  w(l  —  (faa~|  —  cra))  is  the 
probability  that  (12. 9  p  holds.  For  A  G  [0,  |],  note  that  w(-)  is  decreasing,  and  w(l  —  •) 
is  increasing.  Therefore 


w( A)  >  w{\)  >  tc(l  —  A)  for  A  G  [0, 

Suppose  that  \aa]  —  cra  <  cra  —  [aa J.  Then  [aa]  —  cra  <  \.  In  this  case  (12. 15ft  is 
equivalent  to  m,  and  holds  with  probability  at  least  w(  \aa~]  —cra),  which  corresponds 
to  (12 . 1 6[) .  Analogously,  if  |~cra]  —  cra  >  aa  —  L<raJ  then  \aa]  —  oa  >  In  this  case 
(I2.15P  is  equivalent  to  (I2.9j) .  and  holds  with  probability  at  least  w(aa  —  L°aJ),  which 
also  corresponds  to  (I2.16p.  Finally,  note  that 


max {w( -  <ja),w(cra  -  [aaJ)} 


is  minimal  for  [<ra]  —  cra  =  |  and  is  equal  to 


sm 


-2 


7T 


>  — . 


M 2  2 M  ~  it2' 


□ 


Unfortunately,  for  [<ja 1  —  cr a  close  to  \  the  probability  of  the  estimate  (12 . 1 5ji  is 
too  small.  However  in  this  case  we  may  use  Corollary  12.3.31  which  yields  an  estimate 
with  high  probability. 

We  now  turn  to  global  error  estimates,  that  is,  estimates  independent  of  a.  The¬ 
orem  izniof  g  states,  in  particular,  that  \a  —  a\  <  n/M  +  tt2/M2  with  probability 
at  least  8/7 r2.  We  now  improve  this  estimate  by  combining  the  estimates  (12.131) 
and  (12451). 
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Corollary  2.3.5.  The  QS  algorithm  outputs  a  such  that 


3  i r 

\a  —  a  \  <  -  — 
1  “41 


(2.17) 


with  probability  at  least  S/n2.  That  is, 


0wor-pro  I  M  ^ 

’  7T2  /  -  4  1' 


Proof.  Let  us  define 


Sill' 


!(7tA)  sin2(7r(l  —  A))  1 


(2.18) 


ft(A)=maxl__,  ^(1_A)2  f 

Clearly,  h(|~cra]  —  aa)  is  a  lower  bound  on  max{ta([<ja]  —  aa),w(  1  —  ([cxa]  —  cra) ) }  and 
therefore  h(\aa~\  —  aa)  is  a  lower  bound  of  the  probability  of  the  output  satisfying 
(I2.15I).  We  consider  two  cases. 

Assume  first  that  A  =  |"cra]  —  cra  G  [0,  |]  U  [|,1].  It  is  easy  to  see  that  then 
h{ A)  >  8 / 7T2  and  the  estimate  (12.151)  yields 

\a  -  a\  <  ^  min{  |"cra]  -  aa,  aa  -  [craJ }  <  ^  ^ 

with  probability  at  least  8/7 r2. 

Assume  now  that  |"(Ta]  —  cra  G  (|,  |).  Then  we  can  use  the  estimate  (12.131).  which 
holds  unconditionally  with  probability  at  least  8 / 7T2 .  In  this  case,  we  have 

\a-a\<^  max{  \aa~\  -  aa ,  aa  -  L^J }  <  ^  jj- 

These  two  estimates  combined  together  yield  (12.171).  □ 


An  obvious  consequence  of  Corollary  12.3.51  is  that  for  M  large  enough  we  can 
compute  the  value  of  a  exactly  by  rounding  the  output. 


Corollary  2.3.6.  Assume  that 


M  >  —N. 
2 


Then  the  rounding  of  the  QS  algorithm  output  to  the  nearest  number  of  the  form  k/N 
yields  the  exact  value  of  the  sum  a  with  probability  at  least  8 / 7r2 . 
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The  proof  of  Corollary  12.3.51  may  suggest  that  the  constant  /  in  (12 . 1 7jl  can  be 
decreased.  Furthermore  one  may  want  to  decrease  the  constant  |  at  the  expense  of 
decreasing  the  probability  8/7 r2.  These  points  are  addressed  in  the  next  corollary.  We 
shall  see  that  the  constant  |  may  be  lowered  only  by  decreasing  the  probability. 


Corollary  2.3.7.  Define 


7 r 


C{p)  =  inf  \  C  :  | cif  —  a/|  <  C—  V/  G  B^v  with  probability  at  least p 


and 


v(A)  = 


sin2(7rA) 

7T2A2 


for  A  G  [|,  3]. 


Then 


1 

2 

for 

P 

G  [0,4/tt2), 

C(p) 

<  < 

1  —  ir1 

ip) 

for 

V 

G  [4/tt2,  8/tt2], 

AI/n 

< 

for 

p 

e  (8/ 7T2,  1], 

Moreover,  1  - 

-v  \p)  G  [ 

14] 

and 

7T2  1 

16P+4' 

(!- 

-v~l{p)) 

< 

0.0085 

for  p  G  [4/ 7f2 

(2.19) 


7T2].  (2.20) 

Proof.  For  p  G  [0,4/7 r2),  Corollary  12.3.71  is  a  consequence  of  Corollary  12.3.41  For  p  G 
(8 / 7r2, 1] ,  Corollary  12.3.71  trivially  holds  since  |a  —  a|  <  1  =  (M/tt)tt/M.  For  the 
remaining  p’s  we  use  a  proof  technique  similar  to  that  of  Corollary  12.3.51 

Let  p  G  [4/ 7T2,  8 / 7T2] .  It  is  easy  to  check  that  v  is  decreasing  and,  therefore,  v~1(p) 
is  well  defined  and  v~1(p)  G  [|,  |].  We  have  to  show  that  the  estimate 


-a|  <  (!-t.  \p))fj 


(2.21) 


holds  with  probability  at  least  p.  We  consider  two  cases. 

Assume  first  that  A  =  \aa~]  —  aa  G  [0,  v_1(p)]  U  [1  —  v~1(p),  1].  Observe  that  the 
function  h  defined  in  (12. 18f)  can  be  rewritten  as 


h(  A)  =  max{v(A),u(l  —  A)}. 
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It  is  easy  to  see  that  in  this  case,  h( A)  >  p,  and  the  estimate  (I2.15j)  yields 

\a  —  a\  <  min{A,  1  —  A}  <  v~1(p)  <  (l  —  v~1(p))  - T- 

i  1  -  m  1  J  M  ~  v  M 

with  probability  at  least  p. 

Assume  now  that  A  =  \aa~\  —  aa  E  (^v~1(p),l  —  v~l  (p)'j .  Then  we  can  use  the 
estimate  (j2.13j).  which  holds  unconditionally  with  probability  at  least  8 / 7T2  >  p.  In 
this  case,  we  have 

\a-a\<^A  max{A,  1  -  A}  <  (l  -  ir\p))^. 

This  proves  (12^11). 

We  found  the  estimate  (12.201)  by  numerical  computations.  □ 


Figure  2.1:  The  estimate  (I2.19P  of  C(p)  for  p  G  [0,8/7T2] 

From  Figure  12.11  we  see  that  the  estimate  (12.191)  is  almost  linear  on  the  inter¬ 
val  [4/ 7T2 , 8 / 7T2] ,  which  explains  why  the  right  hand  side  of  the  estimate  (I2.20j)  is 
small. 

We  now  find  a  sharp  bound  on  the  worst-probabilistic  error  of  the  QS  algorithm. 


We  show  that  for  large  M  and  N/M  the  bound  obtained  in  Corollarvl2.3.7lis  optimal 
for  p  G  (1  /2, 8 / 7T2] . 


46 


CHAPTER  2.  QUANTUM  BOOLEAN  SUMMATION 


Figure  2.2:  The  function  v  on  [0, 1].  The  two  horizontal  lines  show  4/7 r2  and  8/7 r2 
levels.  The  part  of  the  graph  between  the  arrows  shows  that  v  is  almost  linear. 


Theorem  2.3.8.  For  large  M  and  N/M ,  the  worst-probabilistic  error  of  the  QS 
algorithm  is  given  by 


e""-rr°(M,p)  =  (l-i’-'W)^(l  +  0(M-‘)  +  0(MiV-1))  /orpe(i,J,]. 
Here,  v  is  as  in  Corollary\2.3.1\  and  1  —  v~1(p)  ~  ^2,P  +  \  by  \2.20\). 

Proof.  From  Corollary  12.3.71  it  is  enough  to  show  a  lower  bound  on  the  error.  Define 


7T 


Si  =  Sill 


\W\ 

M 


For  large  M,  we  have 


and  s  2  =  sin2 


Si  —  \  +  0(M  2) 


'7r(1+  \\M]  ) 


M 


and 


ttx  (n(l  +  2\\M])\  (1  +  0(M-1))tt 

s2  —  Si  =  sm  —  sin  1  1  — 


MJ 


M 


M 


There  exist  two  Boolean  functions  /1  and  /2  with  sums  =  a and  a2  =  a f2  such 
that 

| a*  —  Sj|  <  N~l  for  i  =  1,  2. 
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Since  aSi  =  \\M~\  +  (i  —  1)  and  the  derivative  of  aa  for  a  —  ^  is  M/n,  we  have 
<rai  =  \\M\  +  0{MN-Y)  and  aa2  =  \\M\  +  1  +  0(MN^1). 


Obviously,  a*  =  ki/N  for  some  integers  ki  with  k\  <  fc2.  Consider  ax/N  for  x  G  {ki,  k\  + 
1 , ,k2}.  Then  ax/N  varies  from  crai  for  x  —  k i  to  <ra2  for  x  =  k2.  Since  G 

[|,  |),  for  a  positive  and  small  7]  (we  finally  let  rj  go  to  zero),  we  can  choose  x  =  xv 
such  that  for  a*  =  x^/N  we  have 


■a*  :=  \\M\  +v-\p)  +  rj  +  0(MN-1). 


For  large  N/M,  we  then  have 


[aa*\  =  \l M]  ,  [cra*]  =  \\M\  +  1 


and 


va*  ~  K*J  =  v  1(p)  +  q  +  0(M/N ),  |~cra*]  -  cra*  =  1  -  v  x{p)  -  r]  +  0(M/N). 

Let  denote  the  output  for  the  outcome  [cra*],  and  a 2  for  [aa*\. 

Due  to  (12.7|)  and  (12.91) .  of  Theorem  12.3.21  we  have 

| a*  —  a* |  =  ^  (1  -v~\p)  -rj)  (l  +  0(M~l  +  MN~1)) 

\a*-a*2\  =  yl(v-\p)+r])(l  +  0(M-1  +  MN-1)). 

Let  us  write  1  +  o(l)  for  1  +  Oipj2  +  M_1  +  MN _1).  The  probability  of  a2  is  given  by 
(12.101)  and  is  now  equal  to 

sin2(7r(n~1(p)  +  rj))  ,  ,  y.  =  sin2(7rn^1(p))  +  ttt]  sin(27rz>-1(p))  ,  ,  ^ 

(vr(n-1(p)+'7))2  1  1  ’’  7r2v-1(p)2(l  +  271/v-1(p))  {  { 

Since  p  =  v(v~1(p))  =  sin2  (nv-1  (p))  /  (nv~l  (p))2 ,  the  probability  of  a2  is 


Since  cot(t)  <  1/t  for  t  G  [^7T,  |7t],  we  see  that  the  probability  of  a2  is  slightly  less 
than  p  for  small  rj. 


48 


CHAPTER  2.  QUANTUM  BOOLEAN  SUMMATION 


We  are  ready  to  find  a  lower  bound  on  the  worst-probabilistic  error 
ewor-pro  (M.p)  =  max  min  max  I  at  —  df(j)\ 

of  the  QS  algorithm.  Take  the  function  /  that  corresponds  to  a*.  We  claim  that 
the  error  is  minimized  if  A  =  {|_cra*J  ,  [cxa*]}.  Indeed,  [aa*\  must  belong  to  A  since 
otherwise  /i( A,  f)  <  1—  P/(L°a*J)  =  1—  p+o(l)  <pforp  >  \.  The  probability  of  [aa*\ 
is  slightly  less  than  p,  and  so  the  set  A  must  also  contain  some  other  outcome  j. 
If  j  =  then  the  error  bound  is  roughly  (1  —  u-1(p)  —  r])n/M ,  and  the  sum 

of  the  probabilities  of  the  outputs  for  the  outcome  and  |Va.]  is  always  at 

least  8 / 7T2  >  p.  On  the  other  hand,  if  \cra*~\  does  not  belong  to  the  set  A,  then  any 
other  outcome  j  yields  the  output  sin2(7r j /M).  Since  sin2 (ir(j +1) / M) — sm2(nj /M)  = 
sin(7r/M)  sin(7r(2j  +  1)/M),  the  distribution  of  the  outcomes  around  ^  is  a  mesh 
with  step  size  roughly  n/M.  Hence,  if  j  ^  [chu],  the  error  is  at  least  roughly  (1  + 
v~1(p))n/M  >  7r(l  —  v~1(p))/M.  Thus  the  choice  j  =  \cra*~\  minimizes  the  error 
and  for  r]  tending  to  zero,  the  error  is  roughly  (1  —  v^1  (p))tt / M .  This  completes  the 
proof.  □ 

From  these  results,  it  is  obvious  how  to  guarantee  that  the  error  of  the  QS  algo¬ 
rithm  is  at  most  e  with  probability  at  least  p.  Since  \a  —  a\  <  (l  —  v~1(p))tt/AI  holds 
with  probability  p,  it  is  enough  to  take  M  >  (l  —  ir1(p))'7r/£.  By  Theorem  12.3.81  this 
bound  is  sharp  for  small  £  and  large  eN.  We  have 

Corollary  2.3.9.  Forp  G  (|,  the  algorithm  QS (/,  |"(1  —  n_1(p)) 7r/e])  computes  a 
with  the  error  £  and  probability  at  least  p  with  |~(l  —  u-1(p))7r/e~|  —1  quantum  queries. 
For  small  £  and  large  eN ,  the  estimate  of  the  number  of  quantum  queries  is  sharp. 

2.3.2  Average-probabilistic  error 

In  this  section,  we  study  the  average  performance  of  the  QS  algorithm  with  respect  to 
some  measure  on  the  set  IB at  of  all  Boolean  functions  defined  on  the  set  {0, . . . ,  N  —  1}. 
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We  consider  two  such  measures.  The  first  measure  pi  is  uniformly  distributed  on  the 
set  i.e., 


Pl(/)  =  2"JV  V/G  Mn. 

The  second  measure  p2  is  uniformly  distributed  on  the  set  of  results,  i.e., 


P2C f)  = 


N 


(N  +  1) 


•f  k 

if  Clf  =  — 

J  N 


We  want  to  estimate  average-probabilistic  errors 


eavg-pro(M,p)  =  p ,(/)  min  max  | af  -  af(j)\  for  i  =  1,2. 

/SBjv  '  M  J  ~P  J 

For  the  measures  p,,  the  mean  of  the  random  variable  a/  is  clearly  | .  However, 
their  hrst  (central)  moments  are  very  different.  As  we  shall  see,  the  moment  for 
the  measure  pi  is  small  since  it  is  of  order  N^1^2,  whereas  the  moment  for  measure 
p2  is  roughly  Since  the  first  moments  are  the  same  as  the  error  of  the  constant 
algorithm  CLf(j)  =  we  can  achieve  small  error  of  order  N -1/2  for  the  measure  p! 
without  any  quantum  queries,  while  this  property  does  not  hold  for  the  measure  p2. 

We  first  consider  the  measure  pi.  It  is  natural  to  ask  if  the  average-probabilistic 
error  of  the  QS  algorithm  features  the  dependence  of  N  similar  to  that  of  the  constant 
algorithm  h/(j)  =  We  shall  prove  that  this  is  indeed  the  case  iff  M  is  divisible 
by  4. 

We  compute  the  first  moment  (the  error)  of  the  constant  algorithm,  which  is 


We  do  this  only  for  odd  N,  the  case  of  even  N  being  analogous.  We  have 
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Thus 


N 


=a7-pro(0,  1)  =  2~N 


k= 0 


N 


1  _  k 
2  ~  N 


if  N  is  odd, 


if  N  is  even. 


(2.22) 


Using  Stirling’s  formula 


k\  =  V2nk  eek/12k  for  some  6k  G  [0, 1]  , 


we  can  estimate  any  of  the  binomial  quantities  in  (12. 22ft  by 

1  1 


=l/(12(iV-l)) 


y/2n  yj N  —  1 


proving  that 

>(0,1)  =  5^2 


N 


=avg-pro^ 

Pi 


-N  i 


k= 0 


(N\ 

1  k 

\k) 

2  ~  N 

< 


,1/(12(JV-1)) 


\A2n  y/N  —  1 


1  1 


=  (l  +  o(l)).  (2.23) 


V2^VN 

We  are  ready  to  analyze  the  average-probabilistic  error  of  the  QS  algorithm. 


Theorem  2.3.10.  Assume  that  M  is  divisible  by  4  and  let  p  G  (1/2,  8/7r2] .  Then 
the  average-probabilistic  error  of  the  QS  algorithms  with  respect  to  the  measure  pi 
satisfies 


0 avg-pro / 


5(M,p)  <  min  < - ,  \  — 

P1  v  ,F)  -  '  4  M’  V  2v r 


7T 


1  + 


7T 


31/(12(7V-1)) 


4 M2  vUvWf 


<  -7r(l  +  o(l))  min  <  — , 


M’  y/N 


Proof.  The  estimate  e^rls~w°(M,p)  <  ep'(g_pro(M,  8/7T2)  <  ep°r"pro(M,  8/7T2)  is  obvious 
and,  applying  Corollary  12.3.51  we  get 


avg-pro 

fcPl 


M,-vr2 


<  -  — 
~  4  M 
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As  before  denote  aa  =  (M/7 r)  arcsin  i/a.  Let  a  =  \  +  x.  We  are  interested  in  the 
behavior  of  <Ji/2+x  f°r  \x\  Clearly  cq/2  =  \M.  Let  |x|  <  \.  By  Taylor’s  theorem, 
we  have 


M  M 
Vi/2+x  —  ~r  + 


x 


4  7T  2^(1 -&)& 

and  2 (1  —  £,x)£,x  >  \/l  —  4x2.  Assume  additionally  that 

M  |x|  ^  1 

7T  a/1  -  4x2  “  4’ 

which  is  equivalent  to  assuming  that 


for  e  (|,  4  +  x) 


7T 

-  (16MM44^2)TA' 

Since  M  is  divisible  by  4,  then  =  |M  for  x  >  0,  and  |~<7i/2+z]  =  \M 

for  x  <  0.  This  yields 

]\/[  I  gg  I 

min  {  \al/2+x\  -  a1/2+x,  a1/2+x  -  [^1/2+xJ  }  <  —7 y==  = 2  •  (2.24) 

We  claim  that 

[cTl/2+a:]  “  &l/2+x  £  [0,  |]  U  [|,  1]. 

Indeed,  for  x  <  0  we  have 

and  for  i>0we  have 

h‘/2+4  “ ffl/2+*  =  1  “  V  vy -W  6 11  “  f 11  =  l!’ 1]’ 

as  claimed. 

Let  a  =  4  +  x.  By  the  proof  of  Corollary  I2.3.5[  the  error  of  the  QS  algorithm 
satisfies 

7 T 

|h  —  a|  <  —  min  {  |"ni/2+x]  —  ai/2+xi  ai/2+x  ~  l_cri/2+iJ} 


a  —  a  I  < 


y/l  —  4x2 


and  by  (12. 24j)  we  have 
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We  split  the  sum  that  defines  e!Ns~pro(M,8/7r2)  into  two  sums.  The  first  sum  is 
for  /  G  Bat  for  which  a  —  cif  —  \  +  x  with  |x|  <  7t/(16M2  +  47T2)1/2  and  the  second 
sum  is  for  /  for  which  a  —  af  —  \  +  x  with  |x|  >  7t/(16M2  +  47T2)1/2.  We  estimate 
the  error  of  the  QS  algorithm  by  |a;|/\/l  —  Ax2  and  by  the  worst-case  error  37t/(4M) 
for  the  first  and  second  sums,  respectively.  Hence  we  have 

gavg-pro  f  j^j-  _  f  £\  \af  2 1 


E(  f\  '  I  2 1 

pi (/)  7  = 

f  :  \df— 1/2|<7t/(16M2+47t2)1/2  V  ^  ^(cty  2) 


—  . Pl(/)i  w 

f :  \\cif  —  l/2|>7r/(16M2+47r2)1/2 


Since  a/  =  k/N  for  some  integer  k  G  [0,  TV],  we  have 


eavg-pro  M  < 


E 

k:  |fc/iV-l/2|<7r/(16M2+47r2)1/2 

3  7T 
f  4M 


-AT  W  |fc/iV-|l 

W  A  -  i(k/N  -  i)2 

E 

fc:  |fc/Af-l/2|>7r/(16M2+47r2)1/2 


Since  1  —  4(/c/iV  —  |)2  >  1  —  7t2/(4M2)  >  the  first  sum  can  be  estimated  as 


E 

k:  |fc/V— l/2|<7r/(16A42+47r2)1/2 


1V\  \k/N-\ 

k  )  /l  A  [  I  AT 


1-4 (k/N -ky 


'  k:  |fc/iV— l/2|<7r/(16M2+47r2)1/2 


/iV\ 

1 

k 

UJ 

2  ~~ 

N 

The  second  sum  can  be  estimated  by 


E 

k:  |fc/Ar-l/2|>7r/(16Af2+47r2)1/2 


N\  (16M2  +  4tt2)1/2  k  1 
/c  /  7T  N  2 


<  3  \  / 1  +  V  2" 

V  4M2  ^ 

fc:  |fc/iV— 1/2|>7t/ (16M2+47r2)1/2 

Adding  the  estimates  of  these  two  sums,  we  obtain 


(N\ 

k 

1 

\k) 

N 

~  2 

M  <3  /!+  E2' 


/iV\ 

1 

k 

UJ 

2  ~ 

N 
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The  last  sum  is  given  by  (I2.22H  and  estimated  by  (12.23b  Hence 

1 


3avg-pro 

Pi 


TT 


\l  + 


Z7( 


7T 


4 M2  yiV^T 


1/(12(JV-1)) 


which  completes  the  proof. 


□ 


In  the  next  theorem  we  consider  the  case  when  M  is  not  divisible  by  4. 


Theorem  2.3.11.  Assume  that  M  >  A  is  not  divisible  by  A,  and  let  p  G  (1/2,  8/7 r2]. 
Then  the  average-probabilistic  error  of  the  QS  algorithm  with  respect  to  the  mea¬ 
sure  pi  satisfies 


avg-pro  ( M  )  >  (  1  _  JL  _  I  )  M  -  2  exp 

P1  V  ~  AM  V  M  pj\ 

Proof.  Let  M  =  AM'  +  r  for  r  G  {1,  2,  3}.  Let  cra  = 
of  Theorem  12.3.101  for  \x\  <  |  we  have 


N  7T2 


\/(3  >  1. 


(8/3M)2/ 

(M/tt)  arcsin  ^/a.  As  in  the  proof 


M  M 


x 


&1/2+X 


4  vr  2^/(1 -&)& 

and  2a/(1  —  £x)£x  >  y/A  —  Ax2.  Assume  additionally  that 


for  ^G(4,|+x) 


M  HI 


1 

<  — . 


TT  y/\  -  Ax2  4/3’ 
which  is  equivalent  to  assuming  that 


^  -  (16M2/32  +  47T2)1/2 '  (2'25^ 

Tims  for  x  satisfying  (12.251).  we  have 

, r  M  x9(x)  .  , 

Vi/2+x  =  M  +  -  + -  with  0(x)  G  [0, 1], 

4  7T  HI  —  4x2 

and  [ai/2+xJ  =  M'  and  |~ai/2+a:~|  =  M'  +  1. 

From  the  proof  of  Corollary  12.3.31  we  have  /x({|"cra/]  ,  |_cra/J},/)  >  8 / 7r2.  Since 
/i(A,  f)>p>\  then  either  |~cra/~|  G  A  or  [craf\  G  A.  We  then  estimate 

eavg-pro  fM,  %  \  >  ^  Pl(/)nhn{|a/-a/(  |_<7a/J  )|,  \af-df(  \aaf]  )|}.  (2.26) 

V  2  f&N 
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We  now  estimate  the  error  of  the  QS  algorithm  for  /  G  IB at  such  that  a/  =  |  +  x 
for  x  satisfying  (12.251)  and  the  outcome  j  =  |"<ra/]  =  M'  or  j  =  [era/J  =  M'  +  1. 
Denote  the  outcome  by  M'  +  k  for  k  G  {0, 1}.  By  Taylor’s  theorem  we  have 


sm 


M'  +  K 
M 


TT  =  Sill 


7  +  -  lr) )  =  5  +  sin(2^)T)^:(«  -  \t 


TT 


M' 


71 


M 


f°r  £k,t  G  [^7t,  f-7r  +  (7 t/M)(k  —  \t)\.  Since  sin(t)  >  2t/ir  for  t  G  [0,7t/2],  we 
have  |  sin(2£K)T)|  >  1  —  |4/c  —  r\/M.  Consider  the  error  for  the  outcome  M'  +  k 
and  x  satisfying  (12.251).  Then  |x|  <  tt/(4(3M )  and  the  error  can  be  estimated  by 


-  +  x  —  sm 


M'  +  k 


M 


71 


7 r 

x-sin(2^iT)  —  (re-  ±t) 


7T  -  In  //-  \i  11  n 

>  T7t\k~  atW  sin(2^KjT)|  -\x\>  - 


TT  1 4/c  —  r| 


M 


4  M 


1  - 


|4/c  —  r| 
M 


Clearly,  |4/c  —  r|  G  {1,2,3}  and  \Ak  —  t\/M  G  [1/M,  3/M].  Then 


TT 

4/3  M' 


|4/c  -  r|(l  -  |4/c  -  t\/M)/M  >  (1  -  1/M)/M. 


Therefore 

7T  ( \  (  1  \  1  \  _  7T  /  1  1 

-MV4V1_My_4^y  “  4M  V1  “  M  “  ^ 
Hence,  for  /  such  that  a/  =  \  +  x  with  x  satisfying  (12.251)  we  have 

min  (| af  -  af(  [a(lf\  )  |,  \af  -  af(  )  I }  ^  ~  •  (2-27) 

We  are  now  ready  to  estimate  e^rlg'pro(M,p).  First,  by  (12.261).  we  have 

e»«”(M,p)  >  £  Pi(/) 

/:  |a/-l/2|<(16M2/32+47r2)1/2 

x  min  { |a/  —  a/(  [o’a/J  )|,  |a/  -  a/(  f^/l  )  |}- 


1 

2+X 


sm 


M'  +  /c 

M 


■  TT 
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This,  (12.2  71).  and  the  Bernstein  inequality,  ^2k.  \k/N-i/2\>e  2  JV(’^)  <  2e  JVe2/4;  yield 


eavg-pro(M)p)  > 


> 

> 


E 


Pi  (/) 


/:  | ay  —  l/2|<7r(16M2/32+47r2)1/2 

7 r 


7T 

4M 


1  -E  -  - 

M  p 


AM 

71 

AM 

71 

AM 


1  1 

1  ~~  M  ~  P 

1  1 
1~M~p 

l«  Jl  _  i 

M  p 


E 

k:  Ifc/Af— l/2|<7r(16M2/32+47r2)1/2 

Ntt2 


-N 


1  —  2  exp 
1  —  2  exp 


4(16M2/?2  -  4) 
Nn2 


(8pMp 


which  completes  the  proof. 


□ 


Obviously,  in  the  average-probabilistic  setting,  we  should  use  the  QS  algorithm 
with  M  divisible  by  4.  Then  Theorem  12.3.101  states  that  the  error  is  of  order 
min{M-1,  TV-1/2}.  Papageorgiou  |36]  proved  that  for  any  quantum  algorithm  that 
uses  M  quantum  queries  the  error  is  bounded  from  below  by 
cmin{M-1,  IV-1/2}  with  probability  p  £  (1/2,  8/ 7r2] .  Here,  c  is  a  positive  number 
independent  of  M  and  N.  Hence,  the  QS  algorithm  is  optimal  also  in  the  average- 
probabilistic  setting  for  the  measure  pi  as  long  as  we  use  it  with  M  divisible  by  4. 

We  now  turn  to  the  measure  p2-  Clearly,  the  average-probabilistic  error  of  the  QS 
algorithm  is  bounded  by  its  worst-probabilistic  error,  which  is  of  order  M^1  with  prob¬ 
ability  p  £  (1/2,  8 / 7T2] .  It  turns  out,  again  due  to  a  recent  result  of  Papageorgiou  |3B] 
that  this  bound  is  the  best  possible,  since  any  quantum  algorithm  that  uses  M  quan¬ 
tum  queries  must  have  an  error  proportional  at  least  to  M~l .  Hence,  the  factor  IV-1/2 
that  is  present  for  the  measure  p!  does  not  appear  for  the  measure  P2,  and  the  behav¬ 
ior  of  the  QS  algorithm  is  roughly  the  same  in  the  worst-  and  average-probabilistic 
settings. 
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2.3.3  Worst-average  error 

The  QS  algorithm  uses  M  —  1  quantum  queries.  In  this  section,  we  restrict  our 
consideration  to  the  only  interesting  case,  namely,  when  M  is  much  smaller  than  N. 
We  denote 

Sa,f  min  ^  [a*,]  ® af  1  > 

where,  as  before,  craf  =  (M/tt)  arcsin  y/a.  Clearly,  saf  G  [0,  |]  and  saf  =  0  iff  aaj  is  an 
integer.  We  shall  usually  drop  the  subscript  /  and  denote  aa  =  <raf,  sa  =  saf  when  / 
is  clear  from  the  context.  As  in  Section fl.2.1. 11  let  //(•,  /)  denote  the  measure  on  the 
set  {0, . . . ,  M  —  1}  of  all  possible  outcomes  of  the  QS  algorithm.  We  also  define 

Am  —  |  sin2  '■  j  =  0, 1, . . . ,  M  -  l| 

as  the  set  of  all  possible  outputs  of  the  QS  algorithm  with  M  —  1  queries.  Let 

pf(a)  =  /r^j  j  G  {0, 1,...,M-  1}  :  sin2  =  «  J , /)  Vo  G  Am, 

denote  the  probability  of  the  output  a.  Note  that  a  =  sin2 (7 rj/M)  =  sin2(7r(M  — 
j)/M).  Hence,  from  (12.61).  we  clearly  see  that  p/(a)  =  Pf(j)  +  p/(M  —  j)  for  j  ^ 
0,  M/2. 

For  q  G  [1,  00),  we  first  analyze  the  local  average  error 

/  M~l  \l!q  f  \ 

<vs(/,Af)  :=  (  \af -«/(i)l9)  =  (  Pf(a)\af  ~  a\q)  »  (2-28) 

'  j= 0  2  '  a£AM  ' 

for  a  fixed  function  /  G  Bjv- 

2. 3. 3.1  Local  average  error 

The  local  average  error  (12.281)  for  q  >  1  is  estimated  in  the  following  theorem. 


Theorem  2.3.12.  Let  q  G  (l,oo).  Denote  a  =  af.  If  a  a  £  Z,  then  e™9(f,M)  =  0. 
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If  a a  (f  Z,  then 


eqV9(f,  M)q  — 


Sin  7T5n 


f*7T— 7TSa  /M 


M  7T 


'  irsa/M 


sin(x)q  2|  sin(x  +  2#a)|c,Gte 


< 


,n  ,  ^  1  sin(7TSo)  ,  sin2(7rsa) 

(!  +  2(1  -  W) - Wq - +  — gj— 

with  sa  =  [a a J  -  cra  and  sa  =  cra  -  [era] 


2(1  -  <59i2)  +  g  /  sin9  2(a;)<ir  J,  (2.29) 


Proof.  If  <7a  E  Z  then  it  was  shown  in  Section  12.3.21  that  there  exists  a  E  *4m  such 
that  a  =  af  and  p/(a)  =  5ai0/  for  all  a  E  .Am-  Then  ej)vg(/,  M)  =  0  as  claimed. 

Assume  that  cra  ^  Z.  Using  the  closed  form  of  of  p / (j )  as  in  (12. 6[),  we  rewrite  (|2.28l) 
as 


M- 1  .  2 


«'(/,  M)Y=Y 


3=0 


sin  2  {M6  a) 
2  M2 


Sill 


+ 


~  g-q) 

M 


q- 2 


Sill 


7T(j  +  da) 
M 


sm 


9-2 


7r(j  +  da) 
M 


Sill 


7T(j  ~  da) 
M 


We  have 

M— 1 


E 

3=0 


Sill 


7T(j  +  0-q) 
M 


9-2 


Sill 


7T(j  -  0-q) 
M 


M 

E 

9=1 


Sill 


?r(M  -  j  +  0-q) 

M 


9-2 


Sill 


?r(M  -  j  -  crQ) 
M 


Using  the  7r-periodicity  of  |  sinx|,  we  see  that  the  last  sum  is  equal  to 


9-2 


sin 


(  n(j  +  <ra)\ 

\  M  J 


9 


da) 


9-2 


sin 


/  7T(j  +da)\ 

V  M  J 


1Note  that  the  last  integral  is  finite.  This  is  obvious  for  q  >  2.  For  q  £  (1,  2),  the  only  singularities 

are  at  the  boundary  points  and  are  of  the  form  xq~2  for  x  approaching  0.  The  function  xq~ 2  is 
integrable  since  q  >  1. 
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Therefore 


with 

M— 1 

Sm,,  = 

3=0 


eoVg(/)  M)q  — 


sm 


n(j  ~ 

M 


q-2 


Sill 


sin2  (M  9  a 
M2 


7T(j  +  Va) 


-S 


M,qi 


(2.30) 


M 


We  split  SM,q  as 


&M,q  —  $M,q  ~ 


Sill 


7T  l_CTaJ 

M 


M—l 

£ 

3=0 


-Or 


■l  71 3  a 
sm  —  —  9  a 

M 


q-2 


si"  I  ft  +  S‘ 


q-2 


Sill 


7T  l_CTaJ 

M 


+  9  a 


Observe  that  (it /M)S'Mq  is  the  rectangle  formula  for  approximating  the  integral 
f  |  sin(x  «  9a)\q  2 1  sin (x  +  9a)\q  dx. 

J  [0,7r]\[7T  [(Ta} /M,n\(Ta\/M] 

The  error  of  the  rectangle  quadrature  for  k  €  N  and  an  absolutely  continuous 
function  /  :  [a,  b]  — >  M  whose  first  derivative  belongs  to  Li([a,b\)  satisfies 

fc-i  /  T  \ 


f(x)dx-b-1-^'^Tf(a+jb  a 


k 


3=0 
19-2  | 


k 


< 


b  —  a 


k 


I  f'(x)  \  dx. 


(2.31) 


Defining  h(x)  =  |  sin(a;  —  9a)  |y  “|  siii(a;  +  #0)|9  and  Da  =  [0, 7r]  \  [it  \_<ja\  /M,tt  |~cra]  /M] 
and  using  the  error  formula  02.311)  for  the  subintervals  [0, 7 r  |_<raJ  /M)  and  (7r  [cxa]  /M,  7r], 
we  get 

7T 

h(x)  dx 


—S' 
M  M’q 


1 Da 


-M'  \h'(x)\dx- 


1 Da 


Define  H(x)  =  h(x+9a)  =  |sin(a;)|g  2|  sin(;r+2#a)|g  and  Aa  =  [— 9a,  n  —  #a]\[7r(|_craj  — 
cra)/M,7r(aa  -  \aa])/M]  .  We  have 

/  h(x)  dx=  H(x)  dx, 


\h'(x)\  dx=  \H'(x)\  dx. 


' Da 


>Da 


By  the  7r-periodicity  of  the  integrand  H  we  have 


H(x)  dx  = 


da 


H(x)  dx  + 

>''K—'KSa/M 

T-da 


'  irSa/M 


H(x)  dx 


rrr-8a 


H(x )  dx  +  H(x)  dx  = 

J  nSa/M 


i>ir—irsa/M 


I  iTSa/M 


H(x)  dx. 
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Analogously, 


C  pTT—TTSa/M 

/  \H\x)\  dx  =  "  \H\x)\dx. 

'  Aa  Jnsa/M 


For  x  G  [irsa/M,  n  —  nsa/M)  the  sine  is  positive  and 


\H'(x)\  <  \q  —  2|  sin9  3(a;)|  cos(x)|  +  gsin9  2(x). 


It  is  easy  to  check  that  for  q  ^  2  we  have 


^7T—7 TSa/M 


'  nSa/M 


\q  —  2|sin9  3(a;)|  cos(x)|  dx 


/  /M2 

—  \q  —  2|l  /  sin9_3(a;)  dsin(a;) 

V  J irSa/M 

\q~  2 


r*7T— 7TSa/M 


'  7r/2 


sin9  3(x)dsin(x) 


9-2 


/  VTSa 

I  M 


2  -  sin9"2  (2AM  sina-2  ,  „  _ 


M 


From  this  we  get 


r»7T  —  7ISa/M 


'  TTSa/M 


\H\x)\  dx  <  (1-<J,,2)  2  +  sin92(^ 


M 


+  sin9  2  ^  ‘  *'1  J  j  +q  j  sin9  2{x)dx. 


We  then  finally  get 


7T 

MbM« 


fir— TTSa/M 


f7TSa/M 


H(x)  dx 


<-  (i-M  2  +  sin-  ^ 


7T  /  „  .  q_2  (  KSn 


f 


+  sin9"2  (  ^  )  )  +  sin9"2  (2AM  +  g  /  sin9"2  (a;)  da: 


V  M 


A 


7TS, 


V  M 


sm  7TS„  =  sm (ns„)  =  sin  7rs„ 


Observe  also  that 
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Since  sin (x)/[M  sin (x/M)]  <  1  for  x  G  (0, 7r] ,  we  get 


vr  sin(? rsa)  ■  i  \  T  K§a/M  zr <  ^  a 

- TT - SM,q  -  sm(7rsa)  /  H(x)  dx 

M  J-KSa/M 


<  -d  -  W  (  (  ^  )  +  «dn^  (^)  +  -in9- 


7rsin(7rsaJ 

M 


2(1  -  5gt2)  +  q  I  sin9  2(x)dx). 


Using  sin(7rsa/M)  <  tt/M  we  obtain 


vr  sin(7rsa)  .  ,  ^  r  n-jM  ZJ,  ,  j 

- 77 - SM,q  ~  sm(7rsa)  /  H(x)  dx 

M  Jnsa/M 


<  (1  +  2(1  -  +  nSm^Sa'  (  2(1  -  ^,2)  +  q  shP-\x)  dx  ] . 


M 


Finally,  since  sin 2{M6a)  =  sin2(7rsa),  we  complete  the  proof  by  using  the  estimate 


of  SM,q  in  (IPOD. 

Theorem  12.3.121  implies  the  following  corollary. 

Corollary  2.3.13.  Let  q  G  (l,oo).  If  aa  G  Z,  then  e<P9{f^M)  =  0.  If  cra  ^  Z,  then 


□ 


Sill  (7TSa) 


7T 


sin9  2(;r)|  sin(:r  +  26,a)|f,ch 


O 


sin(7rga) 


1  1/9 


,  (2.32) 


with  sa  G  (0,  |];  where  the  factor  in  the  big  O  notation  is  independent  of  f  from  Bat, 
and  also  independent  of  N . 

We  now  consider  the  case  q  =  1  and  present  estimates  of  efve(f,  M )  in  the  following 
lemma. 


Lemma  2.3.14.  Let  a  =  aj.  If  aa  G  Z,  then  e“U£,(/,  M)  =  0.  If  aa  ^  Z,  then 


er  (/,  M) 


sm2(nsa)sm(26a) 

M 


JM,a 


sin2(7Tsa) 

<  - Yl - |cos(20o)|, 


(2.33) 
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where  sa  G  (0,  and 


M- 1 


M  22 


3= 0 


cot 


n(J  +  gg) 

M 


Proof.  The  case  oa  G  Z  can  be  proved  as  in  Theorem  12.3.121  Assume  that  cra  ^  Z. 
Using  (I2.6jh  we  get 


M-l  .  2 


=rs(/.  m)  =  £ 


3=0 


sin \M0a 
2  M2 


( 

sin(7r  (j  +  aa)/M ) 

sin(7r(j  -  <Ja)/M ) 

l 

sin(vr(j  -  cra)/M ) 

1 

sin(vr(j  +  cra)/M) 

As  in  the  proof  of  Theorem  12.3.121  we  conclude  that 

sin2  (M61 


ers(/,M)  = 


M2 


*Sm,i, 


where 


sin(7r  (j  +  <ra)/M) 

M-l 

-  v 

sin(vr (j  -  fcra]  +  sa)/M  +  2 Ba) 

sin(7r(j  -  cra)/M ) 

3=0 

sin(vr(j  -  \aa]  +  sa)/M) 

M-l 

Sm,i  =  y 

3=0 

with  sa  =  \<Ja~\  ~  o~a-  Changing  the  index  j  in  the  second  sum  to  j  —  [cra],  and  using 
periodicity  of  the  sine,  we  get 

M-l 


Sm,i  —  22 

3=0 


sin(7r(j  +  sa)/M  +  2  6a) 


sin(7r(j  +  sa)/M) 


and  consequently 


M-l 


Sm,i  —  22 

3=0 


cos(2 9 a)  +  sin(20o)  cot 


7T(j  +  Sg) 

M 


Using  the  triangle  inequality  twice,  we  obtain 

M-l 


SM,i  -  sin(26»a)  22 
j=o 


cot 


7T(j  +  Sg) 

M 


<  M |  cos(20a)|. 


Let  sa  =  ua  —  [cra\.  Observe  that  sa  —  1  —  sa.  Since  the  cotangent  is  7r-periodic  and 
the  function  |  cot(7r(-)/M)|  is  even,  we  get 


£ 

3=0 


cot 


7T(j  +  Sg) 

M 


£ 

3=0 


cot 


Kti  +  Sg) 

M 


=  ME 


Ad, a- 
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This  and 


sin 2  (MO  a)  =  sin2(7T(ja)  =  sin2(7rsa) 


yield  (12.331).  as  claimed. 


□ 


From  Lemma  12.3. 141  we  see  that  the  sum  Em, a  is  the  most  important  part  of  the  local 
average  error  e^vg(M,  /).  We  now  estimate  Em, a- 

Lemma  2.3.15.  Assume  that  cra  (f  Z  and  M  >  3.  Then 


E  M,a  ~  —  COt  ( 


(  TTSc 


1 

M 


cot 


M  \M 

rir(M— l+sa)/M 
^  J 7r(l  +  Sa)/M 


7t(M  -  1  +  Sa) 


M 


cot  x\  dx 


< 


’•n(M—l+sa)/M 


7 tM 


'  ir(l+sa)/M 


sin2  x 


dx.  (2.34) 


Proof.  This  can  be  shown  by  using  the  error  formula  for  rectangle  quadratures  (12.31  ft . 
Note  that  7tE M,a~  (^ / M)  cot(nsa/M)  —  (n/M)\  cot(7r(M  —  l  +  sa)/M)\  is  the  rectangle 
quadrature  for  the  integral  J7(i+s] um^M  I  cot  x\  dx  with  k  —  M  —  2  >  1.  We  then 
obtain  (j2.34jl  by  using  (I2.3ip.  □ 

We  now  present  the  final  estimate  on  the  local  average  error  e^vg(/,  M). 

Theorem  2.3.16.  Assume  that  /  £  IB at  ond  a  =  aj.  For  M  >  3,  the  average  error 
of  the  QS  algorithm  for  the  function  f  satisfies 
2  sin2(7rsa)  sin(2#a)  In  M 


e7V,  M) 


3tt  T  2  T  ln(7r2) 

<  - -  sin(7rs0).  (2.35) 


7 r  M 

Proof.  For  cxa  €  Z  we  have  sa  =  0,  so  that  (I2.35P  holds  since  e^vg(/,  M)  =  0  by 
Assume  that  ua  (f  Z.  From  Lemmas  12.3. 141  and  12.3. 15l  we  have 


eH/,  M)  -  ,  cot  x|  dx 


7 tM 


< 


Sin  f 7T sa ) 
M 


sin(20, 

M 


'ir(l+sa)/M 

”"cot(^,+ 


cot 


7 T(M  -  1  +  Sa) 
M 


^  /•7r(M-l+sa)/M  2 

—  /  — o —  dx  I  +  |  cos(20a)| 

7T  Jn(l+sa)/M  Sill  X 
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Observe  that 


cot  d  dx  =  In  I  sin-1  ( sin  1  Sa) 


0(1  +sa)/M 

( 7 r(M  -  1  +  s0) 

cot  - — - 

V  M 

rn{M—\+sa)/M  i 


M 


,  7r(l  —  Sa)  \  (  VTSq 

=  COt  M  )  “  COt  (  ~M  )  ’ 


M 


'7r(l  +  Sa)/M  Shl  X 


dx  =  coti’r(1~S<>)i+cot^(1  +  S“) 


<  2  cot 


M 

VTSa 

M 


M 


The  four  formulas  above  yield 


e?\f,  M)  -  sin2(jr»«)sill(2f)«)  hl  (  sin-i  f yiyMi  Sin-1  0(!  -  >») 


7 rM  \  \  M  j  V  M 

sin2(7rsa)  ^  (2  +  2 / 7r)  sin(2#a)  ^7rs 


71  r  1  71  r  COt  77  +  I  COS(26b)|  . 

MV  M  \  M  J  v  ;7 

Observe  that  sin(7rs0)/[M sin(7rsa/M)]  <  1  since  sa  G  (0,  |].  This  and  the  obvious 
estimates  of  sine  and  cosine  yield 


avg ,,  ti  sin2(7rsa)  sin(20a)  /  .  _x  f  n(l  +  sa)  \  .  f  n(l  -  sa) 
e1  &(f,M) - rv - m  |  sm  (  - rv -  I  sm  1 


7 tM 


M 


<  3+  - 


M 

2  \  sin(7rsa) 


(2.36) 


7T  J  M 

Consider  now  the  left  hand  side  of  (12.361).  Remembering  that  M  >  3,  and  since  2x/n  < 
sin  x  <  x  for  x  G  [0,  7t/2],  we  get 


.  _i  (tt(1  +  Sa)\  .  _i  f  7(1  —  Sa) 
111  Sill  - -  Sill  1 


M 


M 


2  In  M 


Thus  by  (12.361)  and  (12.37)).  we  get  the  final  estimate  (12.351). 


<  ln(7r2).  (2.37) 


a 


From  Corollary  12.3.121  and  Theorem  12.3.161  we  get  sharp  estimates  on  the  worst- 
average  error  of  the  QS  algorithm. 

Theorem  2.3.17.  Let  M  >  3.  Then  the  worst-average  error  of  the  QS  algorithm 
satisfies  the  following  bounds. 
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For  q  e  (1,  oo), 


0wor-avg 


(M)< 


1  fl 


\  +9 


M1/^  \  7T 

The  last  estimate  is  sharp,  i.e. 


|  —  /  sm9  2(x)dx\  (l  +  o(l)). 


(2.38) 


'0 


e~vg(M)  =  0 


In  particular,  for  M  —  2  divisible  by  4  we  have 


(2.39) 


0wor-avg 


l  /i 


\  1/9 


(M)>^I^f-/  sin9  2  (a:)  |  cos(:r)|9  da:  )  (l  +  o(l)),  (2.40) 


and  the  ratio  of  the  integrals  in  \2.SSi)  and  \2.fO\)  are  approximately  1  for  q 
close  to  1. 


For  q  =  1, 

e”5(M)  < 

This  estimate  is  sharp,  i.e., 


2  In  M  3n  +  2  +  ln(7r2) 

7T  M  M  7T 


emM)  =  @(M_1  InM). 


(2.41) 


(2.42) 


In  particular,  for  M  —  2  divisible  by  4  we  have 

ewor-avg(M)  >  2  In  M  _  3tt  +  2  +  ln(7T2) 

1  —  77  M  Mtt 

Proof.  Consider  first  the  case  g  e  (1,  oo).  By  Corollary  12.3. 131  we  have 

1/1  rn  \1N 

egVg(/,M)  <  Jo  smq-2(x)dxj  (l  +  o(l))  V/  G  By, 

where  o(l)  is  independent  of  /.  This  yields  f)2.38p. 

The  estimate  (12. 38ft  is  sharp  since  if  we  let  /  be  a  Boolean  function  such  that 
Saf  ~  then  (12.321)  yields  (12.391).  In  particular,  for  M  =  4k  +  2  and  aj  =  1/2  we 
have  9af  =  vr/4,  craf  —  M/4  —  k  +  1/2  and  sa  =  /.  Therefore 

eqVg(f>M)  —  f  smq-2(x)\cos(x)\qdx\  (l  +  o(l)), 
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which  proves  (12.40b  For  q  close  to  1,  the  value  of  [J  sin q~2(x)dx  is  mostly  due  to  the 
values  of  sin9~2(;r)  close  to  0  and  it.  Since  |  cos(a;)|'3  is  then  approximately  equal  to  1, 
the  ratio  of  the  upper  and  lower  bound  integrals  is  about  1. 

For  q  —  1,  the  estimate  (12.411)  follows  directly  from  Theoreml2.3.16l  To  prove  (12.421) 
it  is  enough  to  choose  a  Boolean  /  for  which  the  numbers 

sin2(7rsa/)sin(26>a/)  =  sin2 {M9af)  sin (20a/) 

are  uniformly  (in  M)  separated  from  0,  see  Theorem  12.3.161  More  precisely,  since  af 
can  take  any  value  k/N  for  k  =  0,1,...,  N,  we  take  a  Boolean  function  /  such 
that  |  aj  —  sin2(7r/4  +  7t/(5M))|  <  1/(2 N).  For  sufficiently  large  N,  we  have  0Of  ~ 
7t/4  +  7r / (5 M).  For  large  M  —  4A;  +  (3  with  (3  G  {0, 1,  2,  3},  we  then  have 

sin2(M#a/)  sin(2#a/)  «  sin2  sin  >  c  >  °’ 

for  some  c  independent  of  M. 

In  particular,  for  M  —  2  divisible  by  4  we  take  N  >  M  and  a  Boolean  function  /  e 
Bat  with  aj  =  1/2.  Then 

sa  —  \  and  sin2 (7rsa)  sin(20a)  =  1, 

which  leads  the  last  estimate  of  Theorem  12.3.171  □ 

2. 3. 3. 2  Quantum  summation  algorithm  with  repetitions 

The  success  probability  of  the  QS  algorithm  is  increased  by  repeating  it  a  number 
of  times  and  taking  the  median  of  the  outputs  as  the  final  output,  see  e.g.,  |T5].  We 
show  in  this  section  that  this  procedure  improves  the  worst-average  error  estimate. 

We  perform  2 n  +  1  repetitions  of  the  QS  algorithm  for  some  n  G  {0, 1, . . .}.  We 
obtain  sin2(7rj1/M),  sin2(7rj2/M), . . . ,  sin2(7rj2n+i/M)  and  let  anj  be  the  median  of 
the  outputs  obtained,  i.e.,  the  (n  +  l)st  number  in  the  ordered  sequence.  Let,  as 
before,  Am  =  {sin2(7r j/M)  :  j  =  0, 1, . . . ,  M  —  1}.  For  a  G  Am,  let  Pnj(&)  he 
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the  probability  that  the  median  anj  is  equal  to  a.  This  probability  depends  on  the 
distribution  function  Fj  of  the  original  outputs  from  Am,  which  is  defined  as 


Ff(a) 


I  En'^My<«/)/K)  for  a  >  0, 

I  0  for  a  =  0. 


(2.43) 


It  is  known,  see  [41j  p.  410,  that  the  distribution  of  the  median  anj  is  of  the  form 

'2n\  fFfO)+pfO) 

I Ff  (a) 

We  are  now  ready  to  estimate  the  worst-average  error 

1/9 


/2rA  rFfO+pfO) 

Pnj(oi)  =  (2 n  +  1)  (  )  /  tn(  1  -  t)n  dt,  Va  G  AM- 

\n  J  JFfipt) 


e“vg  (M)  =  suP  max/"  ^  pnJ(a)\af  -  a\A  ,  q  e  [. l,oo) 


of  the  QS  algorithm  with  2 n  +  1  repetitions. 

We  estimate  e"°r_avg(M)  by  using  Theorem  12  of  jl]  which  states  that  the  QS 
algorithm  with  M  queries  computes  df  such  that 

—  i  iiii  C2 

\af  —  a/ 1  >  ci^  with  probability  at  most  — 


for  any  positive  integer  k,  Here  c\  and  C2  are  absolute  constants  and  /  is  any  Boolean 
function  from  Bat.  If 

k 

I  af  ~  ®n,f\  —  Cl  ’ 

then  for  at  least  n  outcomes  o/(ji),  ■  ■  ■ ,  h/(jn)  we  must  have 

k 

I  af  ~  df(ji)  |  >  ci—  for  l  =  1, . . .  ,n. 


But  the  probability  that  this  occurs  is  bounded  by 

(*: ')  ©  ■  ■ 


Prob  {  | a /  —  anj|  >  C\k/M}  <  ck  n, 


It  follows  then  that 
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where  c  which  may  depend  on  n.  We  now  use  the  standard  summation  by  parts. 
Define 

Pk  =  Prob  {  ci(k  —  1)/M  <  \af  —  anj\  <  Cik/M}. 

Then,  by  the  estimate  above  we  get, 


Therefore 


J^Pk  <  cl'"  VZeN. 

k>l 


oo  oo  k 

e™-,vg(M),  <  J2p t(Clk/M)"  =  (d/M)’  £pt  £(I«  -  (i  -  1)«) 
k= 1  k=  1  1=1 

oo  oo  oo 

=  (ci /M)q  ~(l~  !)9)  -  cM_,E|rl"n  <  cM~q 

1=1  k=l  1=1 

for  n  >  q,  with  the  number  c  =  cq^n  depending  only  on  q  and  n.  In  fact,  taking  n  = 
|~g~|  +  1  it  is  easy  to  check  that  cq^n  is  a  single  exponential  function  of  q.  Hence,  by 
taking  the  gth  root  we  have 


wor-avg 

eq,n 


(M)  < 


-1 


with  c 


q,n 


of  order  1.  Therefore  we  have  proved  the  following  theorem. 


Theorem  2.3.18.  The  worst-average  error  of  the  median  of2(\q]  +1)  +  1  repetitions 
of  the  QS  algorithm  with  M  quantum  queries  satisfies 


=  0{M-% 

with  an  absolute  constant  in  the  big  O  notation  independent  of  q  and  M . 

The  essence  of  Theorem  12.3.181  is  that  the  number  of  repetitions  of  the  QS  al¬ 
gorithm  is  independent  of  M  and  depends  only  linearly  on  q.  Still,  it  permits  us  to 
dramatically  improve  the  worst-average  error  of  the  QS  algorithm.  As  we  already 
mentioned  in  the  introduction,  the  bound  of  order  AW1  is  a  lower  bound  on  the  worst- 
average  error  of  any  quantum  algorithm.  Hence,  the  QS  algorithm  with  repetitions 
is  also  optimal  in  the  worst-average  setting. 
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2.4  Simulation 

In  this  section  we  show  how  the  quantum  Boolean  summation  algorithm  for  the 
class  IB tv  can  be  simulated  on  a  classical  computer.  From  the  previous  sections  we 
see  that  the  core  problem  is  the  computation  of  amplitudes  of  the  final  state  since 
they  are  needed  to  compute  the  probabilities  of  all  possible  outcomes.  We  present 
a  MATLAB  procedure  QSsimul  that  computes  all  amplitudes  of  the  final  state  before 
the  measurement.  The  cost  of  QSsimul  is  of  order  (N/e)  log2(l/e). 


2.4.1  QS  algorithm  simulation 


Our  simulation  computes  all  amplitudes  of  the  final  state  of  the  QS  algorithm  for  a 
Boolean  function  /  e  Btv  and  given  error  e  6  (0, 1). 

First,  based  on  the  estimate  (12. 171) .  we  compute  the  number  m  of  qubits  required 
to  provide  a  desired  accuracy  with  probability  at  least  8/7 r2.  During  the  computation 
we  use  an  N  x  M  matrix  fc=0~'’  where  M  =  2m,  with  MN  coefficients  of  the 

quantum  states  1 772)  and  (773)  as  defined  in  Section  [2721  The  coefficients  of  the 
states  |  r]i)  are  defined  with  respect  to  the  computational  basis 

Clearly  for  i  —  1  we  have  ajk  =  (MTV)-1/2.  To  obtain  the  state  (772)  we  need  to 
apply  the  Grover  iterate  operator  A m(Qf)-  We  now  show  how  to  do  this  efficiently. 
We  have 

M—l  N—l  M—l  N- 1 

m = Am(<3)ivm+n|o)  =  -==Am(Q)  y  ij>  E  w  =  -/m  E  b’>  &  E  1*1 

v  j=0  k= 0  V  j= 0  k= 0 


Observe  that 

JV— 1  N—l  N- 1 

Q  a ^  =  Y  -akWnS0WnSf\k)  =  Y  -(-l)f(k)akWn(I  -  2\0) (0\)Wn\k) 


k= 0 


k= 0 
N—l 


k= 0 
N- 1 


E  — (— k)  +  -mYa“(--1'imwM 

k= 0  »  k= 0 

)  | k). 


TV-1  /  N—l 
k= 0  \  j= 0 


k= 0 

fU)  _  (_i  \f(k ), 
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Thus,  the  coefficients  ajk  of  the  state  1 772)  satisfy  the  following  recursive  formula 


(. MN ) 


-1/2 


EiV  - 

v=< 


N-l  „2 


N 


p= 0 


■l)/(p)  -  (-1  )/(fc)aJ_llfc, 


i  =  0, 

j  =  1, M  —  1, 


(2.44) 


for  k  —  0, . . . ,  N  —  1. 

To  obtain  the  state  lr/3)  we  need  to  apply  the  inverse  quantum  Fourier  trans- 


form  <8)  I.  Hence 

M— 1  N—l 

b>  =  {f-1  alk\j)\k )  - 

j= 0  k= 0 

M—l  N—l  M—l 

*  j= 0  k= 0  p= 0 

which  yields  the  final  form  of  the  coefficients  as 


M—l 

p= 0 

We  see  that  the  rows  of  the  matrix  [aj  k]  are  discrete  Fourier  transforms  of  the  re¬ 
spective  rows  of  the  matrix  [a^k]. 

It  is  easy  to  see  that  the  costs  of  the  steps  computing  the  states  \rji)  and  (772)  are 
of  order  MN.  We  can  compute  the  state  (773)  by  using  the  Fast  Fourier  Transform 
algorithm  with  the  cost  of  order  MNm  =  AIN  log2  M.  We  see  then  that  the  total  cost 
of  our  simulation  is  of  order  of  the  cost  of  computing  the  state  (773) .  Assuming  m  = 
[log2(7r/e)]),  see  Corollary  12.3.91  the  total  cost  is  of  order  (N/e)  log2e_1. 

The  simulation  algorithm  outlined  above  has  been  coded  in  GNU  Octave,  see  [10] , 
and  can  be  easily  ported  to  matlab. 

The  program  QSsimul  computes  all  the  coefficients  of  the  final  state  of  the  QS 
algorithm,  stores  numerical  results  and  creates  a  PostScript  file  containing  graphs  of 
the  exact  value  of  the  sum,  the  error  bounds,  possible  outputs  and  their  probabilities. 
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-  QSsimul.ra  - 

function  [outputs,  probabs,  sp,  pmax,  avg,  qubits] =QSsimul (n,  eps,  f) 

# [outputs,  probabs,  sp,  pmax,  avg,  qubits] =  QSiraul(  n,  eps,  f  ) 

# 

# INPUT : 

#n  -  number  of  qubits  for  coding  the  domain  of  the  Boolean  function  f, 

#  i.e.,  the  cardinality  is  2~n 

#eps  -  the  desired  accuracy  0  <  eps  <  1 

#  f  -  Boolean  function  {0,  ..  2~n  -1}  ->  {0,1} 

#0UTPUT : 

#outputs  -  vector  Mxl,  all  possible  outputs  of  the  QS  algorithm 
#probabs  -  vector  Mxl,  corresponding  probabilities 
#sp  -  the  best  estimate 

#pmax  -  the  probability  of  obtaining  an  estimate  with  accuracy  eps 
#avg  -  exact  value 

#qubits  -  the  number  of  qubits  in  the  second  register 
#The  visualization  of  the  final  state  is  in  the  file  out.ps 

qubits=ceil (log2 (pi/ eps) ) ; 

M=2~qubits ; 

N=2~n; 

Reg=ones(M,  N) ; 
outputs=zeros (M, 1) ; 
avg=0 ; 
i=0; 

fo=zeros(N, 1) ; 
for  i=l:N 


f f =f eval (f , i-1)  ; 
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if  ff==l 
avg++ ; 
end 

fo(i)=(-l) ~ff ; 
end 

avg/=N ; 

for  i=2 : M 

sum=(2/N) * (Reg(i-1 , :  )  *f o)  ; 
for  j=l:N 

Reg ( i , j ) =sum-Reg ( i— 1 ,  j)*fo(j); 
end 
end 

Reg=l/ (M*sqrt (N) ) *f f t (Reg) ; 

p=0; 

sp=0; 


probabs=zeros (M, 1) ; 
for  1=1 :M 
for  11=1 :N 

probabs (1) =probabs (1) +abs (Reg (1 , 11) ) ~2 ; 
end 

output s (1)=( sin (pi* (1-1) /M) ) ~2; 
if  (p  <  probabs (1)) 
p  =  probabs (1) ; 
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sp  =  1-1; 
end 
end 

sp=outputs (sp+1) ; 
pmax=0 ; 

for  1=1:M 

if  (abs (outputs (l)-avg)<  eps) 
pmax=pmax+  probabs(l); 
end 
end 

data=[(linspace(0,  M-l,  M) ) ’ , outputs ,  probabs,  avg*ones(M, 1)] ; 
data=  [data,  sp*ones (i , 1)  ,  (avg-eps) *ones (i , 1) ,  (avg+eps)*ones(i, 1)] ; 

gset  output  "outl.eps" 
gset  terminal  postscript  eps 
gset  style  line  1  It  1 
gset  style  line  2  It  2 
gset  style  line  3  It  3 
gset  style  line  4  It  4 
gset  key  outside  bottom 

gset  xlabel  "States  |  j  >  (Total  #  States  =  M) " 

gset  ylabel  "Value  of  the  sum" 

gplot  [0 : M— 1]  [0:1]  data  using  1:2  title  "Possible  outputs"  \ 

with  points, \ 

data  using  1:3  title  "Pr (state  I  j  >)"  \ 
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end 


data  using  1:4  title 


data  using  1:5  title 


data  using  1:6  title 


data  using  1:7  title 


with  boxes  1,\ 
"True  value"  \ 
with  lines  1,\ 
"Estimated  value"  \ 
with  lines  2,\ 
"Error  bounds"  \ 
with  lines  3,\ 

""  with  lines  3 


The  results  of  running  QSsimul  for  three  different  functions  /  are  presented  in 
Figures  12.31 12.41 12.51 


Possible  outputs 
Pr(state  |  j  >) 
True  value 
Estimated  value 
Error  bounds 


Figure  2.3:  QSsimul  output  for  N  =  210,  e  =  0.1,  M  =  32  and  f{k)  —  1  for  k  not 
divisible  by  3  and  f{k)  =  0  otherwise. 
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Possible  outputs 
Pr(state  |  j  >) 
True  value 
Estimated  value 
Error  bounds 


Figure  2.4:  QSsimul  output  for  N  =  210,  £  =  0.1,  M  =  32  and  f(k)  —  1  for  k  not 
divisible  by  8  and  f(k)  =  0  otherwise. 


Possible  outputs 
Pr(state  |  j  >) 
True  value 
Estimated  value 
Error  bounds 


Figure  2.5:  QSsimul  output  for  N  =  210,  £  —  0.1,  M  —  32  and  f{k)  —  1  for  k 
divisible  by  8  and  f(k)  =  0  otherwise. 
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2.4.2  Computation  of  the  final  state  distribution 

Suppose  that  we  want  to  know  the  distribution  of  the  final  state  of  the  QS  algorithm 
without  computing  all  amplitudes  of  the  final  state.  This  problem  can  be  solved 
much  faster  than  the  full  simulation.  First  we  need  to  compute  the  arithmetic  mean 
of  a  Boolean  function  with  cost  of  order  N.  Then,  from  (12. 12ft  we  can  compute 
the  probability  the  final  outcome  |  k).  Thus  the  cost  of  computing  the  final  state 
distribution  is  clearly  of  order  M  +  N,  which  is  of  order  e"1  +  N  for  the  desired 
accuracy  e.  Hence,  we  can  save  a  factor  of  e^1  log2e_1  over  the  cost  of  simulation. 

2.5  Conclusions 

The  results  of  this  chapter  show  the  robustness  of  the  QS  algorithm.  Its  optimality 
in  the  the  worst-probabilistic  error  setting  was  extended  to  two  more  error  settings — 
average-probabilistic  and  worst-average. 

The  average-probabilistic  error  criterion  is  weaker  than  the  worst-probabilistic 
one.  The  difference  is  that  the  probabilistic  error  is  considered  on  the  average  with 
respect  to  the  class  MN  of  input  Boolean  functions.  We  consider  two  measures  on  the 
class  Bjv-  The  first  measure  is  uniform  on  Boolean  functions,  while  the  second  one  is 
uniform  on  arithmetic  means  of  Boolean  functions.  We  show  that  for  the  first  measure, 
the  QS  algorithm  retains  its  optimality  for  a  certain  choice  of  its  parameters.  For  the 
second  measure,  the  average-probabilistic  error  is  essentially  of  the  same  order  as  the 
worst-probabilistic  one,  so  weakening  the  error  criterion  does  not  yield  any  essential 
cost  gain  over  the  worst-probabilistic  error  setting. 

In  the  worst-average  setting  the  QS  algorithm  with  repetitions  is  optimal.  This 
shows  its  superiority  over  classical  randomized  algorithms.  Indeed,  the  worst-average 
error  criterion  with  respect  to  L2  norm  is  analogous  to  the  usual  error  criterion  by 
which  the  error  of  classical  randomized  algorithms  is  considered,  among  them  the 
Monte  Carlo  algorithm.  We  recall  from  [3l]  that  the  complexity  of  the  Boolean  sum- 
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mation  in  the  classical  randomized  setting  is  of  order  e~2,  while  for  the  QS  algorithm 
with  repetitions  we  have  the  optimal  cost  of  order  e^1.  Furthermore,  the  worst- 
average  error  is  a  reasonable  choice  for  algorithms  for  which  we  cannot  compute  the 
actual  a  posteriori  error,  so  the  result  of  an  algorithm  cannot  be  verified.  The  reason 
for  this  is  that  the  worst-average  error  takes  into  account  all  the  outcomes  of  an  algo¬ 
rithm,  while  the  worst-probabilistic  error  deals  only  with  highly  probable  outcomes. 
Finally,  the  worst-average  error  criterion  is  stronger  than  the  worst-probabilistic  one, 
which  can  be  easily  proved  by  using  Chebyshev’s  inequality. 


Chapter  3 
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Chapter  3 

Multivariate  Feynman-Kac  path 
integration 

3.1  Introduction 

In  this  chapter,  we  analyze  the  multivariate  Feynman-Kac  path  integration  problem. 
Although  we  are  mainly  interested  in  the  quantum  setting  we  also  include  the  worst- 
case  deterministic  and  randomized  settings,  which  will  allow  us  to  compare  the  results 
for  all  three  settings. 

We  recall  that  multivariate  Feynman-Kac  path  integrals  are  path  integrals  over 
the  space  of  continuous  functions  from  M+  to  Md,  equipped  with  a  Wiener  measure. 
The  multivariate  Feynman-Kac  formula  is  the  solution  of  the  initial  value  problem 

dz 

—  (u,  t)  —  |Az(u,£)  +  K(u)z(u,  t)  for(u,f)  G  x  (0,  oo),  (3.1) 

z(u,0)=v(u)  (3.2) 

for  the  diffusion  (heat)  equation.  Here  v,V  :  — >  M  are  the  initial  value  function 
and  the  potential  function.  As  usual,  A  denotes  the  Laplacian.  The  solution  z  of  (13.11) 
and  (13.21)  is  given  by  the  famous  Feynman-Kac  formula 
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Here,  C  is  the  set  of  continuous  functions  x  :  M+  — >  M.d  such  that  x(0)  =  0.  The 
path  integral  (13.31)  is  taken  with  respect  to  the  d-dimensional  Wiener  measure  w,  see 
[23,  3D]-  Obviously,  (13. 3 p  only  holds  for  functions  v  and  V  for  which  the  path  integral 
exists.  In  what  follows,  we  assume  that  the  functions  v  and  V  belong  to  a  class  F  for 
which  (13.31)  exists.  This  class  is  precisely  defined  in  Section  [3731 

Various  computational  algorithms,  mostly  randomized,  have  been  developed  for 
the  univariate  case  d  —  1,  where  the  Feynman-Kac  path  integral  is  the  solution  of 
the  diffusion  equation  with  one  space  variable.  A  novel  approach  for  the  univariate 
case  d  —  1  was  proposed  in  [39],  where  a  new  deterministic  algorithm  based  on 
./^-approximation  was  constructed  and  the  complexity  of  the  univariate  case  was 
studied.  Those  results  were  modified  and  generalized  to  the  multivariate  case  in  j2bj 
and  then  improved  in  [29] .  The  multivariate  algorithm  retained  the  structure  of  the 
algorithm  from  [39]  and  is  based  on  uniform  approximation.  We  briefly  discuss  the 
worst-case  setting  in  Sections  13.2.11 13.5.21  and  13.6.11 

As  we  shall  see,  the  quantum  setting  for  multivariate  Feynman-Kac  path  inte¬ 
gration  is  strongly  related  to  the  randomized  setting.  We  present  algorithms  that 
compute  approximations  of  multivariate  Feynman-Kac  path  integrals  in  both  these 
settings.  Both  algorithms  are  based  on  uniform  approximation,  similarly  to  the  de¬ 
terministic  algorithm  from  [29].  We  analyze  the  (informational)  costs  of  these  al¬ 
gorithms,  i.e.,  the  numbers  of  function  evaluations  and/or  quantum  queries  used  to 
compute  approximations  with  a  given  error  bound.  Finally  we  relate  these  costs 
to  the  complexity  of  multivariate  Feynman-Kac  path  integration  in  the  randomized 
and  quantum  settings.  As  in  (29,  39],  the  complexity  is  bounded  from  below  by  the 
complexity  of  multivariate  weighted  integration.  The  upper  bounds  are  provided  by 
the  costs  of  the  algorithms  presented  in  this  chapter.  As  we  shall  see,  the  power  of 
randomization  or  quantum  computation  yields  a  substantial  improvement  over  the 
worst  case  complexity. 

We  now  discuss  the  complexity  results  presented  later  in  this  chapter.  We  know 
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that  deterministic  algorithms  for  the  problem  of  approximating  multivariate  Feynman- 
Kac  path  integrals  with  the  worst  case  error  assurance  are  highly  inefficient  for  large  d, 
see  [29].  This  is  caused  by  the  provable  curse  of  dimensionality  of  this  problem.  More 
precisely,  the  cost  of  computing  an  ^-approximation  by  any  deterministic  algorithm  is 
at  least  of  order  e~d^r.  Here,  r  measures  the  smoothness  of  the  initial  value  and  poten¬ 
tial  functions,  e.g.,  we  consider  the  class  F  of  functions  which  are  r  times  continuously 
differentiable.  Clearly,  the  cost  depends  exponentially  on  d. 

Switching  to  the  randomized  setting  is  one  way  to  vanquish  the  curse  of  dimen¬ 
sionality.  In  this  setting  the  path  integral  is  approximated  by  a  multivariate  integral, 
and  then  this  integral  is  approximately  evaluated  by  a  randomized  algorithm,  e.g.,  by 
the  celebrated  Monte  Carlo  algorithm,  see,  e.g.,  0E3,  also  m  This  yields  a  cost  of 
order  (1/e)2,  so  that  the  dependence  on  d  disappears.  As  we  shall  see,  the  exponent 
of  1/e  can  be  improved  at  the  expense  of  introducing  a  dependence  on  d.  More  pre¬ 
cisely,  for  positive  r  we  construct  an  optimal  algorithm  with  cost  of  order  £~2/(1+2r/d)  ^ 
see  Section  13.6.21 

The  use  of  a  quantum  computer  yields  even  greater  improvement.  Let  us  mea¬ 
sure  the  cost  of  an  algorithm  by  the  number  of  queries  and  function  evaluations  it 
uses.  We  show  in  Section  13.6.31  that  an  optimal  quantum  algorithm  computes  an  e- 
approximation  with  cost  of  order  e_1  with  no  dependence  on  d  and  of  an  optimal  or¬ 
der  £-ldl+r/d')  with  a  dependence  on  d.  Thus,  we  obtain  a  roughly-quadratic  speedup 
over  the  randomized  setting  and  (as  with  the  randomized  algorithm)  an  exponential 
speedup  over  the  worst-case  deterministic  setting. 


3.2  Computational  problem 


We  want  to  compute  an  e- approximation  of  the  path  integral  (13.31)  at  a  given  point 
(u,f)  G  x  [0,  oo)  and  for  arbitrary  functions  v,  V  from  the  class  F.  The  definition 
of  an  ^-approximation  depends  on  the  setting;  this  will  be  made  precise  in  the  next 
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three  subsections.  The  e-approximation  a„ty(u,  t)  is  computed  by  an  algorithm  An 
that  uses  n  function  values  of  v  and  V,  i.e. , 

avy(u,t)  =  An  (u,t,u(ui), . . .  ,v(uk),  l/(ufc+i), . . .  ,h(un)) . 

In  the  quantum  setting,  n  denotes  the  number  of  quantum  queries  and  classical  func¬ 
tion  evaluations. 

3.2.1  Worst-case  deterministic  setting 

In  the  worst  case  setting,  the  error  of  the  algorithm  An  is  defined  as 

eWOT(An)  =  SUp  \zvy{u,t)  -  avy(u,t)\. 
v,veF 

We  also  want  to  determine  the  complexity,  i.e.,  the  minimal  number 

nwor(e,  F)  =  min {n  :  3  An  such  that  eWOT(An)  <  e} 

of  function  values  that  are  needed  to  compute  an  ^-approximation  in  the  worst  case 
setting.  The  worst-case  deterministic  setting  is  analyzed  in  [TE,  25], 

3.2.2  Randomized  setting 

In  this  setting  we  use  randomized  algorithms  and  replace  the  worst  case  error  assur¬ 
ance  by  an  expected  one.  A  randomized  algorithm  An  depends  on  a  random  element  u ; 
chosen  from  some  probability  space  fh  More  precisely,  we  compute 

Gh,v(u,  t- u)  =  AUt u  (u,  t,  «(uWii), . . . ,  u(uW)fc),  V (uWifc+i), . . . ,  V (uUi„J) ,  (3.4) 

with  n  =  E(na;).  This  means  that  we  allow  a  random  choice  of  a  mapping  An^  and 
sample  points  uu^,  as  well  as  the  number  nu  of  sample  points,  whose  expected  value 
is  fixed  and  equal  to  n. 

We  measure  the  randomized  error  of  the  algorithm  An  with  respect  to  the  L2 

erand(An)  =  sup  (M(zvy(u,t)  -  avy(u,t-,uj))2)1/2  . 
v,v&f 


norm,  i.e., 
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To  make  the  notation  more  compact  we  shall  skip  the  index  u  if  the  context  clearly 
indicates  that  we  deal  with  a  randomized  algorithm. 

As  before,  we  want  to  determine  the  complexity,  i.e. ,  the  minimal  expected  number 

nraild(e,  F)  =  min {n  :  3  An  such  that  erand(An)  <  e} 

of  function  values  needed  to  compute  an  ^-approximation  in  the  randomized  setting. 
The  randomized  setting  is  analyzed  in  [27,  28]. 

3.2.3  Quantum  setting 

In  the  quantum  setting  we  use  quantum  algorithms  with  randomized  quantum  queries; 
we  also  assume  that  we  can  perform  function  evaluations  and  arithmetic  operations 
on  a  classical  computer.  These  classical  operations  are  used  to  prepare  an  input 
for  a  quantum  algorithm  and  to  transform  the  outcome  of  a  quantum  algorithm  to 
an  approximation  of  the  exact  solution.  We  will  be  interested  in  minimizing  the 
total  number  of  quantum  queries  and  function  evaluations  needed  to  compute  an  e- 
approximation. 

We  base  our  analysis  on  the  simplified  quantum  model  of  computation  for  con¬ 
tinuous  problems  with  randomized  queries  from  Section  11.1.11  We  refer  the  reader  to 
pEJ  DEI  M,  E21  E2j  for  more  detailed  information.  We  recall  that  for  a  given  class  F  of 
input  functions  /  :  D  — >  C,  we  want  to  approximate  the  solution  operator 

S  :  F^G, 

with  G  being  a  norrned  space  whose  norm  is  denoted  by  ||  ■  ||g-  A  quantum  algorithm 
may  use  a  classical  algorithm  Ps  with  s  classical  function  evaluations  to  transform  a 
given  input  function  /  G  F  into  /  =  Ps(f )  :  D  — >  C,  which  is  then  used  as  an  input 
to  a  quantum  algorithm. 

In  this  chapter  we  only  use  randomized  quantum  queries  as  defined  in  [52],  A 
quantum  algorithm  with  randomized  queries  UnjW(f)  :  TLn  — *■  Hn  is  a  unitary  operator 
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of  the  form 

Un,uj{f)  Qn  Q  f,  uj  Qn—  1  '  '  '  Q  iQ  f.ui  Qoi 

with  unitary  operators  Q0, . . .  ,Qn  and  a  quantum  query  QftW,  for  some  /  G  PS(H). 
The  query  QjtU1  depends  on  a  random  element  u  G  hi.  This  permits  the  computation 
of  approximate  values  of  /  at  randomized  points,  as  explained  in  Section  11.1.11  see 
also  (52]. 

After  performing  the  computation,  we  obtain  a  final  state 

IV’/.w)  =  UnjU(f)  |0)  =  QnQf,uiQn- 1  '  '  '  Q  lQ  f ,  uQ  o\®)  ■ 

We  then  measure  the  final  state  to  obtain  an  outcome  j  G  {0, 1, . . . ,  2k  —  1}  with 
probability 

Pf,u(J)  =  \(^fjj)\2- 

Knowing  the  outcome  j  we  compute  the  final  result  on  a  classical  computer,  and  the 
quantum  algorithm  An  yields 

A  n,w(f,j)  = 

for  some  (j). 

The  error  of  a  quantum  algorithm  with  randomized  quantum  queries  An  is  de¬ 
fined  as 

e^nt(An)  =  sup  (E Eq  || S(f)  -  AntUJ(Ps(f)  J)\\g)1/2  ,  (3-5) 

where  E  is  the  expectation  over  the  probability  space  fl,  and  Eq  is  the  expectation 
with  respect  to  distribution  of  the  quantum  algorithm  outcomes.  We  shall  skip  the 
index  uj  if  it  is  clear  from  the  context  that  we  deal  with  randomized  queries. 

Similarly  to  the  other  settings,  we  want  to  determine  the  (information)  complexity, 
i.e. ,  the  minimal  number 

nquant(£,  F)  =  min {s  +  n  :  3  Ps  3  An  such  that  equant(An,  P8)  <  e} 

of  random  quantum  queries  and  classical  function  evaluations  needed  to  guarantee 
that  the  error  does  not  exceed  e. 
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Remark  3.2.1.  We  briefly  comment  on  the  quantum  error  setting  defined  by  19. <51 ). 
Let  us  concentrate  for  a  moment  on  the  randomness  introduced  by  a  quantum  algo¬ 
rithm,  leaving  aside  randomized  queries.  So  far,  the  literature  dealing  with  continuous 
problems  in  the  quantum  setting  has  mainly  considered  probabilistic  error.  That  is, 
instead  of  taking  an  expectation  with  respect  to  all  possible  outcomes  of  a  quantum 
algorithm,  as  Eq  in  /1 3.  <51).  we  want  an  error  estimate  such  that 


(E||S(/)-/l„,p,(/l,„||J)1/2<e 


to  hold  with  a  given  (high)  probability,  for  any  f  £  F.  Obviously  these  two  ways 
of  measuring  the  error  of  a  quantum  algorithm  are  related.  We  choose  to  study  the 
average  error  for  simplicity  and  two  other  reasons.  The  average  error  is  probably 
more  natural  when  we  cannot  verify  the  result  of  an  algorithm.  Moreover,  113.5 1)  is  a 
stronger  error  criterion  than  that  above. 

The  multivariate  Feynman-Kac  path  integration  problem  in  the  quantum  setting 
is  defined  by  taking  /  =  ( v ,  V),  with  F  x  F  as  the  input  function  class  and  S(f)  = 
zvy(u,  t).  We  shall  make  more  assumptions  on  the  input  function  class  in  the  next 
section. 

3.3  The  function  class  Bp 

To  assure  the  existence  of  the  path  integral  (13.31) .  we  need  to  choose  an  appropriate 
class  of  input  functions  F,  see  also  [29]: 

1.  To  make  the  path  integral  (18.31)  well  defined  we  assume  that  for  every  u  €  M.d, 
the  functional  Lu  :  F  — >  M  defined  by  Luf  =  /( u)  is  continuous,  and  for 
arbitrary  a,  t  £  M+  we  have 


(3.6) 


By  Fernique’s  theorem,  see  e.g.,  [25],  condition  (13.61)  holds  if  there  exists  a  £ 
(0,2)  such  that  ||Tx||,f  =  0(||x||“)  for  ||x||  approaching  infinity,  see  [39]  for 
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details.  Here  and  elsewhere  in  this  chapter,  ||x||  =  denotes  the 

Euclidean  norm  in  Md. 

2.  We  assume  that  F  is  continuously  embedded  into  Loc(Rd).  That  is,  F  C  Loc(Rd ) 
and  there  exists  a  positive  K  such  that 

ll/IU^)  <  K\\f\\F  v/  G  F.  (3.7) 

This  assumption  permits  us  to  relate  the  multivariate  Feynman-Kac  path  inte¬ 
gration  problem  to  uniform  approximation.  By  uniform  approximation,  we 
mean  the  approximation  of  functions  from  F  in  the  norm  of  Loc(Md),  i.e. , 
given  e>0we  want  to  End  a  function  /  e  L00(Md)  such  that 

II  /  /||  Z/oo(Rd)  <£||/||f- 

3.  We  assume  that  we  can  compute  a  uniform  approximation  /  of  the  function  / 
from  the  class  F  by  a  linear  algorithm 

nApp 

/  =  Ci,  u*  e  Ci  e  ^oo(Md)  (3.8) 

i=l 

that  uses  uapp  function  evaluations,  where 

?rApp  =  nApp(£,  F)  =  0(e~a{F) )  as  £  ->  0,  (3.9) 

for  some  positive  a(F).  The  asymptotic  constant  in  (13.91)  may  depend  on  the 
dimension  d.  Usually  the  exponent  cx(F)  depends  on  the  smoothness  and  on 
the  number  of  variables  of  functions  from  F,  see  Section  13.81 

We  stress  that  these  assumptions  are  not  overly-restrictive.  It  is  known  that 
algorithms  of  the  form  (j3.8j)  are  optimal  for  the  uniform  approximation  problem, 
see  [23].  Moreover  the  number  of  function  evaluations  often  depends  on  £  by 
an  expression  similar  to  (13.9)1.  see  also  Section  13781 
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4.  We  restrict  the  norms  of  the  initial  value  and  potential  functions.  Namely,  we 
assume  that  ||v||F  <  (3  and  ||I4||/  <  B  for  given  positive  f3,B.  In  other  words, 
the  pair  (v,  V )  belong  to  the  class 

Bf  =  {(/i,/2)  G  F  x  F  :  \\M\p  <  B ,  ||/2||F  <  /?}.  (3.10) 


3.4  Feynman-Kac  formula  as  a  series  of  multivari¬ 
ate  integrals 


In  this  section  we  briefly  recall  some  results  from  [2S]  that  are  needed  for  our  analysis. 

Without  loss  of  generality  we  can  assume  u  =  0  in  (13.3(1.  Then  we  can  express 
the  path  integral  as  a  series  of  multivariate  integrals 


S(v,  V)  :=  2(0,  t)  -  y  St(y,  V), 


(3,11) 


k= 0 


where 


„  k 

Sk(v,V)  =  /  v(zk+1)T\V(zi)gk(z1,...,zk+l)dz1...dzk+1,  (3.12) 

jR(k+l)d  ^ 


with 


9k(A I;  ■  ■  ■  j  Zfc+i)  /  fk  (tii  •  •  •  j  4>  t,  Zi, . . . ,  zfc+i)  dt i  . . .  c?4  (3.13) 

J  0<ti<— <tfc<t 


and 


fk  (4,  . .  .,tk,t,  zi, . . .  ,zfc+i)  =  ((27r)fc+1fi(t2  -  4)  •  •  •  (t  -  4)) 


-d/2 


,  1  /^llZl||2  ,  II z2  Z1 II 2  ,  ,  ||Zfc+l  -  Zfc||2 

x  exp  --  - — —  +  - - —  H - b  - - — 

2  \  4  4-4  t  -  tk 

Note  that  the  integral  (13.12(1  depends  on  the  input  functions  v  and  V  only  through 


the  product 


hk( zi,  •  •  • ,  zfc+i)  v (zfc+i)  RvM- 


i= 1 


(3.14) 
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Moreover,  the  weight  functions  gk  can  be  computed  in  advance,  albeit  with  difficulty. 
Let  us  recall  also  that 

tk 

IMIli(r(h-i)<*)  =  for  (3‘15) 

and  so  the  norm  of  the  weight  function  gk  decreases  super-exponentially  as  k  goes  to 
infinity. 

3.5  Approximating  one  term  of  the  series 

In  this  section,  we  present  algorithms  approximating  one  term  of  the  series  (13.111). 
To  make  the  notation  more  clear,  we  define  a  weighted  integration  operator 

h(f)  —  /  /(zi, . . . ,  zfc+i)  gk(zu  . . . ,  zfc+1)  dz1...  dzk+1, 

jR(k+l)d 

where  /  :  R(k+1P  M  is  an  integrable  function.  We  can  then  rewrite  one  term  of  the 
series  (13.111)  as 

Sk(v,V)  =Ik(hk). 

In  all  three  settings — worst-case  deterministic,  randomized  and  quantum — we 
shall  use  deterministic  uniform  approximation  of  the  function  hk-  In  the  worst-case 
setting  we  use  uniform  approximation  directly  as  the  main  building  block,  whereas 
in  the  randomized  and  quantum  settings,  we  will  apply  uniform  approximation  as  a 
preprocessing  step  that  will  achieve  variance  reduction. 

3.5.1  Uniform  approximation  by  Smolyak’s  algorithm 

Smolyak’s  algorithm  is  a  powerful  tool  for  approximating  tensor  product  problems. 

k 

/ - A - s 

For  y*,  G  Fk  :=  F  ®  ®  F,  Smolyak’s  algorithm  is  of  the  form 

n(s,k) 

Ue,k(}Pk)  ^  ^  ^Pk(^-i,£, I?  •  •  •  ?  U.z,£:j/c)0,£,A:? 
i= 1 


(3.16) 
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for  some  u ji£)J-  G  and  £i,e,fc  G  L00(M(fc+1)d).  It  is  proved  in  [29,  Lemma  2]  that 


II Tk  -  ^,A:(^fc)llL00(R(fc+1)d)  <  ^ll^fcllFfe, 


(3.17) 


where 


n(e,  k )  <  Co 


fci  +  c2 


In  l/e\ 

k-  1 ) 


(«(F)+l)(fc-l) 

e"“ (F) 

+ 


(3.18) 


for  some  c;  G  M.  Here  a+  denotes  max{a,0},  and  the  right  hand  side  of  (13.18))  is 
dehned  to  be  CoS-"*^  when  k  —  1.  We  shall  use  Smolyak’s  algorithm  to  approximate 
the  functions  hk  dehned  by  (13.14)). 


3.5.2  Deterministic  algorithm 

We  first  consider  approximating  Sk  (v ,  V )  by  the  algorithm 

$et(v,V)  =  Ik(hk,e), 

where  hkyE  =  UEjkhk  and  hk  is  dehned  by  (13.14)).  Using  (I3.16p.  we  can  rewrite  (pfet  as 

n(e,fc+ 1) 

0£et(u  V)  =  hfe(uj)£ii, . . . ,  Uj!£,fe+i)4(Ci,£,fc+i).  (3.19) 

We  stress  that  4(Ci,e,fc+i)  does  not  depend  on  the  input  functions  v  and  V  and  can 
be  precomputed. 

We  note  that  the  error  of  the  algorithm  <^et  satishes 

\Sk(v,V)-^\v,V)\<e^A-,  (3.20) 

Indeed  by  (13. 10)1 .  (13.15)1  and  (13.171)  we  get 
I Sk(v,  V)  -  </>£et(w,  V)\  <  Ik(  | hk  -  UE,k+1hk | ) 

B 1 3 ^ 

<  II hk  —  t4)fe+ihfe||Loo(R(fe+i)d)  ||5,fc||L1(R(fc+i)d)  <  £  — ^ — ; 

moreover,  the  total  number  of  function  evaluations  used  by  the  algorithm  (p^et(v,  V) 
is  n(e,  k  +  1),  as  dehned  in  (13.18)). 
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3.5.3  Variance  reduction 

The  idea  underlying  variance  reduction  is  as  follows.  First  we  compute  hk,£  = 
U£jk+i{hk)  using  n(e,  k  +  1)  function  values.  Then  we  compute 

n(e,fc+l) 

hk{hk,e)  ^  hk(Uij£j ij  -  -  -  ,  ^i,e,k+l)^k(Ci,£,k+l) • 

i=  1 

Observe  that  the  functions  Ci,e,k+ l  do  not  depend  on  the  input  functions  v  and  V,  and 
so  the  integrals  Ik(Ci,e,k+i)  can  be  precomputed.  We  also  stress  that  hk.e  and  Ik(hk,e ) 
are  deterministic.  We  shall  use  randomized  or  quantum  algorithms  to  approximate 
the  multivariate  integrals 

hk,e)- 

Since  the  error  depends  on  the  norm  ||  hk  —  hk,e  ||.L00(iR(fc+1)d)>  which  is  now  small,  we 
can  do  this  efficiently.  We  present  the  details  in  the  following  two  sections. 

3.5.4  Randomized  algorithm 

To  make  our  formulas  shorter,  we  define 

fk,e  -=  hk  hk:e- 

We  use  a  randomized  algorithm  of  the  form 


/rand / 


v)  =  h(hkf) + 


(3.21) 


Here 


QTd(f) 


l 

m 


3= 1 


(3.22) 


denotes  the  classical  Monte  Carlo  algorithm  with  m  randomized  sample  points.  Ran¬ 
domized  sample  points  are  chosen  with  respect  to  the  density  gk/\\gk\\L1(s.(k+1)d)- 
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Using  the  well-known  error  formula  for  the  classical  Monte  Carlo  algorithm,  we 
conclude  that 


(E  (4(40  -  r))2)1/2  =  (E  (4(4,0  -  «"d(A0)2) 

1 


2\V2 


(Var  (/M))1/2,  (3.23) 


with 

Var(A,«)  =  4(/i2J-(4(A,e))2. 
Clearly,  from  (I3.17P  and  then  from  (I3.10p  and  (13 .151).  we  get 

tk 
k\ 


Ur  ,1  M  „  /  IHIfIIVII^^  B(3k  tk 

(Var (fk}£))  <  —  ||/fc>£||Loo(R(fe+i)d)  <  e - rn - <  £ 


k\ 


k\ 


This  yields  the  error  estimate 


(e(4(40  -  ^'“d(«,r))2)1/2  < 

v  7  y/m  k\ 

and  the  total  number  of  function  evaluations  being 


(3.24) 


(3.25) 


n(e,  k  +  1)  +  m. 


(3.26) 


3.5.5  Quantum  algorithm 

The  structure  of  the  quantum  algorithm  is  similar  to  the  randomized  one  and  has  the 
form 

v)  =  +  «""‘(A0,  (3.27) 

with,  as  before,  fk^  =  hk  —  hk,e-  Here,  we  use  a  quantum  algorithm  Q ^nt,  with 
k  randomized  quantum  queries,  that  approximates  the  classical  Monte  Carlo  algo¬ 
rithm  (13.221) .  In  [IH],  the  problem  of  approximating 

1  m 

f(xi ) 

3=  1 

was  analyzed  for  Boolean  functions  /.  By  reducing  the  summation  problem  for 
bounded  real  functions  to  the  summation  problem  for  Boolean  functions  as  in  na, 
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we  see  that  a  result  similar  to  that  of  ra  holds.  From  [19]  and  (13. 1 7j)  we  conclude 
that 

m  2\  1 /2 

M  ^  «T(/m)) )  =o(hiA,«ii1„(I».+,M))  =o(Ni3> 


3  = 1 


By  integrating  over  randomized  sample  points,  we  obtain 


quant  (  f  M  2^  1/2  ^  (e  Bf3kt 


EEq| ) 


=  0  - 


k  -tk 


(3.28) 


k  k\ 

The  total  number  of  randomized  quantum  queries  and  function  evaluations  is 

n(e,  k  +  1)  +  k. 

We  stress  that  this  number  does  not  depend  on  m,  which  is  only  used  for  the  definition 
of  the  Monte  Carlo  algorithm  providing  sample  points  for  the  quantum  algorithm. 
We  now  estimate  the  total  error  as 


(EE „(/*(/»*)  -  «>—(*>. V)f)'n  <  (E  -  QT\hf)Y) 


A  1/2 


+  (EEq|  QTd(fk,e)-Q%T(f^)n 


A  V2 


This,  by  (|3.25j)  and  (13.281).  yields 

(EE,(4(At)  -  <#.?“>,  V)ffn  =  O 

Letting  m  =  k2  we  get  the  error  bound 


£  B(3k  tk  e  B/3k  tk 


m  k\ 


k  k\ 


(EE,(4(fe)  -  <P^(v,V))2)1/2  =0(1 


£  2  Bf3ktk 


n  k\ 


(3.29) 


using 


n(e,  k  +  1)  +  k  (3.30) 

function  values  and  quantum  queries.  For  the  sake  of  convenience  we  denote 


a  quant  i quant 

T£,K  r£,m,/C 


with  m  =  n2 . 
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3.6  Complete  algorithms 


Based  on  the  previous  two  sections,  we  are  ready  to  present  algorithms  computing 
an  ^-approximation  of  multivariate  Feynman-Kac  path  integral  S(v,  V).  We  simply 
approximate  consecutive  terms  of  the  series 

OO 

S(v,V)  =  Y.S^v’V'l 

k= 0 

by  the  algorithms 


or  0™nd 


/quant 

Or  0  quant  • 

£k 


Here,  the  accuracies  £det,  £™nd  and  e)(uaut  in  the  corresponding  settings  are 

k\ 


4et  =  ^ 


B/3k  tk2k+l  ’ 


rand  _  2/{a{F)+2) 

6k  ~  £ 


quant  _  l/(a{F)+l) 

£k  ~  £ 


k\ 


B/3k  tk2k+l  ’ 
k\ 


(3.31) 

(3.32) 

(3.33) 


Bpk  tk 2k+2  ’ 

with  the  numbers  of  randomized  sample  points  rrik  and  quantum  queries  nk  being 


mk  =  \s-MF)/HO+ 2)1  ;  Kk  =  ^£-a(F)/(a(F)+ 1)-| 


(3.34) 


Since  mk  and  Kk  are  independent  of  k,  we  shall  drop  the  indices  and  write  m  =  mk 
and  n  =  nk. 

Remark  3.6.1.  We  see  from  the  definitions  \3.31\).  \3.32fi  and  \3. 33)  that  for  k  going 
to  infinity,  we  have  £det,  erknd  and  £^uant  also  tending  to  infinity  super- exponentially. 
Then,  by  ([ 3 .  7|)  and  then  by  l\3.20\).  1)3. 251 )  and  \3. 291).  we  see  that  for  k  =  0(lne-1) 
the  deterministic  zero  algorithm  provides  sufficient  accuracy  in  all  three  settings.  See 
also  the  following  three  subsections. 

We  are  now  ready  to  present  the  final  forms  of  the  algorithms  and  analyze  their 


errors  and  costs. 
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3.6.1  Deterministic  algorithm 

The  deterministic  algorithm  is  of  the  form 

jVdet 

offer)  =  n 

k= 0 

with  the  finite  sum  limit  N^et.  From  Remark  13.6.11  we  see  that  for  £^et  >  Kk+1,  the 
zero  algorithm  yields  a  sufficient  accuracy  so  we  can  define 


N^et  =  min {k  e  N  :  ^et  >  Kk+1} 


and  N^et  =  0(ln£  x)  .  The  error  of  the  algorithm  <fUet  can  be  estimated  from  (13.201) 
and  (I3.3ip  by 

OO 

| S(v,  V )  -  <£>*>,  V)\  <  I h(hk)  ~  V)\<e.  (3.35) 

k= 1 

The  number  n(<f>^et)  of  function  evaluations  used  by  <fUet  satisfies 


n 


(<f>£)  =  0(e-a{F)~s )  V<5  >  0. 


(3.36) 


The  proof  is  given  in  [29]  and  is  similar  to  the  one  from  [39].  First,  we  have 

OO 

i(*f )  <  +  a 


n( 


k= 0 


Using  (I3.18P  and  (13.311),  we  obtain  the  bound 

111  1  /  E^V  \  ^k^k^k+l  \  aH) 


»(«f)=0  1  +  E  '.+ 


k= 1 


c2- 


k  )  ,  V  k\ 

and  it  can  be  shown,  similarly  to  [39],  that 

In  1  /  E^e^  \  (aU)+i)k  /  pk^2k+i  \  aV) 


OO 

(  0+0- 

k= 1 


-  +  (t)  =°(^ 


for  all  8  >  0. 
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3.6.2  Randomized  algorithm 


The  randomized  algorithm  approximating  S(v,  V )  is  of  the  form 


Ni‘ 


&r>,v)  =  5] 


i  rand 


k= 0 


with  jV5and  defined  as 


iVerand  =  min  {fceN:  4and  >  mKk+1}. 

Again  A£rand  =  0(ln£_1)  .  As  in  Section  [3.6.11  by  (13.25(1  and  (I3.321).  we  can  prove 
that  the  error  of  the  algorithm  <f>£and  satisfies 

1  /2 

(E (S(v,  r)  -  V^))2) 1/2  <  £  (E(4(ht)  -  V)f )  <  £  (3.37) 

fc=l 

The  number  n(<f>£and)  of  function  evaluations  used  by  the  algorithm  <f>£and  satisfies 

n  (<4and)  =  O  (£-MO/0(f)+2)-6^  ys  >  Q  (3.38) 


The  proof  is  based  on  an  argument  similar  to  that  from  Section  13.6.11  By  the  bounds 
(13. 18j) .  (13.261)  and  (13.34jl  we  can  estimate  n(<f>£and)  by 


n  ($™nd)  =0  ( 1  +  Ba(F)  +  J2  (  Cl  + 


k= 1 


c2- 


lnl/4andX  {a(F)+1)k 


X 


pkj-k^k+l  \  aT) 


k\ 


T~2a(F)/(a(F)+2) 


Similarly  to  Section  13.6.11  we  can  show  that 

111  1  /erand  \  (a(F)+Lk  /  Qk^k^k+X  \  aH) 


k= 1 


C2- 


=  0(e-s), 


for  all  8  >  0,  which  proves  (13.38]). 
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3.6.3  Quantum  algorithm 


The  quantum  algorithm  <f>quant  is  defined  as 

^y-quant 

z '  £k  >K 

k= 0 

with  A'quant  satisfying 

^quant  =  min  £  N  .  £rand  >  ^  ^k+l 

and  ./Vquant  =  0(ln  e_1)  .  As  before,  we  can  easily  prove  that  the  error  of  the  algorithm 
is  of  order  e.  Indeed,  by  f!3.29p  and  (13.331)  we  get 

oo  1  /2 

(EE,(5(b,  V)  -  »?"“*(«,  V)ffn  <  E  (EE,U*(ft*)  -  >,  I"))2)  '  =  0(e). 

fc=0 

The  number  of  function  evaluations  and  quantum  queries  rz  (<t>^fuant)  of  the  algorithm 
cp quant  can  pe  estimated  as 

n  (<f>f  ant)  =  O  r)-5)  v5  >  0  (3.39) 


As  in  Section  13.6.21  the  proof  follows  from  the  estimate 


n  (<f>f  ant)  =  O  (  Ba W  +  J2  (  ci  + 


k= 1 


C2 


In  1  y,£<iuan*  \  (°f-0+i  A 


x 


(td  ak+kryk+l  \  a(F) 

1  *  j  +l)  £-a(F)/{a(F)+ 1) 


and  the  fact  that 


k= 1 


Em  + 


c2- 


ln  1  /£quant  \  ("(F)+1)fc  /  Q/]kik2k+1  \ 


k 


PAM  =  °«o 


for  all  5  >  0. 

Obviously  we  can  obtain  the  error  estimate 


(EEq(5(u,  O)  -  $fant(u,  0))2)1/2  <  £ 


(3.40) 


by  redefining  e  modulo  a  factor  that  would  increase  the  asymptotic  constant  in  the 
estimate  of  the  number  of  function  evaluations  and  quantum  queries  (13. 381). 
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3.7  Complexity  of  multivariate  Feynman-Kac  path 
integration 

The  analysis  of  the  complexity  of  the  multivariate  Feynman-Kac  path  integration  in 
randomized  and  quantum  settings  is  based  on  that  presented  in  |29]  and  [32] . 

3.7.1  Lower  bounds 

Lower  bounds  for  our  problem  complexities  are  provided  by  the  complexities  of  mul¬ 
tivariate  weighted  integration  problem.  By  this  problem,  we  mean  an  approximation 
of  the  integration  operator  /  :  F  — ►  M  defined  by 

1(f)  =  (2vr t)~d/'2  I  /( u)  exp(— ||u||/(2t))  du  V  /  G  BF. 

JRd 


1.  Let  ddet  be  a  deterministic  algorithm  that  uses  n  function  values  and  approx¬ 
imates  the  integration  operator  /.  We  say  that  this  algorithm  computes  an 
e- approximation  of  the  weighted  integral  if 

\I(f)-A^(f)\<e  V/GF. 

We  denote  by  nf^ ^(e,  BF)  the  minimal  number  of  function  values  needed  to  com¬ 
pute  an  e- approximation  in  the  worst-case  deterministic  setting  for  the  class  BF. 

2.  Let  d^and  be  a  randomized  algorithm  that  uses  n  function  values  and  approxi¬ 
mates  the  integration  operator  I.  We  say  that  this  algorithm  computes  an  e- 
approximation  of  the  weighted  integral  if 

(E  (/(/)-  AZT'U))  )  <e  V/eSf. 


We  denote  by  nj^(e,  BF)  the  minimal  number  of  function  values  needed  to  com¬ 
pute  an 

e- approximation  in  the  randomized  setting. 
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3.  Let  kLquant  be  a  quantum  algorithm  that  uses  n  randomized  quantum  queries  and 
approximates  the  operator  I.  We  say  that  kLquant  computes  an  e- approximation 
of  the  weighted  integral  if 

/  o\  1/2 

(EE,(/(/)-^rnt(/))  )  <£  V/eSF.  (3.41) 

We  define  Si?)  as  the  minimal  number  of  quantum  queries  needed  to 

compute  an  ^-approximation. 

Since  £(1;,  0)  =  I(v),  as  in  [25],  we  can  reduce  multivariate  Feynman-Kac  path 
integration  to  multivariate  integration  with  Gaussian  weight  by  taking  V  =  0.  More¬ 
over,  (13.1  Op  and  (13. 7p  imply  that 

n?NT(e,Si?)  <  ndet(e,BF), 
n™nd(£,BF)<nrand(£,BF), 
n^\e,BF)<n^\£,BF). 

3.7.2  Upper  bounds 

Obvious  estimates  of  the  complexity  of  the  multivariate  Feynman-Kac  path  integra¬ 
tion  are  provided  by  the  (information)  cost  of  the  algorithms  derived  in  Section  13.61 
Thus,  by  (13. 361).  (13.381)  and  (I3.39p  we  get 

ndet(e,  BF)  =  O  (e-«(^))-«)  , 
nrand(e,  BF)  =  O  (£-WWWj  , 
nquant(e,SF)  =  O  (e-db/btb+GJ) 

for  all  5  >  0,  where  01(F)  is  the  exponent  of  the  uniform  approximation  problem 
complexity  for  the  space  F  containing  the  class  BF  defined  by  (j3.9j) . 

From  the  previous  two  sections,  we  can  determine  when  the  deterministic,  ran¬ 
domized  and  quantum  algorithms  presented  here  are  almost  optimal.  This  happens 
for  the  classes  of  input  functions  for  which  worst-case  deterministic,  randomized  and 
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quantum  complexities  of  the  integration  problem  defined  in  Sectionl3.7.1lare  of  orders 
£~2a(F)/(AF)+2)  and  £~a(F')/(a(F)+1')  respectively.  We  present  two  examples  of 
such  classes  in  the  next  section. 

3.8  Examples 

In  this  section  we  present  two  examples  of  function  classes  F  satisfying  the  assump¬ 
tions  from  Section  13.31  and  compute  lower  and  upper  bounds  of  of  the  complexities  of 
the  multivariate  Feynman-Kac  path  integration.  For  both  examples,  the  algorithms 
presented  in  this  chapter  are  almost  optimal. 

3.8.1  The  Sobolev  space  of  compactly  supported  functions 

Let  F  be  a  class  of  d  variate  r  times  continuously  differentiable  functions  whose 
supports  are  contained  in  a  cube  [ a,b}d  C  Md.  Thus  F  is  a  subclass  of  the  Sobolev 
space  W^d([a,b}d)  with  the  norm 


where  a.  =  [ay, . . . ,  ari]  G  and  /*“)  =  /dai  ■  ■  ■  dad. 

Clearly,  assumptions  1  and  2  of  Section  13.31  are  satisfied.  From  00E3],  we  also 
know  an  optimal  uniform  approximation  algorithm  that  satisfies  assumption  3  with 
the  exponent  ct(F)  =  d/r.  Based  on  this  algorithm,  we  can  construct  almost  optimal 
algorithms  and  <F^uant  computing  e- approximations  of  the  multivariate 

Feynman-Kac  path  integral  (13. 3p  for  (v,V)  G  Bp  in  the  sense  of  (13. 35f) .  (13.371)  and 
(13.401).  with  the  number  of  function  evaluations  and/or  quantum  queries  being  roughly 
of  order  £~d/r,  £-2/(1+2r/d)  and  £~1/(1+r/c0.  y\fe  see  that  the  cost  of  the  deterministic 
algorithm  depends  exponentially  on  the  dimension  d,  whereas  for  the  randomized  and 
quantum  algorithms  the  exponents  of  £-1  are  at  most  2  and  1  respectively.  Thus,  the 
curse  of  dimensionality  present  in  the  worst-case  deterministic  setting  is  vanquished 
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in  both  the  randomized  and  quantum  settings.  This  corresponds  to  the  case  when 
the  exponent  of  e_1  is  2  in  the  randomized  setting  and  1  in  the  quantum  setting. 
We  can  decrease  these  exponents  to  2/(1  +  2 r/d)  and  1/(1  +  r/d ),  at  the  expense  of 
introducing  a  dependence  on  the  number  d  of  variables. 


3.8.2  Periodic  functions 

Although  this  example  was  considered  in  [29],  we  repeat  the  details  for  the  reader’s 
convenience.  Following  pf5j  we  consider  the  class  T  of  27r-periodic  functions  /  : 
[0,  2'a]d  — >  satisfying  the  condition 


V/gT  Vj  =  l,...,d  3ipj  G  L00([-27r,27r]‘i) 

1  f2n 

/(x)  =  —  J  (pj(x1,...,xj-t,...,xd)Fr(t)dt,  (3.42) 


where  r  >  0  and 

00 

Fr(t)  =  1  +  2  k~r  cos  [k  t - —  j  . 

k= o 

The  norm  in  the  class  T  is  defined  as 


ll/ll 


$  — 


1 

d 


d 

^  1  \\Tj  llicx)([—  27T,27r]d)) 

j= 1 


where  the  c pj  are  functions  from  the  representation  f)3.42j)  of  the  function  /.  The 
paper  [45]  provides  a  linear  algorithm  UE  that  computes  a  uniform  e-approxinration 
of  functions  from  the  class  T,  so  that 


ll/-^/IU«([o)2W]-)<e||/lk  V/GT, 

with  a  cost  of  order  e~d^r . 

Let  F  denote  the  class  of  functions  /  :  M  that  are  periodic  extensions  of 

functions  from  T.  Let  ||/||f  :=  ||/|[o,27r]dlk-  The  problem  of  uniform  approximation 
for  the  class  F  can  be  obviously  solved  using  the  algorithm  mentioned  above,  with 
the  same  cost  as  for  the  class  T.  Similarly  to  the  previous  example,  we  have  to 
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check  the  assumptions  of  Section  13.31  It  is  easy  to  see  that  for  /  6  F,  z  e  Md,  and 
arbitrary  j  E  { 1,2, ... ,  d},  we  have 

/(z)  <  ||  /||  L°°(Md)  —  \\f\  [0,27r]d  II  L°° ([0,27r]d)  A  C  ||  <fij  ||  L°° ([-2-k ,2ix\d)  i 

with  C  =  (27t)^1  \Fr(t)  \  dt.  Hence 

I/(z)|<||/||l-Cr-)<C7||/||f 

and  so  function  evaluations  are  continuous.  The  remaining  assumptions  follow  im¬ 
mediately. 

Finally  we  can  construct  the  algorithms  <h^et,  <|uand  and  ^quant  ag  -n  ^jie  previous 
example.  Since  the  uniform  approximation  exponent  a(F)  is  the  same  for  both  ex¬ 
amples  the  algorithms  approximating  multivariate  Feynman-Kac  path  integrals  for 
the  space  of  periodic  functions  have  properties  similar  to  the  algorithms  from  the 
previous  example. 

3.9  Conclusions 

The  results  of  this  chapter  provide  estimates  of  the  complexity  of  multivariate  Feynman- 
Kac  path  integration  in  the  worst-case  deterministic,  randomized  and  quantum  set¬ 
tings.  We  also  present  optimal  algorithms  in  these  three  settings.  Although  this 
problem  considerably  differs  from  ordinary  multivariate  integration,  the  complexity 
estimates  that  we  obtain  are  similar  to  those  for  multivariate  integration. 

In  Section  11.2.2.21  we  indicated  that  one  of  our  goals  is  to  compare  the  complex¬ 
ity  of  multivariate  Feynman-Kac  path  integration  to  the  complexity  of  general  path 
integration  considered  in  |48j.  This  paper  deals  with  path  integrals  over  a  separable 
Banach  space  A",  with  respect  to  a  zero  mean  Gaussian  measure  /i,  with  eigenvalues 
A i  =  0(Fk),  k  >  1.  For  a  functional  /  :  A"  — >•  M  being  s  times  Frechet  differentiable, 
the  general  path  integral  is  defined  as 

S(f)=  [  f(x)fj,(dx). 

J  x 


(3.43) 
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We  stress  that  in  [48]  evaluations  of  the  functional  /  from  (13 ,43ft  are  used,  while  for 
multivariate  Feynman-Kac  path  integration  we  do  not  assume  that  we  can  compute 
the  values  of  the  integrand  in  (13. 3p .  Instead  we  can  only  compute  the  values  of  the 
input  functions — initial  condition  and  potential.  We  compare  the  results  from  [48] 
with  the  complexities  of  multivariate  Feynman-Kac  path  integration  for  the  input 
function  class  considered  in  Section  13.8.11 

The  cost  of  the  worst-case  deterministic  algorithm  of  [48]  is  of  order  e~£  7(s)/fc;  with 
p(s)  =  1  +  8S!i.  This  estimate  is  significantly  larger  than  the  complexity  order  e~d!r 
for  multivariate  Feynman-Kac  path  integration.  Since  the  deterministic  algorithm 
presented  in  this  chapter  is  optimal,  its  cost  order  is  also  e~d!r  and  we  see  that  in  the 
deterministic  setting  the  special  structure  of  the  Feynman-Kac  path  integral  provides 
a  significant  improvement. 

The  randomized  algorithm  for  general  path  integration  of  [48]  has  the  cost  order 
£-2-7 (s)/fc,  wiiereas  tiie  complexity  of  the  multivariate  Feynman-Kac  path  integra¬ 
tion  (with  input  functions  as  in  Section  13.8ft  in  the  randomized  setting  is  of  order 
2/(1+2 r/d)'  "phis  js  aiso  the  order  of  the  optimal  randomized  algorithm  presented  in 
this  chapter.  Thus,  our  algorithm  has  a  polynomial  gain  for  moderate  r/d.  The  same 
can  be  said  about  the  quantum  algorithm  of  [48],  whose  cost  order  is  e^1,  while  our 
algorithm  has  the  cost  of  order  e~1T1+rld\  Thus  for  moderate  d/r  the  special  struc¬ 
ture  of  the  Feynman-Kac  path  integral  yields  an  improvement  in  both  the  randomized 
and  quantum  settings. 
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Chapter  4 

Open  problems 


The  analysis  of  Boolean  summation  and  multivariate  Feynman-Kac  path  integration 
from  Chapters  [2]  and  [3]  leaves  some  open  problems  and  possible  extensions.  In  this 
chapter  we  present  three  research  problems  that  seem  to  be  interesting  and  challeng¬ 
ing. 


4.1  Average- average  error  for  quantum  Boolean 

summation 

The  quantum  Boolean  summation  was  analyzed  in  Chapter  [2j  We  recall  that  the  per¬ 
formance  of  the  QS  algorithm  is  measured  with  respect  to  the  input  Boolean  function 
class  Byv  and  the  final  state  distribution.  This  corresponds  to  the  four  settings — worst- 
probabilistic,  average-probabilistic,  worst-average  and  average- average.  We  studied 
the  worst  performance  with  respect  to  the  class  Mjy  and  the  probabilistic  and  average 
performance  with  respect  to  the  final  state  distribution.  We  also  studied  the  aver¬ 
age  performance  with  respect  to  the  class  IB  at  and  the  probabilistic  performance  with 
respect  to  the  final  state  distribution. 

What  remains  to  be  studied  is  the  average-average  error,  which  measures  the 
average  performance  with  respect  to  both  a  distribution  p  on  the  Boolean  function 
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class  Btv  and  the  distribution  of  the  QS  algorithm  final  state.  We  define  the  average- 
average  error  as  follows 


eavg-avg(Qg)  _  £  ^P(/)P/(j)|S(/)-QS(/,j)|« 


,avg-avg 


As  with  the  worst-average  error  we  measure  the  average  performance  in  the  Lq  norm, 
we  allow  q  G  [l,oo].  Similarly  to  the  average-probabilistic  error  we  may  consider 
two  measures  on  the  set  Bjy:  pi — the  uniform  distribution  over  Boolean  functions, 
and  P2 — the  uniform  distribution  over  the  means.  We  stress  that  the  average-average 
error  setting  for  the  measure  p2  and  q  =  1  provides  a  lower  bound  for  the  complexity 
of  quantum  summation  and  integration  when  randomized  queries  are  used.  Lower 
bounds  for  the  general  quantum  Boolean  summation  problem  were  established  in  [36] . 
We  expect  that  for  uniformly  distributed  Boolean  functions  (the  measure  pi),  we  shall 
observe  that  the  cost  of  the  QS  algorithm  will  on  M,  as  in  the  worst-probabilistic 
setting.  For  the  uniformly  distributed  means  (the  measure  p2)  the  results  are  probably 
quite  different. 

4.2  Multivariate  Feynman-Kac  path  integration  in 
the  quantum  setting  with  deterministic  queries 

We  used  the  concept  of  randomized  query  developed  in  [52J  to  study  the  multivariate 
Feynman-Kac  path  integration  in  the  quantum  setting  in  Chapter  [3]  It  would  be 
natural  to  consider  deterministic  (bit)  queries  used  in  studies  of  almost  all  continuous 
problems  mentioned  in  Section  11.11  However,  as  we  mentioned  before,  the  analysis 
of  the  multivariate  Feynman-Kac  path  integration  in  the  quantum  setting  leads  to 
some  technical  difficulties.  That  is  why  our  results  has  been  so  far  established  only 
for  randomized  queries.  It  would  be  important  to  determine  whether  the  power  of 
deterministic  queries  is  comparable  to  that  of  randomized  queries. 
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4.3  Multivariate  Feynman-Kac  path  integration  in 
finite-order  weighted  spaces 

The  notion  of  finite-order  weighted  spaces,  see  e.g.,  [331,  3S  ED],  has  recently  been 
studied  for  multivariate  problems.  Finite-order  weights  are  used  to  model  continuous 
multivariate  problems  for  which  d-variate  functions  can  be  decomposed  as  sums  of 
functions  of  fewer  variables.  The  weights  describe  the  relative  importance  of  each 
group  of  variables.  Multivariate  problems  over  finite-order  weighted  spaces  are  often 
tractable.  That  is,  the  minimal  number  of  function  values  needed  to  compute  an 
e-approximation  is  polynomial  in  e_1  and  d.  The  finite-order  weighted  structure  of  a 
function  space  seems  to  also  be  very  promising  for  multivariate  Feynman-Kac  path 
integration.  The  integrands  that  are  used  in  applications  often  have  this  kind  of 
structure.  Moreover,  uniform  approximation,  which  is  used  in  optimal  algorithms  for 
all  three  settings  discussed  in  Chapter  [3l  is  tractable  for  finite-order  weighted  spaces, 
see  (32]  •  Therefore,  it  would  be  important  to  also  study  multivariate  Feynman-Kac 
for  finite-order  weighted  spaces  in  the  quantum  model  of  computation. 
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